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I*  INTRODUCTION 


This  report  summarizes  the  results  of  a  two-year  investigation  of 
optical  waveguide  scattering  phenomena.  It  builds  on»  and  in  some  cases  cor¬ 
rects,  information  obtained  during  the  first-year  investigation,  described  in 

Ref.  1,  The  program  was  funded  by  the  Air  Force  in  connection  with  their 

(2) 

on-going  effort  to  develop  an  integrated  optical  spectrum  analyzer.  Air 

Force  personnel  have  long  recognized  that  the  signal  processing  capabilities 
of  an  integrated  optical  spectrum  analyzer,  as  measured  by  its  dynamic  range, 
was  likely  to  be  limited  by  scattering  events  in  the  waveguide  carrier. 

The  objective  of  the  present  program  has  been  to  study  these  events,  obtain  in¬ 
formation  regarding  the  sources  of  scattering,  and  develop  fabrication  proce¬ 
dures  that  minimize  the  number  or  effectiveness  of  these  sources. 

Owing  to  the  difficulty  of  the  problem,  only  one  of  several  candidate 
waveguides  has  been  considered.  That  has  been  the  Ti-dif fused  LiNbO^  waveguide, 
currently  the  waveguide  of  choice  for  the  spectrum  analyzer  application.  At 
the  conclusion  of  the  research,  we  are  able  to  state  with  confidence  that  this 
class  of  waveguide  is  indeed  useful  for  spectrum  analyzers  having  a  dynamic 
range  of  40  dB,  as  required  by  currently  envisioned  device  applications. 

Unfortunately,  this  conclusion  cannot  be  based  on  the  actual  fabrica¬ 
tion  and  testing  of  a  waveguide  commensurate  with  this  level  of  performance. 
Spurious  scattering  associated  with  prism  input  coupling  has  constituted  a 
source  of  noise  that  has  prevented  this.  We  have  developed  methods  for  dealing 
with  this  noise  and  have  derived  an  encompassing  if  not  rigorous  theory  to  use 
in  interpreting  our  results .  This  theory  allows  us  to  reach  conclusions  regard¬ 
ing  the  number  and  effectiveness  of  scattering  centers  in  our  waveguides.  From 
this  we  can  predict  the  performance  of  waveguides  fabricated  using  altered  dif¬ 
fusion  conditions.  Our  opinion  that  LiNbO^  is  an  eminently  suitable  substrate 
for  an  integrated  optical  spectrum  analyzer  is  based  on  the  highly  favorable 
outcome  of  these  predictions. 

Progress  toward  this  conclusion  has  resulted  from  achievements  fully 
described  in  the  body  and  appendices  of  this  report.  Section  II  contains  a 
listing  and  discussion  of  the  physical  basis  for  nine  potential  sources  of 
scattering  relevant  to  LiNb03  waveguides.  They  range  from  those  associated 
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with  substrate  preparation,  such  as  polishing,  to  those  associated  with  dif~ 
fusion  phenomena,  such  as  lithium-titanate  compound  formation  at  the  diffused 
waveguide  surface. 

Section  III  is  a  compendium  of  experimental  methods  that  may  be  used  to 
study  the  scattering  centers  itemized  in  Sec.  IT,  These,  include  diagnostic 
methods  which  ideally  produce  direct  evidence  of  individual  scattering  centers, 
as  well  as  optical  waveguide  experiments  from  which  a  picture  of  the  scattering 
centers  is  pieced  together  using  indirect  evidence.  We  have  included  in  the 
discussion  our  experimental  results  for  those  methods  which  were  only  used  in 
the  program  briefly.  This  included  virtually  all  diagnostic  tools,  which  were 
^rankly  not  very  useful  in  providing  information  about  the  elusive  scattering 
centers  in  Ti  diffused  LiNbO^. 

In  Sec.  TV,  we  present  the  theoretical  analysis  necessary  for  the  inter¬ 
pretation  of  our  optical  waveguide  experiments.  The  section  begins  with  the 
derivation  of  a  formula  for  calculating  spectrum-analyzer  dynamic  range  in  terms 
of  scattering  cross  section,  scattering-center  density,  and  type  of  scattering 
involved  (Rayleigh,  Rayleigh-Gans-Debye,  or  Mie) .  Then  we  move  to  the  deriva¬ 
tion  of  formulae  for  scattering  cross  sections.  By  inserting  in  these  formulas 
the  values  of  the  scattering  parameters  estimated  for  the  various  mechanisms  of 
Sec.  II,  we  are  able  to  predict  the  spectrum -analyzer  dynamic  range  associated 
with  each  mechanism.  The  results  of  this  section  show  that  in-plane  scattering 
caused  by  surface  compounds  formed  during  diffusion  is  likely  to  be  an  important 
limiter  of  spectrum-analyzer  dynamic  range. 

Section  V  constitutes  the  principal  experimental  section  of  the  report, 
in  which  results  consistent  with  scattering  from  both  surface  roughness  and  sur¬ 
face  compounds  are  obtained.  The  section  outlines  the  experimental  methods  that 
were  used  to  obtain  in-plane  scattered-energy  distributions,  and  the  analytical 
methods  that  were  used  to  interpret  them.  Included  is  a  discussion  of  scatter¬ 
ing  associated  with  prism  input  coupling  and  of  experimental  attempts  to  polish 
two  opposite  waveguide  edges  for  end-fire  input  and  output  coupling.  Then  we 
present  experimental  results  for  seven  different  waveguide  samples  used  in  the 
program,  fabricated  by  various  diffusion  treatments.  The  results  are  used  to 
predict  very  good  performance  levels  for  waveguides  having  smaller  values  of 
surface  index  change  An  and  larger  values  of  diffusion  depth  1)  than  those  we 
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fabricated.  These  conditions  may  be  achieved  using  longer,  hotter  diffusion 
treatments  (say  8  h  at  1000®C),  and  the  use  of  such  conditions  is  the  principal 
recommendation  of  this  report. 

Section  VI  summarizes  the  results  and  conclusions  of  the  program,  while 
Appendices  A-E  expand  on  various  items  mentioned  in  the  body  of  the  report.  The 
most  important  of  these  is  a  theoretical  analysis  of  scattering  enhancement 
associated  with  the  use  of  prism  input  coupling.  This  enhancement  is  produced 
by  the  interaction  of  the  evanescent  field  of  incoming  light  reflected  at  the 
base  of  the  prism  with  scattering  centers  at  the  surface  of  the  waveguide.  An 
unfortunate  finding  is  that  the  enhancement  factor  increases  for  waveguides 
having  index  and  depth  parameters  conducive  to  better  scattering  performance. 

The  result  is  that  good  and  bad  waveguides  tend  to  produce  similar  in^plane 
scattered  energy  distributions  when  obtained  using  prism  input  coupling.  This 
feature  obscures  the  dependence  of  waveguide  quality  on  the  diffusion  treatment 
and  caused  us  to  report  erroneously  in  Ref.  1  that  this  dependence  was  minimal. 

We  now  believe  that  the  best  Ti-diffused  LiNbO«  waveguides  will  be 

o 

formed  from  thinner  Ti  films  (say  150  A)  using  longer  diffusion  treatments  (say 
8  h)  at  higher  temperatures  (say  i000°C) .  Some  post  diffusion  polishing  may 
be  helpful,  especially  if  inhomogeneous  mixtures  of  lithium-titanate  compounds 
remain  at  the  waveguide  surface  following  the  heat  treatment.  This  is  more 
likely  for  thicker  Ti  films,  and  may  also  be  influenced  by  the  diffusion  atmos¬ 
phere,  oxygen  or  argon,  though  this  has  not  been  studied. 

Post-diffusion  polishing  may  also  be  helpful  in  reducing  surface  rough¬ 
ness  caused  by  the  Ti  diffusion,  though  generally  surface  roughness  levels  are 
too  low  to  produce  significant  scattering  except  in  waveguides  having  relatively 
large  values  of  the  ratio  An/D  (say  ^  0.01  pm  ^) ,  If  such  waveguides  are  re¬ 
quired  because  of  system  considerations,  we  point  out  that  good  quality  wave¬ 
guides  can  result  from  heavily  diffused  samples  by  substantial  post-dif f usion 

1-  K- 
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II.  SOURCES  OF  SCATTERING  IN  Ti  DIFFUSED  LiNbO^  WAVEGUIDES 

In  this  section  we  discuss  potential  sources  of  scattering  in  Ti:LlNb03 
waveguides  that  we  have  been  able  to  Identify.  Not  all  of  them  have  been  studi¬ 
ed  in  detail,  and  not  all  of  them  are  viewed  to  be  significant.  But  it  is  ap¬ 
propriate  to  generate  a  list  to  quantify  the  scope  of  the  scattering  problem  and 
to  use  as  a  reference  when  discussing  the  work  of  the  program. 

The  sources  of  scattering  are  divided  into  two  groups:  fabrication-re¬ 
lated  sources  and  diffusion-related  sources.  The  former  are  associated  with  the 
human  role  in  preparing  waveguides;  the  latter  are  related  to  the  physics  of  the 
diffusion  process,  and  would  be  expected  to  be  the  more  difficult  sources  to 
eliminate. 


FABRICATION-RELATED  SOURCES  OF  SCATTERING 


Surface  Roughness 

Surface  roughness  has  been  the  potential  source  of  waveguide  scatter¬ 
ing  most  often  treated  in  the  literature.^  In  the  context  of  this  report, 

it  refers  to  the  deviation  of  the  waveguide  surface  from  an  ideal  plane.  The 
simplest  model  of  a  rough  surface  is  indicated  in  Fig.  1.  The  average  deviation 
of  the  surface  height  from  the  ideal  plane  is  xero,  while  the  rms  average  devia¬ 
tion  is  termed  a.  The  parameter  U  is  the  lateral  distance  over  which  the  sur¬ 
face  height  function  is  correlated.  It  is,  in  simpler  terms,  the  average  wave¬ 
length  of  the  surface-roughness  undulations.  One  may  think  of  the  waveguide 
thickness  as  being  modulated  by  the  roughness  fluctuations.  This  causes  propor¬ 
tional  fluctuations  in  the  waveguide  mode  index.  These  fluctuations  distort  the 
wavefronts  of  guided  light  and  cause  scattering. 

A  perhaps  more  realistic  view  of  a  polished  surface  is  shown  schemati¬ 
cally  in  Fig.  2.,  where  the  effects  of  subsurface  polishing  damage  are  indicated. 
Additional  features  not  shown  in  Fig.  1  are  voids  and  cracks  caused  by  imperfect 
polishing,  in  which  successively  smaller  abrasives  are  not  used  long  enough  to 
remove  the  damage  caused  by  the  previous  abrasives.  Subsurface  damage  phenomena 
are  known  to  exist  for  polished  glass  and  probably  exist  for  LiNbO-^,  though  the 
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phenomena  may  be  very  different  because  of  the  crystalline  nature  and  different 
mechanical  properties  of  LiNbO^  relative  to  those  of  glass.  We  emphasize  that 
Fig.  2  is  schematic  in  nature.  We  know  of  no  work  that  has  been  performed  to 
study  the  nature  of  polishing  damage  in  LiNbO^. 


Crystalline  Imperfections 

Nassau,  Levinstein,  and  Loiacono  have  stated  that  solid -phase  in¬ 
clusions,  gaseous  inclusions,  prominent  low-angle  grain  boundaries,  twin  planes, 
and  dislocations  have  been  revealed  in  pulled  LiNbO^.  They  add  that  all  except 
dislocations  can  be  eliminated  by  careful  crystal  growth.  They  present  photo¬ 
graphs  showing  10-|jm  sized  triangular  and  hexagonal  inclusions  which  they 
describe  as  metal  particles  from  the  crucible. 

We  have  not  done  a  careful  study  of  these  phenomena,  but  our  observa¬ 
tions  tend  to  be  supportive.  Samples  used  in  the  program  have  been  selected- 
acoustic-grade  material  from  Crystal  Technology,  optically  polished  to  ’^/lO  on 
both  sides  to  facilitate  inspection.  The  most  often  observed  imperfections  are 
macroscopic  refractive-index  striations  that  run  for  centimeter  lengths  parallel 
to  the  optic  axis.  The  manufacturer  has  been  cooperative  in  providing  us  with 
crystals  that  do  not  exhibit  this  phenomenon,  which  may  reflect  imperfect  poling 
or  fluctuations  in  crystal  stochiometry . 


Titanium  Nonuniformity 

Fluctuations  in  the  titanium  concentration  in  diffused  waveguides  will 
Lo  associated  with  refractive-index  fluctuations  that  can  scatter  light.  Some 
potential  sources  of  Ti  nonuniformity  are  diffusion  related  and  are  discussed 
below.  One  source  that  is  fabrication  related  has  to  do  with  the  deposition  of 
a  Ti  film  for  diffusion  on  an  imperfectly  cleaned  LiNbO^  substrate.  At  the  high 
temperature  required  for  dif fusion,  dust  and  surface  residues  will  oxidize  and, 
possibly,  mar  the  uniformity  and  diffusion  characteristics  of  the  overlaying  Ti 
film.  These  variations  in  the  film  later  produce  variations  in  the  diffused 
layer  which  scatter  light. 
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Figure  3  shows  a  Nomarskl  Micrograph  of  the  edge  of  a  Ti-diffused 
LiNbO^  waveguide-  The  center  of  the  waveguide  is  defect  free.  The  artifacts 
presumably  result  from  poor  cleanliness  at  the  sample  edge.  The  micrograph  re¬ 
veals  circular  structures  of  5-10  ym  diameter  that  could  be  areas  where  the  Ti- 
film  exploded  away  at  the  high  diffusion  temperature.  Also  noted  are  smaller 
reflective  structures  where,  for  unknown  reasons,  the  Ti  was  not  diffused  into 
the  surface. 
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Fig.  3.  Nomarski  mic rograph  showing  artifacts 
near  the  edge  of  a  Ti-diffused  LiNb03 
waveguide.  (500X) 
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DIFFUSION  RELATED  SOURCES  OF  SCATTERING 


Li-Ti-0  Compound  Formation 

Burns,  Klein,  West  and  Plew  have  published  evidence  that  the  diffused 
Ti  concentration  in  Ti:LiNbO^  waveguides  shows  a  peak  extending  to  about  0.3  pm 

(7)  ^ 

below  the  surface.  This  peak  is  imposed  upon  the  anticipated  Gaussian 

diffusion  profile  extending  2  pm  below  the  surface.  The  authors  postulate  that 
the  excess  Ti  concentration  is  bound  in  the  form  of  Li-Ti*0  compounds  such  as 
Li2TiO^  or  Li2Ti^0^.  These  compounds  form  a  dilute  mixture  with  the  LiNbO^ 
host.  Unless  this  mixture  is  homogeneously  distributed  in  the  plane  of  the 
waveguide,  it  will  constitute  a  source  of  waveguide  scattering. 

Microdomains 


Microdomains  in  gro\gn  crystals  are  generally  prevented  by  the  use  of 

/O') 

proper  poling  techniques.^  However,  Ohnishi  has  reported  that  domain  re¬ 

versal  at  the  positive  dipole  (+c)  face  of  LiNbO^  can  occur  at  the  elevated 

(9) 

temperatures  used  to  make  waveguides  by  diffusion.  Also,  Venables  reports 

that  microdomains  can  be  generated  in  LiTaO^  crystals  by  the  process  of  grind¬ 
ing  and  polishing  thin  wafers ^^^^oth  of  these  mechanisms  could  conceivably  in¬ 
troduce  microdomain  scattering  centers  into  LlNbO^  waveguides. 

Growth  Strains  and  Misfit  Dislocations 


Ramaswamy  and  Standley  have  reported  the  observation  of  dislocations 

in  Nb-dif fused  LiTaO^  caused  by  elastic  strain  due  to  the  lattice  mismatch  be- 

^  (11) 

tween  LiNbO^  and  LiTaO^.  Boyd,  Schmidt,  and  Storz  report  that  Ni  diffusion 

^  ^  (12) 
into  a  LiNbO^  rod  produced  strains  sufficient  to  fracture  the  crystal. 

Sugii,  Fukama,  and  Iwasaki  found  a  0.1%  lattice  contractions  caused  by  Ti 

( I  3) 

diffusion  into  LlNbO^.  This  resulted  in  misfit  dislocations  and  cracks  in  the 
diffused  layer  that  could  serve  as  a  source  of  waveguide  scattered  light. 
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LiNb^Og  Formation 


Svaasand»  Eriksrud,  Nakken,  and  Grande  have  observed  the  gradual  con¬ 
version  of  LiNbO^  to  LiNb-Oo  upon  heating  the  crystal  to  the  vicinity  of 

/  I JO 

810°C.  The  latter  phase  occurs  in  the  form  of  small  crystallites  that  were 

demonstrably  effective  scatterers  of  light.  Some  generation  of  LiNb^Og  crystal¬ 
lites  could  occur  during  waveguide  fabrication  as  the  substrate  is  heated  to  or 
cooled  from  the  diffusion  temperature,  giving  rise  to  optical  waveguide  scat¬ 
tering.  In  this  connection,  note  that  the  outdif fusion  process  is  described  by 
the  chemical  reaction  3LiNbO« Li«0  +  LiNb^Oc.  Outdif fusion  during  waveguide 
formation  may  thus  be  conducive  to  scattering  by  separated  phase  material,  pro¬ 
vided  this  material  is  not  homogeneously  distributed  in  the  plane  of  the  wave¬ 
guide.  Figure  4  shows  the  surface  degradation  that  can  result  from  a  heat  treat¬ 
ment  specifically  intended  to  produce  the  separated  phase  LiNb^O-.. 

3  o 

Globular  Metal  Films 

It  is  known  that  vacuum  evaporation  of  thin  layers  of  Sn  onto  heated 
substrates  can  produce  globular  films  containing  particles  of  the  order  of  the 
wavelength  of  visible  light  in  size.^^^^  If  this  phenomena  were  to  occur  for 
the  case  of  Ti  films  evaporated  on  LiNbO^  substrates,  the  resultant  diffused 
layer  could  contain  index  inhomogeneities  capable  of  scattering  light. 


Statistical  Fluctuations 


Even  if  perfection  is  achieved  in  the  fabrication  of  a  uniform  Ti  film 
and  in  the  subsequent  diffusion  process,  a  finite  level  of  scattering  will  still 
result  from  statistical  fluctuations  in  the  concentration  of  Ti  throughout  the 
waveguide.  This  is  the  lowest  level  of  scattering  which  is  theoretically 
possible. 
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III.  METHODS  FOR  STUDY  OF  SCATTERING  SOURCES 


In  this  section  we  will  present  a  compendium  of  methods  that  may  be 
used  to  study  sources  of  waveguide  scattering.  We  have  employed  many  of  these 
methods  to  varying  degrees,  and  we  will  indicate  when  appropriate  our  experi¬ 
ences  with  them.  Our  results  with  the  methods  that  we  employed  most  often  will 
however,  be  covered  in  later  sections  in  more  detail. 

We  classify  methods  for  studying  scattering  sources  as  being  either 
direct  or  indirect,  depending  on  whether  the  method  provides  actual  information 
about  the  source  or  information  that  must  be  inferred  from  a  related  observa¬ 
tion.  Microscopy  in  its  varied  forms  would  be  an  example  of  a  direct  method, 
while  etching  studies  would  be  an  example  of  an  indirect  method.  Direct  method 
can  be  further  divided  into  topographical  approaches,  which  look  for  surface 
structures,  and  compositional  approaches,  that  look  for  spatial  variation  in 
chemical  make  up. 


DIRECT  METHODS 


Topographical  Approaches 


Noroarski  Microscopy 

The  microscope  is  the  most  obvious  tool  for  examining  surface  topo¬ 
graphy,  but  in  the  case  of  highly  polished  surfaces  showing  little  structure  it 
leaves  much  to  be  desired  as  a  means  for  resolving  small  surface  height  differ¬ 
ential.  An  adaptation  known  as  dif ferential-interfercnce-contrast  microscopy 
(rr  Nomarski  microscopy)  has  the  capability  for  converting  phase  information 
imposed  on  an  optical  beam  by  reflection  from  a  nearly  perfect  surface  to 
amplitude  information,^^  ^  This  makes  it  possible  to  view  topographical  struc¬ 
tures  as  small  . several  Angstroms  provided  the  slopes  are  sufficiently  steep. 

Our  main  use  of  the  Nomarski  microscope  throughout  this  program  has 
been  to  evaluate  the  topography  of  diffused  waveguide  surfaces  directly  after 
diffusion.  We  usually  have  found  a  granular  texture  such  as  that  shown  in 
Fig.  5  which  we  spt^culated  was  one  result  of  the  Li-Ti-0  compound  formation 
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Fig,  5,  Nomarski  micrograph  showing  the  granular 
surface  texture  associated  with  Ti  dif¬ 
fusion  (1000  A  Ti  film,  26  h  diffusion, 

500X  magnification) ,  Vertical  line  shows 
the  demarkation  between  the  diffused  region 
(right)  and  the  undiffused  region  (left). 


described  by  Bums  and  coworkers.  Other  have  noted  a  change  in  this  granu¬ 

larity  with  increasing  diffusion  time,  with  the  appearance  of  the  surface  eenex> 

(17) 

ally  becoming  smoother.  We  have  tentative  evidence  that  this  is  the  case 

but  have  not  performed  a  careful  study.  Our  approach  has  been  to  achieve  wave¬ 
guide  quality  enhancement  by  polishing  away  the  granular  layer.  Our  experi¬ 
ments  are  summarized  in  Ref.  1. 


TIR  Microscopy 

A  different  optical  microscopic  technique  has  been  employed  by  H.  E. 

/  1  Q\ 

Bennett  and  coworkers  at  the  Naval  Weapons  Center.  As  shown  in  Fig.  6,  a 

laser  beam  is  reflected  from  a  polished  surface  at  an  angle  greater  than  the 
critical  angle.  The  evanescent  field  is  scattered  by  surface  structures  and 
detected  with  a  conventional  optical  microscope.  This  is  similar  to  dark  field 
microscopy.  We  tested  the  method  and  found  that  it  was  very  sensitive  to  sur¬ 
face  dust,  but  we  were  not  sufficiently  encouraged  by  our  observations  to  use 
the  technique  in  place  of  Nomarski  microscopy. 

Multiple -Beam- Interference  Microscopy 

We  also  tested,  briefly,  the  use  of  multiple-beam-interference  methods 
for  studying  surface  topography.  These  methods  generally  employ  a  conventional 
microscope  to  image  the  wavefront  formed  by  multiple  reflections  between  a  test 
surface  and  a  reference  surface.  The  two  surfaces  are  coated  with  a  highly 
reflective  layer  to  increase  the  number  of  reflections  of  light  waves  trapped 
between  them.  This  amplifies  the  phase  distortion  associated  with  surface 

o 

topography  to  the  extent  that  5  A  roughness  can  be  observed  with  2  pm  horizontal 
(19) 

resolution.  However,  the  time  available  during  this  program  to  establish 

this  measuring  capability  was  not  sufficient  to  achieve  the  reported  sensitivity. 
Using  a  commercial  device  for  measuring  the  thickness  of  thin  films,  called  an 

o 

Angstrometer ,  we  achieved  a  100  A  vertical  resolution  during  preliminary  investi¬ 
gations  (see  Fig,  7).  This  was  not  adequate  to  resolve  the  surface  roughness  of 
our  LiNbO^  samples,  and  we  did  not  pursue  the  method  further. 
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EYEPIECE 


Fig.  6.  Total-internal -ref lection  microscopy. 
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The  final  microscopic  technique  that  we  tested  was  that  of  the  scan¬ 
ning  electron  microscope  (SEM) ,  Actually  we  employed  this  equipment  in  several 
modes  of  operation  during  the  program,  but  the  one  we  describe  here  is  a  shad¬ 
owing  method  for  looking  at  the  topography  of  relatively  smooth  surfaces.  In 
this  technique  a  sample  of  LiNbO^  is  coated  with  gold  at  glancing  incidence 

1°)  to  allow  shallow  topographical  structures  to  shadow  the  surface.  Fig.  8 
shows  the  resultant  SEM  micrograph  of  a  dusty  surface,  in  which  the  shadowing 
effect  is  clearly  evident.  However  the  technique  becomes  less  useful  as  the 
slope  of  the  surface  irr<=‘[^,ularities  becomes  small  in  comparison  to  the  angle 
of  incidence  of  the  evaporated  beam,  and  this  was  apparently  the  case  for  clean 
portions  of  the  sample  examined.  There  was  no  variation  in  surface  contrast 
that  could  obviously  be  attributed  to  the  topography  of  the  polished  surface. 


Talystep  Prof ilometrv 


The  last  topographical  approach  that  was  employed  to  study  surface 

roughness  was  the  Talystep  instrument  operated  by  Drs.  Ted  Vorburger  and  Clayton 

Teague  of  the  National  Bureau  of  Standards.  This  instrument  has  a  maximum  mag- 
6 

nification  of  10  ,  making  possible  the  detection  of  rms  surface  roughness  as 

o 

small  as  1  A.  At  this  level  of  magnification  the  instrument  was  used  with  a 
chisel- shaped  stylus  with  approximate  dimensions  0.1  pm  x  1  pm. 

The  instrument  was  used  to  examine  the  topography  of  a  polished  LiNbO^ 
substrate  purchased  from  Crystal  Technology.  No  processing  of  the  substrate  was 
performed  prior  to  the  examination.  Ted  Vorburger  of  NBS  made  three  1.5-nira  long 
scans  of  the  surface  both  parallel  and  perpendicular  to  the  optical  axis.  He 
reported  that  the  surface  was  among  the  smoothest  ever  examined  at  NBS,  with 

o 

less  than  3  A  rms  roughness  and  an  autocorrelation  length  of  about  36  pm.  The 
autocorrelation  length  may  have  been  slightly  longer  when  measured  parallel  to 
the  optical  axis.  There  was  possibly  an  additional  short  wave  correlation 
length  of  about  1  pm,  but  this  was  too  close  to  the  lateral  resolution  of  the 

o 

system  to  be  sure.  For  that  matter,  a  significant  fraction  of  the  3  A  rough¬ 
ness  could  have  its  origin  in  system  noise.  The  complete  NBS  report  is  Appen¬ 
dix  A  of  this  report.  Surface  evaluation  was  repeated  after  an  interval  of 
73  days  with  similar  results. 
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Fig.  8 


SEM  shadowgraph  of  a  dusty  LlNbO« 
surface.  (2000X) 
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The  sensitivity  of  the  Talystep  instrument  makes  it,  in  our  opinion, 
the  best  tool  for  surface-roughness  studies  in  this  type  of  work.  Unfortu- 
nately,  awareness  of  this  instrument  and  its  capabilities  came  too  late  to  make 
further  use  of  it  during  the  program. 


Other  methods  for  direct  observation  of  waveguide  scattering  centers 

look  for  variations  in  the  composition  of  the  waveguide  materials  rather  than 

for  variations  in  waveguide  topography.  One  such  method  is  e-beam  x-ray  analy- 

in  wliicli  a  micron-sized  electron  beam  is  used  to  bombard  a  sample  with  the 

objective  of  creating  x-rays  characteristic  of  the  materials  encountered . ^ 

This  tecl^nique  has  been  used  to  measure  Ti-concentration  profiles  in  planar  and 
(  20  21) 

channel  waveguides.  ’  However,  in  a  scattering  program  the  objective  is 
to  measure  nonuniformities  in  the  Ti  concentration,  rather  than  the  average  con¬ 
centration,  and  a  much  greater  sensitivity  is  required  of  the  technique. 

We  employed  e-beam  x-ray  analysis  to  look  for  Ti  nonuniformities  in  a 
diffused  waveguide,  in  an  undiffused,  unoxidized  Ti  film,  and  in  an  undiffused, 
oxidized  Ti  film.  Both  films  were  on  LiNbO^  substrates.  The  only  situation  for 
which  we  observed  a  Ti  nonuniformity  was  in  the  case  of  the  undiffused  film  that 
had  been  oxidized  at  600°C  for  1  hour.  The  initial  film  was  deposited  on  an  im¬ 
perfectly  clean  surface  and  served  to  decorate  the  surface  residues.  After  ther¬ 
mal  oxidation,  further  decoration  occurred,  and  new  imperfections,  such  as  micron 
sized  bubbles,  were  observed.  Tt  might  be  anticipated  that  the  oxidation  process 
would  cause  vaporization  of  organic  residues  which,  in  turn  would  explosively  re¬ 
move  some  of  the  overlaying  Ti  film.  We  did  not  acquire  any  evidence  for  this, 
even  when  we  scanned  the  electron  beam  through  highly  decorated  regions  of  the 
oxidized  film.  However,  we  did  obtain  one  scan  of  the  720  pm  length  which  shows  a 
S-pm-wide  region  of  excess  Ti.  This  scan  is  shown  in  Fig.  9.  It  is  not  known 
if  the  excess  Ti  is  associated  with  the  oxidation  process  or  with  the  initial 
Ti -depos i t ion  process.  A  scan  of  an  unoxidized  Ti  film  revealed  no  similar 
feature,  although  In  measurements  of  this  type  it  must  be  remembered  that  a 
millimeter-sized  scan  of  a  micron-sized  beam  does  not  cover  much  surface  area. 
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e-beam.  X-ray  analysis  data  showing  the  presence  of  a 
region  of  excess  Ti  on  an  oxidized  but  undiffused  Ti 
f ilm/LiNbO^  substrate . 


The  chances  of  seeing  a  5  um  particle  using  a  1  mm  scan  are  good  only  if  the 

2  -1  “2 

particle  density  is  large  compared  to  (0.005  x  1  mm  )  =  200  mm  .  We  did  not 

make  enough  scans  of  the  surface  to  fully  test  the  potential  of  this  method  of 
measurement.  In  retrospect  this  would  have  been  worth  doing. 

Secondary  Ion  Mass  Spectroscopy 

Another  technique  for  looking  at  compositional  variations  in  waveguide 
materials  is  Secondary  Ion  Mass  Spectroscopy  (SIMS).^^^  In  this  technique  the 
compositional  profile  is  determined  in  depth  as  material  is  ejected  from  the 
surface  and  analyzed  using  mass-spectroscopic  techniques.  This  method  does  not 
usually  offer  high  spatial  resolution  in  the  plane  of  the  waveguide,  but  it 
does  so  in  depth.  In  this  sense,  SIMS  is  complimentary  to  e-beam  x-ray  analy¬ 
sis  which  averages  over  about  2-ym  depth  but  affords  good  lateral  resolution. 

The  most  important  use  of  SIMS  in  studying  optical  waveguides  was  re¬ 
ported  by  Burns  and  coworkers.  They  observed  an  anomalously  large  Ti  con¬ 

centration  in  the  top  0,3  pm  of  a  diffused  waveguide.  They  attributed  this  to 
the  formation  of  Ll-Ti-0  compounds  during  diffusion.  We  have  speculated  that 
the  granular  structure  shown  on  the  waveguide  surface  of  Figure  3  is  a  further 
indication  of  this  phenomenon.  The  importance  of  this  Li-Ti-0  surface  contami¬ 
nation  to  waveguide  scattering  is  inferred  from  the  fact  that  waveguide  quality 
degrades  when  thicker  Ti-films  are  diffused,  but  this  should  ultimately  be 
tested  by  a  series  of  measurements  in  which  waveguide  quality  is  compared  to 
Li-Ti-0  contamination  as  measured  using  SIMS. 


INDIRECT  METHODS 


Waveguide  Scattering 

We  now  move  to  a  discussion  of  methods  for  studying  scattering  sources 
in  which  information  must  be  inferred  from  observations  which  are  indirectly 
influenced  by  scatterers.  The  method  that  we  have  used  the  most  is  in-plane 
waveguide  scattering.  Measured  scattering  levels  can  be  related  by  theory  to 
characteristic  scattering  parameters.  The  parameters  can  then  be  compared  to 
those  which  would  be  expected  for  various  sources.  The  information  obtained  by 
this  procedure  is  model-dependent  and  should  therefore  be  used  to  corroborate 
direct  information  whenever  this  is  possible.  We  will  have  much  to  say  about 
waveguide  scattering  measurements  and  their  interpretation  in  subsequent  sec¬ 
tions  of  this  report. 

Electric-Field  Induced  Scattering 

Occasionally*  modifications  of  the  standard  in-plane  waveguide  scatter¬ 
ing  experiment  can  be  performed  to  obtain  more  concrete  information  about  par¬ 
ticular  scattering  sources.  For  example,  we  have  applied  electric  fields  to 
LiNbO^  waveguides  with  the  expectation  of  enhancing  scattering  from  any  micro- 
domains  that  may  be  present.  Our  first  observations  reported  in  Ref.  1  showed 
a  positive  effect.  More  recently  we  have  determined  that  the  electric-field 
enhanced  scattering  that  we  observed  was  due  to  fringing-f ield  effects  at  the 
electrodes  used  to  apply  the  electric  field  to  the  waveguide.  This  conclusion 
is  based  on  the  observation  that  field-enhanced  scattering  decreases  as  the 
electrode  gap  becomes  significantly  wider  than  the  guided  beam. 

Indif fused  and  Outdif fused  Waveguide  Comparison 

Another  use  that  can  be  made  of  the  in-plane  waveguide  scattering  ex¬ 
periment  is  to  compare  the  scattering  from  both  Tl-indlf fused  LiNbO^  waveguides 
and  Ll20-outdif fused  LiNbO^  waveguides.  Any  excess  scattering  in  the  former 
may  be  ascribed  to  a  source  that  is  related  to  the  T1  diffusant.  Similar  scat¬ 
tering  levels  suggests  a  source  of  scattering  common  to  the  host  crystal,  such 
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as  surface  roughness  or  the  occurrence  of  the  separated  phase  LiNb^Og.  The 
situation  is  complicated  by  scattering  differences  associated  with  prism  coupl¬ 
ing  and  by  the  fact  that  Ti-diffusion  can  influence  surface  roughness  and 
LiNbgOg  formation.  We  will  discuss  this  further  in  connection  with  our  experi¬ 
mental  results  in  Sec. V. 

Polishing  Studies 

Waveguide  scattering  can  be  used  in  "before  and  after"  studies  in 
which  the  waveguide  is  repeatedly  modified  and  examined  to  see  if  the  modifica¬ 
tion  produced  a  reduction  in  scattering.  Our  best  example  of  this  is  the 
improvement  in  waveguide  quality  that  results  from  polishing  the  surface. 

o 

This  improvement  is  most  pronounced  when  thick  700  A)  Ti  films  are  diffused 
to  make  the  waveguide.  We  have  attributed  the  improvement  with  polishing  to 
removal  of  Li-Ti-0  compounds  at  the  surface  as  described  in  Ref.  7. 

Reflected  Light  Scattering 

Tn  addition  to  the  use  of  waveguide  scattering  to  obtain  information 
about  sources,  we  have  employed  conventional  beam  scattering  from  waveguide  sur¬ 
faces  to  obtain  information  specifically  about  surface-roughness  parameters.  The 
results  will  be  discussed  in  detail  in  Appendix  E. 

Thermal  and  Chemical  Decoration  Techniques 

There  are,  in  addition  to  these  optical  approaches  based  on  light  scat¬ 
tering,  several  non-optical  approaches  for  obtaining  indirect  information  about 
scattering  sources.  These  are  generally  thermal  or  chemical  techniques  for 
decorating  flaws  that  might  not  otherwise  be  visible.  The  most  well  known  of 
these  is  wet  chemical  etching. 

We  have  mentioned  several  etching  studies  of  LiNbO^  used  as  a  means  for 
observing  microdomains,  Our  etching  experiments  gave  no  clear  indication 

of  the  presence  of  microdomains,  but  we  did  observe  some  interesting  effects. 

The  most  significant  is  that  Ti  diffusion  makes  the  LiNbO^  surface  more  etch 
resistant.  Figure  10  shows  the  results  of  an  experiment  in  which  a  sample  con¬ 
taining  a  'V'  130  pm  wide  Ti-dif  fused  channel  was  etched  in  hot  The 
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material  surrounding  the  channel  has  been  removed  to  a  greater  extent  than  the 
material  in  the  channel. 

Figure  11  shows  the  results  of  an  etch  in  hot  concentrated  H«SO.  of  a 

o  2  4 

sample  that  was  coated  with  150  A  Ti  and  diffused  to  form  a  waveguide.  Al¬ 
though  the  sample  was  coated  uniformly,  the  resulting  waveguide  surface  showed 
regions  of  high  and  low  granularity.  When  etched,  the  regions  of  low  granu¬ 
larity  showed  the  characteristic  pattern  that  we  believe  to  be  associated 
with  the  host  material.  The  regions  of  high  granularity  were  etch  resistant. 

This  result  is  supportive  of  the  idea  that  Li-Ti-0  compounds  are  re¬ 
sponsible  for  etch  resistance  and  that  surface  granularity  is  an  indication  of 
their  presence.  Nonuniform  distribution  of  the  compounds  can  cause  light  scat¬ 
tering.  The  granular  texture  of  the  diffused  surface,  both  before  and  after 
etching,  suggests  that  the  compounds  are  indeed  distributed  nonuniformly  in 
space.  The  characteristic  size  of  the  granular  texture  is  on  the  order  of  1  um 
before  etching.  After  etching,  the  characteristic  size  estimated  from  Fig.  11 
is  more  like  20  pm. 

Figure  12  shows  the  results  of  an  etching  experiment  carried  out  using 
1  part  HF  and  2  parts  HNO^  at  90®C  for  70  min.  The  sample  was  a  scrap  piece  of 
LiNbO^  of  uncertain  history.  It  may  or  may  not  have  been  Tl-dlf fused,  although 
notes  on  the  experiment  do  not  mention  any  particular  etch  resistance  that  may 
have  been  associated  with  the  presence  of  diffused  Ti . 

The  occurrence  of  etch  pits  or  hillocks  similar  to  the  large  struc- 

(9) 

tures  in  Fig.  12  have  been  noted  by  others,  usually  in  z-cut  crystals.  Most 
often  they  are  isolated  and  occur  at  the  terminus  of  dislocations.  In  Fig.  12 
the  etch-induced  structures  are  contiguous.  Moreover,  the  crystal  is  y-cut 
rather  than  z-cut.  This  may  further  effect  their  interpretation. 

One  tantalizing  way  to  interpret  Fig,  12  is  to  associate  the  size  and 
distribution  of  etch  structures  with  the  size  and  distribution  of  scattering 
centers.  This  is  reasonable  if  the  scattering  centers  give  rise  to  the  etch 
structures.  The  size  of  the  larger  etch  structures  are  15  pm.  Interestingly, 
this  is  about  the  size  of  the  scattering  centers  as  we  measure  them  using  in¬ 
plane  waveguide  scattering. 
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11.  Results  of  etching  a  uniformly  Ti-diffused 
LiNbO^  sample  in  hot  concentrated  H2S0^. 
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Fig.  12.  Results  of  etching  I.iNbO_  in  hot 
HF/HNO^. 


IV.  EFFECTS  OF  SCATTERING  SOURCES  ON  WAVEGUIDE 
AND  SPECTRUM  ANALYZER  PERFORMANCE 


In  this  section  we  use  simple  models  for  sources  of  scattering  to  pre¬ 
dict  their  effect  on  waveguide  quality  and,  in  particular,  on  the  dynamic  range 

(2) 

of  a  spectrum  analyzer.  These  calculations  are  helpful  in  identifying  t3rpes  of 
scatterers  that  are  likely  to  be  a  particular  problem  in  limiting  spectrum  ana¬ 
lyzer  pertormance  and  they  provide  us  with  guidelines  to  use  in  interpreting 
our  experimental  results  of  the  next  sections. 

CALCULATION  OF  SPECTRUM-ANALYZER  DYNAMIC  RANGE 


Figure  13  shows  the  front  and  back  focal  region  of  a  spectrum  ana¬ 
lyzer  design  that  we  will  use  in  calculating  dynamic  range.  The  object  of  the 
calculation  is  to  predict  the  scattered  signal  at  any  detector  relative  to  the 
unscattered  signal  at  the  central  detector.  We  consider  scattering  events  that 
occur  in  both  the  diffracted  and  undiffracted  beams  both  before  and  after  the 
lens.  But  we  do  not  consider  scattering  events  that  occur  before  the  acoustic 
wave  in  the  front  plane  of  the  lens.  To  further  simplify  the  calculation,  we 
take  the  lens  to  be  infinitely  thin  and  positioned  at  the  center  of  the  real, 
distributed  lens. 

The  specifications  for  the  spectrum  analyzer  shown  in  Fig,  13  were 
provided  by  the  Air  Force;  however,  we  can  base  our  derivations  on  a  general 
spectrum  analyzer  for  which  F  is  the  lens  focal  length,  W  is  the  unfocused  beam 
width,  6y  is  the  detector  width,  and  Ay  is  the  separation  of  a  particular  detec¬ 
tor  from  the  central  detector.  The  z-axis  is  the  direction  of  the  diffracted 
beam  and  z  is  the  distance  from  the  detector  array. 

To  keep  the  scattering  mechanism  as  general  as  possible  for  the  pre¬ 
sent,  we  describe  it  by  the  differential  scattering  cross  section  a((()),  defined 
as  follows:  Tf  I  is  the  power  per  unit  width  W  of  the  waveguided  beam,  then 

AP  =  Ia((t))A^  (1) 
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Fig,  13.  Geometry  used  In  calculating  spectrum  analyzer  performance. 
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is  the  power  scattered  in  plane  at  an  angle  (p  into  an  angular  sector  A4>.  Note 
that  as  a  result  of  the  waveguide  geometry,  a(4))  has  dimensions  of  length.  This 
is  in  contrast  to  the  case  for  scattering  in  a  3-dimensional  geometry,  where  a 
has  dimensions  of  area  and  I  is  power  per  unit  area  of  the  incident  beam.  If 
there  are  Z  scattering  centers  per  unit  area,  then  the  total  power  scattered 
into  the  sector  A(j)  is 


=  ZLW  lo(i>)A(t)  ,  (2) 

where  the  product  LW  is  the  scattering  area. 


Contribution  to  Scattered  Power  from 
the  Focused  Beams 


To  apply  these  concepts  to  a  calculation  of  spectrum-analyzer  dynamic 
range ,  consider  the  contribution  to  scattered  signal  from  scattering  events  that 
occur  between  the  center  of  the  waveguide  lens,  assumed  thin,  and  the  detector 
array.  Figure  14  defines  the  scattering  parameters  for  a  strip  located  be¬ 
tween  z  and  z  +  dz,  where  z  is  the  distance  from  the  array.  For  rays  scattered 
at  height  y  and  incident  on  the  array  at  height  y’,  the  scattering  angle  is 

(p  =  tan  ^[(y’-y)/z]  +  tan  ^[y/z]  .  (3) 

The  total  scattered  power  from  the  strip  into  a  detector  extending  from 

y*  =  Ay  -  6y/2  to  y’  =  Ay  +  6y/2  is  obtained  by  integrating  over  all  scattering 

angles  subtended  by  the  detector,  all  initial  ray  heights,  and  all  thin  strips. 

For  convenience,  we  convert  the  <()  integration  to  a  y’  integration  using 

d(|)  =  (9(})/3y  ’  )dy  ’  .  We  have,  then 

F  ztanO  Ay+6y /2 

AP  =  /  dz  /  dy  EI(y,z)  /  dy'  a[(()(y,y '  ,z)  ]  (3iJ>/3y')  .  (4) 

o  -ztan0  Ay-6y/2 

In  this  expression,  0  is  the  half-angle  defined  by  the  converging  beam.  For 
most  cases  of  interest,  the  differential  scattering  cross  section  will  not  vary 
much  with  y*  across  the  span  of  the  detector  element.  Consequently,  it  may  be 
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replaced  by  its  value  at  y*  *  Ay  and  taken  outside  the  y*  Integration.  The  y' 
integral  is  then  just  the  difference  in  the  scattering  angle  (j)  when  evaluated 
at  the  edges  of  the  detector.  We  are  left  with 

F  ztan0 

AP  =  /  dz  /  dy  II(y,z)a[4i(y,Ay,z)]  (5) 

o  -ztan0 

X  {tan  ^[(Ay  +  6y/2  -  y)/z] 

-tan  ^  f (Ay  -  6y/2  -  y)/z]}  , 

It  is  desirable  to  obtain  an  analytic  approximation  to  this  integral,  one  which 
at  least  contains  the  correct  dependences  on  the  parameters,  though  it  may  be 
off  by  numerical  factors  on  the  order  of  two.  Examination  of  Fig.  14  sug¬ 
gests  that  the  angle  subtended  by  the  detector,  given  by  the  quantity  in  curly 
brackets  in  Eq .  (5>,  is  approximately  described  by  the  value  obtained  when 
Ay  and  y  are  zero.  This  approximation  is  best  for  detectors  near  the  optical 
axis  which  are  of  particular  interest  since  they  determine  the  spectrum  analyzer 
dynamic  range. 

We  also  simplify  Eq .  (5)  by  using  the  geometrical  optics  form  for 
the  beam  intensity,  given  by 

I(y»z)  =  I^W/2ztan0  in  -ztanO<y<z tan0 ,  (6) 

=  0  elsewhere. 

In  this  expression,  I^W  is  the  total  power  in  the  beam,  assumed  to 
be  approximately  uniformly  distributed  in  y  at  each  value  of  z.  Finally,  the 
scattering  angle  ((),  given  by  Eq ,  (  0  with  y*  =  Ay,  varies  only  slightly  with 

y  for  a  beam  with  a  small  value  of  0  such  as  we  are  considering.  A  good  ap¬ 
proximate  value  for  ((>  is  therefore  the  value  on  axis, 

(ti(y,Ay,z)  =  tan  ^  (Ay/z)  .  (7) 

In  the  limit  of  the  validity  of  these  approximations  the  y-integration 
of  Eq.(S)is  easily  carried  out  with  the  result 
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(8) 


•  1  *  1 

AP  =  2  dz  altan  (Ay/z)]  tan  6y/2z  , 

We  now  identify  two  cases  of  interest;  narrow-angle  scattering,  associated  with 
scattering  centers  of  characteristic  size  a  large  compared  with  the  wavelength 
(Rayleigh-Gans-Debye  (RGD)  limit)  and  wide  angle  or  dipole  scattering,  associ¬ 
ated  with  scattering  centers  smaller  than  the  wavelength  (Rayleigh  limit).  In 
the  latter  case  we  may  take 


o(())) 

-  a(0)  for  TM  modes  , 

(9) 

a((J)) 

2 

-  a(0)cos  for  TE  modes 

In  the  former  case,  we  have  shown  in  Reference  1  that 

a((J>)  =  o(0)  [1  +  ,  (10) 

where  k  =  2it/A,  A  being  the  wavelength  in  the  material. 

We  can  now  evaluate  the  scattered  power  by  inserting  these  expressions 
into  Eq.  (^)»using  Eq.(7)  for  ({).  The  simplest  case  is  that  of  Rayleigh  scattering 
of  TM  modes.  The  differential  scattering  cross  section  is  constant  because  the 
in-plane  scattering  direction  is  always  normal  to  the  mode  polarization.  The 
z-dependence  of  the  Integrand  of  Eq.(8)is  that  of  the  term  tan  ^  6y/2z  alone. 

The  value  of  the  integral,  taken  between  the  center  of  the  lens  and  the  diode 
array,  is 


AP  =  i:i^W5ya(0)£n(2F/6y)  .  (11) 

This  is  an  approximation  to  the  actual  integral.  It  is  valid  in  the  limit 

2F/6y>>l,  which  is  well  satisfied  in  all  examples  of  practical  interest. 

Equation  11  does  not  show  a  dependence  on  Ay  owing  to  our  use  of 

2tan  ^6y/2z  to  describe  the  angle  subtended  by  the  detector.  For  off-axis  detec- 

2  2 

torSf  a  good  approximation  to  the  subtended  angle  is  5yz/(z  +  Ay  ) .  When  this 
is  used  in  Eq. (8)ln  place  of  2tan~^5y/22,  the  resulting  formula  for  scattered 
power  is  found  to  be 
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(12) 


J  2  ?  1/2 

AP  =  Zl^V6ya(0)inl(F  +  Ay  )/Ay  ]  ' 

One  can  phenomenologically  combine  Eqs, (11)  and  (12)  to  obtain  an  expression  for 
scattered  power  that  reduces  to  Eq. (11)  when  Ay  is  zero,  and  yet  does  not  differ 
significantly  from  Eq.(l2)when  Ay  is  not  zero.  The  combined  expression  is 

2  2  2  2  1/2 

AP  =  EI^W6ya(0)  to[(F'^  +  Ay  )/(6y  /4  +  Ay  )]^^  .  (13) 

This  result  is  valid  in  the  limit  of  Rayleigh  scattering  of  TM  modes. 

In  considering  the  Rayleigh  limit  for  TE  polarized  light,  the  inte¬ 
grand  of  Eq,(8)is  reduced  from  the  TM  case  by  the  factor 

cos^<j)  =  z^/(z^  +  Ay^)  ,  (14) 

This  factor  results  from  the  obliquity  of  the  polarization  of  the  scattered 
light  relative  to  the  polarization  of  the  incident  light.  For  on-axls  scatter¬ 
ing,  the  obliquity  factor  is  unity  and  the  scattered  power  is  the  same  as  that 
given  by  Eq.(ll)  for  TM  polarization. 

For  off-axis  scattering,  the  obliquity  factor  forces  the  integrand  of 
Eq. (8)  to  vanish  near  z  =  0.  For  this  reason  we  can  take  tan  ^6y/22  =  6y/2z  in 
Eq.  (8)  without  introducing  much  error  to  the  integral,  even  near  2  =  0.  The  re¬ 
sult  of  integrating  over  z  is  exactly  the  result  of  Eq.  (12)  obtained  for  TM 
polarization.  We  conclude  that  the  modified  expression  of  Eq.  (13) is  valid  in 
the  Rayleigh  limit  for  both  TE  and  TM  polarization.  Differences  associated  with 

obliquity  can  be  determined  by  using  the  more  precise  substitution 
-1  2  2 

tan  6y/2z  (z(Sy/2)/(z  +  Ay  )  in  Eq.  (8), but  these  difference  are  found  to 

modify  Eq.  13  by  a  small  factor  of  about  {2  ^ which  we  ignore. 

We  now  consider  the  case  of  weak  scattering  centers  large  compared 
with  a  wavelength  (RGD  scattering — later  we  will  consider  aspects  of  the 
more  general  case  of  Mie  scattering).  In  the  RGD  limit,  the  z-dependent  part 
of  the  integrand  of  Eq.  (S)is  [1  +  k^a^(tan  ^Ay/z)^]  ^  tan  ^6y/22.  Since  kil»l, 
the  Integrand  is  largest  for  large  z  except  for  the  on-axis  case  Ay  =  0.  This 
suggests  that  we  may  replace  the  inverse-tangent  functions  by  their  arguments 
except  when  Ay  =  0.  In  the  Ay  =  0  case,  small  values  of  z  make  Important 
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contributions  to  the  integral  which  are  not  correctly  accounted  for  by  the 

small-angle  approximation.  If  we  ignore  this  problem  temporarily,  the  approxi- 

2  2  2  2 

mate  integrand  has  the  z  dependence  z/(z  4*  k  a  Ay  ) .  The  resulting  expression 
for  scattered  power  is 

2  22  21/2 

AP  2:i^K6ycrCO)£n[(F'^  +  k  a^Ay  )  '  /kaAy]  ,  (15) 


This  expression  is  infinite  at  Ay  =  0.  However,  we  can  replace  the  logarithmic 

2  2  2  2  1/2 

denominator  kaAy  by  the  term  (k  a  Ay  +  6y  /4)  .  This  is  similar  to  the  phe¬ 

nomenological  approach  used  in  proceeding  from  Eq.  (12)  to  Eq.  (13).  The  added 
term  removes  the  infinity  and  gives  us  the  correct  result  when  Ay  =  0.  At  the 
same  time  it  is  negligible  for  all  off-axis  values  of  Ay  that  will  be  of  inter¬ 
est.  The  appropriate  expression  is  then 

2  222  2  2221  /2 
AP  II^W6yo(0)Jln[(F'^  +  k  a  Ay  )/(5y  /4  +  k  a^Ay  ) ]  '  .  (16) 

The  results  of  Eqs.  13  and  16  give  the  detected  scattered  power  from 
converging,  SAW  diffracted  rays  in  the  spectrum-analyzer  configuration  of  Fig. 
13.  The  same  formulas  may  be  used  to  calculate  the  scattered  power  from 
undiffracted  rays  once  appropriate  changes  for  the  terms  and  Ay  are  incor¬ 
porated.  1'  n  is  the  acoustooptic  diffraction  efficiency  and  is  the  dif¬ 
fracted  power  ler  unit  beam  width,  the  undiffracted  power  per  unit  beam  width  is 
Simil^rly,  if  ij;  is  the  angle  between  the  diffracted  and  undiffracted 
beams  (ij^  *  3®  in  Fig.  13).  and  Ay  is  the  value  of  detector  position  relative 
to  the  optical  axis  of  the  diffracted  beam,  then 

I  I  /n 

o  o 

Ay  Ftani})  -  Ay  ,  (17) 

are  the  substitutions  to  make  in  Eqs.  (13)  and  (16)  to  describe  the  detected 
scattered  power  from  the  undiffracted  beam.  For  the  geometry  of  Fig.  13,  the 
value  of  F  tanii>  is  943  urn,  while  the  largest  value  of  Ay  is  that  for  the  +50th 
detector,  400 
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♦  V  «  • 


& 


Contribution  to  Scattered  Power  from 
the  Unfocused  Beams 


Light  scattered  before  the  lens  at  an  angle  (p  is  focused  by  the  lens 
onto  the  detector  array  at  position  y*  =  Ftancf).  (As  before,  we  begin  by  consid¬ 
ering  only  the  contribution  of  the  SAW  diffracted  beam.)  The  scattering  angle 
subtended  by  each  detector  is  approximately  6y/F.  The  power  intercepted  from  a 
single  scattering  event  is  I^o((()) 6y/F.  Since  the  total  number  of  scattering 
centers  in  the  front  focal  region  is  EWF,  the  total  scattered  power  detected  is 

AP  =  n^WSyo((P)  (18) 

The  angle  dependence  of  a(cl))  is  that  of  Eq.(9),for  Rayleigh  scattering,  and 
Eq.  (10)  for  RGD  scattering. 

This  expression  for  scattered  power  from  the  SAW  side  of  the  lens 
applies  to  the  undiffracted  beam  when  is  replaced  by  l^/n  and  (p  is  replaced 
by  !]/-(().  From  Fig.  13,  the  on-axis  detector  receives  light  from  the  undif- 
frncted  beam  that  has  been  scattered  at  =  3° .  Consequently  the  <p  dependence 
of  Eq .  (18) is  not  significant  for  Rayleigh  scattering,  which  is  wide  angle, 
though  for  RGD  scattering  it  could  result  in  a  gradual  decrease  in  scattered 
power  across  the  detector  array. 

Total  Contribution  to  On-Axis  Scattered  Power 


The  total  scattered  power  is  obtained  by  summing  the  contributions 
from  the  diffracted  and  undiffracted  beams,  both  before  and  after  the  beams 
pass  through  the  lens.  At  first  glance,  one  might  anticipate  that  the  dominant 
contribution  comes  from  the  portion  of  the  undiffracted  beam  scattered  on  the 
detector-array  side  of  the  lens.  The  undiffracted  beam  is  more  intense  than  the 
diffracted  beam  by  a  factor  of  l/n,  while  scattering  on  the  array  side  of  the 
lens  is  detected  in  greater  amounts  because  of  the  larger  acceptance  angle  of 
the  detector  elements. 

For  this  contribution  the  power  intercepted  by  the  central  detector  is 
calculated  from  either  Eq .  (13) or  (16) , depending  on  the  type  of  scattering,  after 
making  tlie  substitutions  and  Ay  ->  Ftani|>-Ay  (=  Ftaniji  for  the  central 
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detector  element).  We  summarize  the  results  using  an  approximation  valid  when 
tani//<<l,  as  is  the  case  in  practice: 

AP  «  Z(I^/n)W6ya(0)  (19) 

9  2  1/2 

X  J2.n[  (1/Mtan\^)  (1  +  M^tan  \fr)  ]  , 

where  M  =  ka  for  the  RGD  limit  (ka»l)  and  M  =  1  for  the  Rayleigh  limit  (ka<<l). 
Taking  =  3*^,  k  =  27t(2 . 2) / .633  pm  =  21.8  pm  and  a  =  10  pm,  we  find  that  the 
logarithmic  term  above  has  the  value  2.95  in  the  Rayleigh  limit  and  0.0038  in 
the  RGD  limit. 

In  the  Rayleigh  limit,  then,  scattering  from  the  post-lens  region  is 
only  about  three  times  as  great  as  that  from  the  pre-lens  region,  where  the 
logarithmic  term  in  Eq,  (19) is  replaced  by  unity  in  accord  with  Eq.  18.  With 
n  -  10%,  the  total  on-axis  scattered  power  from  the  undiffracted  beam  in  the 
Rayleigh  limit  is  close  to 

AP  39.5  2:i^W5ya(0)  .  (20) 

The  numerical  factor  39.5  is  the  result  of  (2.95  +  1.0)/ti. 

In  the  RDG  limit,  the  value  of  the  logarithmic  term  in  Eq.  (19)  is  small 
enough  to  suggest  that  other  contributions  to  scattering  dominate  that  from  the 
undiffracted  beam  on  the  output  side  of  the  lens.  The  contribution  from  scat¬ 
tering  of  the  undiffracted  beam  on  the  input  side  of  the  lens  is  an  expression 

2  2  2  —1 

similar  to  Eq.  (19)  except  that  the  logarithmic  term  is  replaced  by  [1  +  k  a  ] 
(see  Eqs.  (10)  and  (18)).  In  the  numerical  example  chosen,  this  factor  is  found  to 
be  0.0076,  up  from  0,0038  but  still  small.  If  n  =  10%,  the  total  scattering 
contribution  from  the  undiffracted  wave  is 

AP  0.114  II^W6ya(0)  ,  (21) 

where  0.114  =  (0.0076  +  0.0038) /n. 

In  view  of  this  small  value,  let  us  now  consider  the  on-axls  scatter¬ 
ing  contribution  from  the  SAW-dif fracted  wave.  Taking  the  sum  of  Eqs.  (16) and  0-8) 
accounts  for  scattering  generated  both  before  and  after  the  lens.  For  the  cen¬ 
tral  detector,  we  have 
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The  quantity  in  brackets  is  9.88  for  the  parameters  of  Fig.  13.  This  is  two 
orders  of  magnitude  greater  than  the  scattering  contribution  from  the  undif-- 
fracted  wave.  This  imbalance  is  caused  by  the  assumed  large  size  of  the  scat¬ 
tering  centers,  a  =  10  pm,  which  results  in  a  highly  localized  scattered  energy 
distribution.  If  a  =  1  pm,  the  contributions  from  the  diffracted  and  undiffract¬ 
ed  beams  are  comparable.  Similarly,  if  n  =  1%  instead  of  10%,  the  relative  un¬ 
til  ffracted-beam  contribution  increases  on  order  of  magnitude  and  may  exceed  the 
dif fracted-beam  contribution. 

Note  that  the  result  of  Eq.  (22)  is  a  valid  expression  for  the  dif- 
fracted-beam  contribution  to  scattering  in  both  the  Rayleigh  and  RGD  limits. 
Consequently  it  may  be  added  to  Eq,  (20)  to  obtain  the  total  on-axis  scattered 
power  in  the  Rayleigh  limit  from  diffracted  and  undiffracted  beams,  both  before 
and  after  the  lens: 

AP(Rayleigh)  =  49,4  ET^W6ya(0).  (23) 

In  the  RGD  limit,  we  had  from  Eq.  (22) 

AP(RGD)  =  9.9  ri^J^iSya(O)  .  (24) 

The  numerical  coefficient  will  always  be  at  least  as  large  as  1  +  iln2F/6y,  con¬ 
tributed  from  the  diffracted  beam,  and  it  may  be  larger  if  the  angular  spread  of 
scattering  is  as  large  as  the  angle  ip  between  diffracted  and  undiffracted  beams. 
Generally  this  will  be  the  case  for  Rayleigh  scattering.  In  the  following  dis¬ 
cussion  we  will  write 

AP  =  QII^W6ya(0)  (25) 

where  Q  10  for  the  large -particle  RGD  limit  and  Q  ~  100  for  the  small- particle 
Rayleigh  limit. 
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Calculation  of  Dynamic  Range 


The  dynamic  range  of  a  spectrum  analyzer  is  the  amount  that  the  peak 
diffracted  signal  must  be  attenuated  to  equal  the  noise  associated  with  scatter¬ 
ed  light.  For  a  uniform  beam  of  width  W  and  power  per  unit  width  the  power 
per  unit  width  in  the  Fourier  transform  plane  is 

I(y')  =  Ip(W^/FX)  sinc^(y'W/FX),  (26) 

where  sincz  =  sinTTz/irz.  The  width  of  the  central  peak  of  the  Fourier  transform 
distribution  is  2FX/W.  For  a  spectrum  analyzer  having  the  geometry  of  Fig.  13 
F  =  18  mm  and  W  =  7,2  mm.  With  X  =  (0,633/2.2)  ym  =  0.29  ym  the 
width  of  the  central  peak  is  about  1.4  ym.  This  is  considerably  less  than  the 
5  ym  detector  aperture.  The  detected  power  is  therefore  essentially  the  total 
power  in  the  diffracted  but  unscattered  incident  beam: 


The  dynamic  range  is 


P  -  I  W 
o  o 


P  /AP  =  lQ)';6vc(0)] 
o 


The  factor  Q5y  in  this  expression  is  somewhat  under  control  of  the  design  engi¬ 
neer.  Provided  that  scattering  from  the  undiffracted  beam  is  negligible,  we 


Q>Sy  =  [1  +  Cn2F/.W](^y  .  (29) 

This  goes  to  zero  as  i^y  goes  to  zero.  However,  6y  =  5  ym  as  used  in  Fig.  13 
is  probably  close  to  the  state  of  the  art.  In  that  case,  improvement  in  dynamic 
range  is  afforded  by  keeping  F  as  small  as  possible.  F  =  18  mm  in  Fig.  13  is 
also  a  reasonable  small  value.  We  conclude  that  the  design  of  Fig,  13  is 
close  to  optimum  from  the  point  of  view  of  maximizing  system  dynamic  range. 


0¥' 


CALCULATION  OF  DIFFERENTIAL  SCATTERING 


CROSS  SECTIONS 


In  order  to  further  analyze  spectrum-analyzer  dynamic  range  it  is 
necessary  to  develop  expressions  for  the  on-axis  value  of  the  differential 
scattering  cross  section  a(0).  This  problem  was  addressed  in  Ref.  1  from  a 
rigorous  point  of  view,  though  ultimately  approximations  were  introduced  that 
limited  the  analysis  to  the  Rayleigh-Gans-Debye  limit;  that  is,  ka>>l  and 
k^dn<<l,  where  dn  is  the  effective  index  perturbation  of  the  scattering  parti¬ 
cle.  In  this  work  we  prefer  to  limit  the  rigor  of  the  analysis  but  to  pre¬ 
serve  its  physical  content  as  much  as  possible.  The  objective  is  to  produce 
simple  but  physically  sound  formulas  for  use  in  predicting  spectrum-analyzer 
performance  in  a  variety  of  scattering  situations.  In  addition,  we  hope  to 
extend  the  results  of  Ref.  1  to  include  the  Rayleigh  limit  of  scattering  by 
small  particles. 


Differential  Scattering  Cross  Sections  in 
the  Rayleigh  Limit 


Figure  15  shows  a  scJiematic  drawing,  4)f  dipole  scatters  distri¬ 
buted  throughout  a  waveguide  layer  of  depth  D.  If  we  ignore  the  effects 

(22) 

introduced  by  the  waveguide  geometry,  the  total  scattering  cross  section  is 

“tot  "  (8’f/3)k^a^(ni^-l)“/(m^+2)^  ,  (30) 

where  a  is  the  scattering  center  radius,  m  =  (n4'(Sn)/n,  and  n  is  the  refractive 

index.  The  ratio  of  the  terms  containing  ra  is  approximately  equal  to 
2  2 

(4/9)6n  /n  when  the  index  perturbation  satisfies  ^n<<n.  However,  this  is  a 
fair  approximation  even  when  the  scattering  centers  are  voids  in  the  waveguide 
layer,  having  n+6n  =  1.0  and,  for  LiiChOj,  n  =  2.2,  Therefore  we  use  the 
approximation  for  all  5n  of  interest,  writing 

“tot  ""  (8tt/3)  (4/9)k'''a'’6n^/n^  .  (31) 

The  differential  scattering  cross  section  is  reduced  by  a  factor  (87i/3)  and 

2 

multiplied  by  the  polarization  factor  cos  (j),  where  •'^/2-<J)  is  the  angle  between 


4] 


angle  of  the 


the  scattering  direction  and  the  incident  polarization*  In  the  problem  we  are 
considering,  that  of  in-plane  forward  scattering,  t})  =  0®  and  we  have 


^dif  ”  (^/9)k^a^6n^/n^  (32) 

as  the  relevant  differential  scattering  cross  section.  Note  that  a  ^  and 

^  tot 

^dif  have  dimensions  of  area,  whereas  a(0)  of  the  previous  subsections  had 

dimensions  of  length.  This  disparity  will  be  remedied  when  we  take  into  account 
the  reduced  dimensionality  of  the  waveguide  geometry. 

For  a  single  scattering  center,  the  amount  of  light  scattered  into  a 
solid  angle  A::  is 

AP  =  ,  (33) 

where  D  is  the  waveguide  depth  and  is  the  appropriate  power  per  unit  cross 

sectional  area  in  the  waveguide  layer.  We  arc  interested  in  the  solid  angle 
subtended  by  a  waveguide  detector.  Since  light  scattered  into  the  substrate 
is  presumed  not  to  be  detected,  we  need  consider  only  solid  angles  that  are 
hounded  by  the  acceptance  angle  of  the  waveguide.  This  angle,  AB,  is  defined 
l’>v  Fig.  1  and  is  found  using  Snell  *s  Law  Lo  he  given  by 

A0  ==  2(2An/n)^'^^  ,  (34) 

where  An  is  the  index  perturbation  of  the  waveguide,  assumed  small  in  compari¬ 
son  to  the  substrate  index  n.  If  the  detector  subtends  an  angle  A(()  measured 
in  the  plane  of  the  waveguide,  the  included  solid  angle  is  AQ  =  A0A<|),  or 

ASil  =  2 (2An/n)  ^  •  (35) 

The  detected  power  is  obtained  from  Eq.(fO: 

1/2 

AP  =  2f(I  /D)o  ._(2An/n)  A^f>  (36) 

o  all 

A  factor  f  has  been  added  to  account  for  the  fact  that  only  a  fraction  of  rays 
emitted  within  the  acceptance  angle  of  the  waveguide  is  trapped  in  waveguide 
modes . 

By  comparison  with  Eq.(l),the  in-plane  differential  cross  section 
to  be  used  in  the  previous  formulas  of  this  section  is 


a(0)  -  2f (0^^^/D) (2An/n) 


(37) 


1/2 

This  has  the  correct  dimension  of  length.  £q.  Q8) for  dynamic  range  depends  on 
the  product  Ea(0)  where  I  Is  the  number  of  scatters  per  unit  waveguide  area. 

If  p  Is  the  volume  density  of  scatters  In  the  waveguide,  then  I  -  pD  and 

Ea(0)  =  po^^^(2An/n)^''^  .  (38) 

We  have,  for  simplicity,  taken  f  =  1/2.  Eq*  (38)  can  be  cast  in  terms  of  the 
attenuation  coefficient  of  the  waveguide,  a  =  pa  ,  using  o  =  a  /(87t/3). 
We  find  for  the  spectrum  analyzer  dynamic  range, 

P^/*^P  =  [Q5y(3/8TT)a(2An/n)^^^]~^  (39) 

For  Q  =  100,  6y  =  5  um,  a  =  1  dB/cn,  An  =  0.007  and  n  =  2,2  we  find 

P^/AP (Rayleigh)  =  39,6  dB  .  (40) 

If  waveguide  losses  can  be  reduced  to  0.1  dB/cm,  the  dynamic  range  improves  to 
nearly  50  dB,  This  result  is  very  encouraging,  but  is  is  based  on  the  assump¬ 
tion  that  only  Rayleigh  scattering  exists  in  the  waveguide.  Other  scattering 
present  can  degrade  the  level  of  performance,  as  we  now  discuss. 

Differential  Scattering  Cross  Sections 
in  the  Rayleigh-Cans^-Debye  Limit 

We  consider  scattering  centers  which  are  much  larger  than  a  wave¬ 
length  (Mie  scattering),  yet  which  are  optically  soft  (RGD  limit).  This 
means  that  the  field  inside  the  scattering  center  is  essentially  the  same  as 
the  incident  field,  and  the  effect  of  the  scattering  center  is  simply  to 
perturb  the  wavefront  of  the  incident  light  by  an  amount  small  compared  to  a 
wavelength*  These  requirements  are  satisfied  when  ka>>l  and  k^a6n<<l,  where 
<Sn  is  the  index  perturbation  associated  with  the  scattering  center. 

Figure  16  defines  the  geometry  of  the  problem  we  are  considering. 

If  is  the  incident  wavefront,  the  wavefront  at  the  trailing 

edge  of  the  scattering  center  is 
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ns  analysis  In  the  Rayleigh-Debye  Limit. 


(41) 


z  +a/2 
o 

A  =  A  exp  { Ik  n  (z  -<l/2)+Tk  /dzn  (y,z)} 

Og  O  O'*  g 

z  -all 
o 


where  k  “  2Tr/A  , 
o  o* 

mode  index  in  the 
may  be  written 

A  =  A 


is  the  free-space  wavelength,  and  n^(y,2)  is  the  perturbed 

scattering  region.  If  we  take  n  (y,z)  ^  <5n  (y,z),  the  field 

g  o  6 

2  -Kl/2 
o 

e  {1  +  [exp(  ik^Jd26n^(y,z))  -  1]}  .  (42) 

2  -a/ 2 
o 


The  term  in  square  brackets  describes  the  perturbed  wavefront.  Its  variation 
along  the  y  axis  depends  on  the  shape  and  index  of  the  scattering  center  and 
in  general  will  be  quite  complicated.  To  keep  the  analysis  simple  and  yet 
preserve  its  important  features,  we  take  6n^(y,z)  to  be  a  rectangular  function, 

6n  =  constant  for  z  -a/2<z<z  'Ki/2 

o  o 


y^-a/2<y<y^+a/2  , 


6n^  =  0  elsewhere. 


(^3) 


We  also  assume  k  Ci6x\  <<1.  The  perturbation  field  is  then 

o  g 

A  =  A  e  ^^^^o  Ci6n  in  v  -a/2<v<v  '+'n/2, 

o  o  g  ■  o  u 


=  0  elsewhere- 


(44) 


The  scattered  field  in  the  far  zone  is  obtained  through  a  straightforward 

(23) 

application  of  Fraunhofer  diffraction  theory: 


K/  \  A  ik(z  •fa/2)  ..  .  r-\/  \ 

A(y,z)  -  A  e  o  a6n  [iA(2-z  )] 

*  o  o  o 

all 


Ikr  f  ^iky'((>^,,, 

-all 


e  /  '•'dy'. 


(45) 


The  parameters  of  the  problem  are  defined  with  respect  to  Fig.  17:  (y,z)  are 
the  coordinates  of  the  observation  point,  and  are  such  that  y''>y^,  z>>z^;  r  is 
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the  distance  of  the  observation  point  from  the  scattering  center  and  4^  is  the 
scattering  angle,  assumed  to  be  small.  The  integral  in  Eq.  (45)  has  the  value 
a  sinc(a(j)/A) .  Finally,  the  distance  r  is  given  by 

r  -  R-y^sln(t)-z^cos4>  ,  (46) 

where  R  is  the  distance  of  the  observation  point  from  the  origin  of  coordinates. 
Using  these  results  to  substitute  in  l-:q .  (4  5),  we  have,  in  polar  coordinates, 

A(R,()>)  =  A  a^6n  sinc(a(t./X) 

o  o  g 

^-iky^sin((>+ikz^(l-cos4>)  ,  (47) 

The  last  term  on  the  right  contains  phase  information  associated  with  the 
location  of  the  scattering  center.  The  other  terms  are  common  to  all  scattering 
centers  in  the  scattering  region.  We  may  replace  A(R,c())  by  the  field  A^(R,<()) 
for  the  ith  scattering  center  located  at  coordinates  (y^^jZ^^).  The  total 
scattered  field  is  then 

E(R,4>)  =  Z  A,  (R,<^)  (48) 

i  ^ 


and  the  scattered  intensity  is  proportional  to 


EE*  =  A  ^(k  iX^6n  )^(A2)  ^sinc^ (a^)/ \) S (4 )  , 
o  O 

SU)  =  |j;e-iky.sin^+ikz.(l-cosO|^ 


(49) 


The  summation  term  may  be  written 


s(o  =  N+E 


(50) 


where  N  is  the  total  number  of  scattering  centers.  The  maximum  value  is  S(0)  =  N 
However  the  terms  in  the  sum  will  begin  to  add  incoherently  at  angles  for  which 
the  maximum  values  of  k(y^-yj )sin<{)  or  k(z^-Zj )  (l~cos(J>)  approach  unity.  Since 
this  will  occur  at  angles  comparable  to  or  less  than  -  1/kW,  the  beam  diffrac¬ 
tion  angle,  the  coherent  contribution  to  scattering  cannot  be  resolved  experi¬ 
mentally.  At  larger  values  of  (p  than  1/kW^  the  exponents  of  some  terms  in  Eq.  50 
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are  large  and  the  value  of  S(((>)  quickly  falls  toward  an  average  value  of  N. 
Coherent  addition  could  occur  at  nonzero  values  of  <p  if  the  scattering  centers 
were  arrayed  in  a  periodic  lattice,  but  for  random  scattering  center  locations 
we  have  not  been  able  to  observe  even  partial  evidence  of  coherent  scattering 
when  the  sum  of  Eq«  50  is  generated  on  a  computer.  Nor  have  we  observed  any 
enhanced  scattering  at  nonzero  scattering  angles  in  our  experimental  work. 

This  does  not  mean  that  S((())  exhibits  no  significant  fluctuations  with 
()>.  Consider  a  short  scattering  path  length  L  so  that 

S((())  ^  N  4*  Z  2cos [k(y  .~y , ) sln<j>]  (51) 

i>j  ^  ^ 

We  have  cast  the  sum  in  terms  of  the  real-valued  cosine  function  and  thereby 

2 

reduced  the  number  of  terms  to  (N  -N)/2.  For  a  given  value  of  sin(()  -  4>>>l/kW, 
at  least  some  arguments  of  the  cosine  functions  in  Eq .  (57)  will  be  large  com¬ 
pared  to  unity.  If  we  assume  that  the  arguments  populate  the  interval  0<arg< 
(kWsin{|)>>l)  in  a  manner  which  is  uniform  and  random,  the  cosine  function  assigns 
to  each  argument  a  positive  or  negative  number  in  the  range  0  to  1  or  0  to  -1 
with  approximately  equal  probability.  We  simplify  the  problem  by  averaging  the 
values  of  the  cosine  function  in  these  ranges,  so  that  each  term  in  the  sum  is 

taken  to  have  the  value  +4/7t  or  This  reduces  the  problem  to  that  of  the 

2 

well-known  random  walk  in  one  dimension.  We  have  (N^-N)/2  terms  that  are 

valued  either  +4/7t  or  -4/7t  in  a  random  fashion,  such  as  might  be  determined  by 

flipping  a  coin.  The  average  value  of  such  an  ensemble  is  0.  Hence  <S(<|))>  =  N. 

However,  statistical  fluctuations  are  such  that  values  on  the  order  of  S((|))  = 

2  1/2 

Nt[N  “N)/2]  might  well  be  observed  at  typical  values  of  (p.  Thus,  S  fluc- 

-1  /2  -1 /2 

tuates  between  ~  N(l+2  )  and  S_  N(l-2  )  for  large  N.  The  ratio 

S^/S_  =  5.83  shows  that  fluctuations  in  the  scattered  signal  about  the  average 
could  easily  amount  to  7.7  dB. 

Let  us  define  NR((j))  to  be  the  randomly  fluctuating  component  of  S (({>), 

so  that 

SU)  =  N[l+R((l))]  .  (52) 

When  this  relation  is  substituted  into  Eq .  (49),  we  have 

EE*  =  A  ^N(k^<X^6n  )^(Az)  ^sinc^ (^4)/0  [  1+I1((|>)  ]  (53) 

o  o  g 
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Tho  dependence  of  scattering  is  contained  in  the  final  two  terms.  The 

2 

average  variation  of  the  scattered  intensity  is  contained  in  the  sine  (a^/X) 

term.  However  this  function  itself  has  structure  that  results  from  our 

unphysical  assumption  of  a  single  particle  size  a  for  scattering  centers,  when 

2 

a  distribution  of  sizes  is  more  realistic.  The  envelope  of  the  sine  function 

would  therefore  appear  to  be  a  more  physically  meaningful  quantity  in  Eq.  (53) 

The  envelope  has  a  value  of  unity  at  =  0  and  a  value  near  (-rra^/A)  at  large 

2  -1 

<t>.  A  Lorentzian  distribution,  [l+CTrac^/X)  ]  ,  has  both  these  functions.  In 

(3) 

fact,  the  derivation  of  Ref.  1,  based  on  the  approach  of  Marcuse  and  very 
different  in  spirit  from  the  derivation  here,  shows  a  Lorentzian  variation 
with  namely  [ l+(2'rrCl())/X)  ]  ,  Note  that  the  coefficient  of  <()  is  a  factor  of 

two  greater  than  that  obtained  using  our  heuristic  analysis.  The  discrepancy 
originates  in  the  statistical  treatment  used  in  Ret,  1  to  account  for  the  ran¬ 
dom  distribution  of  scattering  centers.  The  same  treatment  results  in  a  scat¬ 
tered  amplitude  twice  that  found  in  Eq ,  (47),  Incorporating  both  these  features 
to  preserve  consistency  with  Ref  1»  we  have 

EE*  =  4A^^N(k^a^6n^)^(Az)"^L(ka(t))[l+R(().)] 

L(ka(fi)  =  [l+(ka()))^]"^  .  (54) 

This  expression  is  readily  converted  into  experimental  intensity 

2 

parameters  using  the  relation  I  /I  =  EE*/A  .  If  the  scattered  light  is 

scat  o  o  ® 

measured  using  a  detector  of  width  5y,  the  detected  power  is  AP  =  I  6y  = 

2  scat 

I  EE*6y/A  •  From  Eq .  (54)  we  have 
o  o 

AP  =  4l^N(k^a^<5n^)^A“^(6y/z)L(ka(|>)[l+R(((>)]  (55) 

The  term  (Ay/z)  is  just  A(f),  the  angular  acceptance  of  the  detector,  while  N, 
the  number  of  scattering  centers,  is  "LW,  where  L  and  W  are  the  length  and  width 
of  the  scattering  area,  respectively.  Making  these  substitutions  in  Eq.(55), 
we  f  ind 

AP  =  ELW  I^  0(4))  A(j) 

a((|))  =  (2/ir)  ^nk^^a^5n^\(ka(ti)[l+R((ti)]  (56) 
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The  latter  expression  is  the  differential  scattering  cross  section  in  the 
Rayleigh-Gans-Debye  limit. 

Let  us  use  the  results  of  Eq.  (56) to  estimate  the  dynamic  range  of  a 

spectrum  analyzer.  The  formula  is  that  of  Eq.  (28),  with  Q  -  10  for  the  case  of 

-2 

RGD  scattering.  As  a  worst-case  situation  let  us  take  1  «  a  ;  that  is, 
scattering  centers  are  contiguous.  There  is  no  space  between  them.  Further, 
let  us  take  6y  =  5  pm,  a  =  b  =  10  pm,  n  =  2,2,  k^  =  27r/0.633  pm,  and  (j)  =  0  =  R(0). 
The  resulting  on-axis  dynamic  range  is 


-7  -2 

P  /AP  =  1.46  X  10  6n 
o 


(57) 


-5 


To  achieve  a  40  dB  dynamic  range  or  better,  it  is  necessary  that  6n  -  0.38  x  10 

-2  ^ 

This  is  to  be  compared  with  An  ^  10  for  the  index  change  at  the  surface  of 
LiNbO^  Ti-dif fused  waveguides. 


Calculation  of  6n 


Thus  far  we  have  made  no  assumption  about  the  origin  of  5n  .  Hence  we 

R 

cannot  comment  on  its  magnitude.  In  this  subsection  we  calculate  in  terms 
of  waveguide  and  scattering  parameters  for  several  physical  models  of  scattering 
beginning  with  surface  roughness. 

Surface-Roughness  Scattering.  We  model  the  effects  of  surface  roughness 
on  scattering  by  considering  the  rough  surface  to  be  associated  with  a  randomly 
fluctuation  waveguide  thickness  D.  If  n^(n)  is  the  effective- mode  index,  the  rms 
fluctuation  of  n^  associated  with  a  surface  having  rms  roughness  a  is 

6ng  =  (3ng/9D)a  •  (58) 

Our  problem  reduces  to  an  analysis  of  the  waveguide  dispersion  equation  to  obtain 
a  formula  for  3n^/3D  in  terms  of  k^.  An  =  n^-n2»  and  where  sur¬ 

face  (substrate)  refractive  index. 

Let  us  base  our  discussion  on  the  analysis  of  Hocker  and  Bums  for 
(24) 

They  derive  a  set  of  universal  curves  which 


diffused  optical  waveguides, 
relate  a  normalized  mode  index 


u  /  2  2.,,  2  2. 

b  »  (n^  -n2  )/(n^  -n2  ) 


(59) 
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to  a  normalized  waveguide  depth, 

V  =  .  (60) 

The  desired  quantity  is 

3n  /3D  =  On  /3b)  Ob/3V)  (BV/BD)  ,  (61) 

8  8 

The  second  term  on  the  right,  3b/3V,  is  just  the  slope  of  the  universal  curve, 
while  the  third  term  is 


3V/3D  =  k  k  (2n„An)^^^  . 

o  1  ^  o  Z 

The  first  term  is  obtained  by  inverting  the  derivative  of  Eq.  59: 


3n  /3b  =  (n  ^-n^^)/2n 
g  1  z  g 


An 


Consequently,  Eq.  61  reduces  to 


Bn  /3D  =  k  An(2n„An) 
g  o  2 


1/2 


(3b/9V) 


(62) 


(63) 


(64) 


To  complete  the  solution  to  the  problem,  one  works  from  a  knowledge  of  waveguide 
fabrication  procedures  to  select  the  appropriate  universal  curve  b(V)  and  the 
correct  operating  point  b  =  b(V^).  The  slope  of  the  curve  at  the  operating 
point  b’ (V^)  is  inserted  into  Eq,  (64)to  specify  the  derivative  Bn  /3D,  Using 
the  universal  curve  for  a  Gaussian  index  profile,  and  taking  V  =4.45,  the 


largest  value  possible  for  a  single  mode  waveguide,  we  have  9b/ 3V 


0.10.  With 
-7 


k  =  2tt/0.633  ym,  n„  =  2.2  and  An  =  0.007,  the  derivative  9n  /3D  is  1.23  x  10 

oZl  ^  ^ 

A  .  The  largest  value  of  the  derivative  for  a  single-mode  guide  occurs  with 

V  =  2.70,  and  is  only  2.6  times  the  above  value.  We  have  measured  the  slopes 

used  in  these  calculations  directly  from  a  set  of  universal  curves  provided  in  a 

preprint  of  Ref.  (7)  obtained  from  the  first  author.  It  appears  that  values  of 


-5 


can  be  obtained  with  surfaces  as  rough  as  100  A  rms. 


-7 


6n  on  the  order  of  10 

g 

Using  Eq.  57  for  dynamic  range  In  conjunction  with  our  result  =  al.23  x  10 

o  o  8 

A"i,  we  find  a  =  31  A  associated  with  a  dynamic  range  of  40  dB.  This, 

o 

however,  is  not  a  particularly  good  polish.  An  average  polish  having  a  =  10  A 

o 

is  consistent  with  a  dynamic  range  of  50  dB,  and  a  3.2  A  surface  roughness  is 
consistent  with  an  60  dB  dynamic  range.  These  results  offer  encouragement  with 
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regard  to  the  possibility  of  reducing  surface-roughness  scattering  to  negligibly 
important  levels. 

Before  moving  on,  we  present  a  derivation  of  6n  =  cr(3n  /3D)  for  a 
strongly  as5rmmetric  slab  waveguide.  For  strong  asymmetry  the  superstrate 
index  n^  is  not  close  to  n^  or  n2  in  magnitude.  As  a  result,  terms  dependent  on 
can  be  eliminated  from  the  model  dispersion  equaticn.  For  TE  modes,  the 
dispersion  equation  can  be  simplified  to 

V  =  .  (65) 

The  derivative  of  this  equation  with  respect  to  V  is 

1  =  (l/2)(db/dV)(l-b)"^[b“^'^“-(l-b)“^^^tan"^(b"^-l)^^^]  .  (66) 

The  last  item  in  brackets  can  be  replaced  by  its  equivalent,  V,  using  Eq.(65). 
Combining  this  result  with  Eqs.  (61),  (62),  (63),  and  (58),  we  obtain 

6n  =  k  a  2An(2n-An)^‘'^(l-b)(b)^''^ 

go  I 

x[k^D(2n2Anb)^'^^+l]‘^  (67) 

One  may  wish  to  insert  this  result  directly  into  Eq.(56)to  obtain  an  expression 
for  differential  scattering  cross  section  (;(4))  in  terms  of  waveguide  and  rough¬ 
ness  parameters  alone.  Such  an  expression  was  derived  in  Ref.  1  using  a  more 
rigorous  approach.  If  the  two  agree,  it  adds  credence  to  the  present  results. 

We  begin  with  Eq.0‘^0of  Ref.  1  and  insert  the  results  of  Eqs.  (9)  ,(23)  ,  and  (20)  of 
that  reference  to  obtain 

dP/d((>  =  P(2iT/3o)(k  ^/16TT^)(n,^-n  ^)^a^E  ^(0)(3  /2u)p  P)^ 

O  .1.  o  o  o  o 

x4Na^H(4))  (68) 

In  our  current  notation,  we  have  replaced  P  by  I^  on  the  right  hand  side  above, 
where  I^  is  power  per  unit  waveguide  width.  This  is  done  to  avoid  confusion 
with  dP,  which  was  and  is  an  increment  of  power.  Also  B  =  n,,k  ~  n^k  is 

the  magnitude  of  the  wavevector  of  the  waveguide  mode,  is  the  vacuum 
permeability,  oj  is  2ttv  ,  v  is  the  optical  frequency,  ^^^(0)  is  the  waveguide  field 
at  the  surface,  and  N  =  )!LW  is  the  number  of  scattering  centers.  Our  current 
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definition  of  the  differential  scattering  cross  section  is 

a(il>)  -  (I^SLW)"^(dP/d((i)  (69) 

This  is  different  from  that  of  Ref.  1  by  the  inclusion  of  the  scattering  center 

density  T  in  the  denominator.  Eq.  (68) gives  us  the  information  necessary  to 

write  an  expression  for  a((j)).  However,  it  is  useful  to  first  express  the  field 

at  the  surface  E  (0)  in  terms  of  physical  waveguide  parameters  .  ^  For  the 

o 

strongly  asymmetric  waveguide  we  are  considering,  the  field  is  approximately 
given  by 

E^^(O)  =  4(l-b)(nj^^-n^^)"^(up^I^/6^)k^b^^^(2n2An)^'^^ 

x(Vb^^^+l)"^  .  (70) 

Using  Eq.  (67),  this  result  may  be  rewritten  as 

E  ^(0)  =  (wp  I  /6  )2(n.^-n  ^)  (6n  /aAn)(2n,An)  (71) 

O  OOO-LO  g  ^ 

Upon  inserting  this  expression  into  Eq .  (68)  and  then  inserting  Eq .  68  into 
(69) ,  we  find 

a(0)  =  (2/'ir)n2k^^<x‘^6n^^  (72) 

This  agrees  with  our  heuristic  result  of  Eq ,  56.  We  conclude  that  the  deriva¬ 
tions  of  this  section  are  theoretically  sound.  In  particular,  the  critical 
assumption  that  the  effect  of  large  scattering  centers  can  be  expressed  in 
terms  of  the  mode  index  perturbation  6n^  appears  correct. 

Volume  Scattering .  In  this  subsection  we  consider  first  the  model 
for  large  volume  scattering  centers  that  was  introduced  in  Ref.  1.  The  scat¬ 
tering  centers  are  assumed  to  have  an  index  different  from  their  surroundings 
by  an  amount  5n. 


This  model  for  volume  scattering  might  be  expected  to  describe  fluc¬ 
tuations  in  the  Ti  concentration  or  in  the  stoichiometry  of  the  host  LiNbO^ 
material.  One  source  of  scattering  that  we  have  considered  is  that  of  gaps 
in  the  waveguide  layer  caused  by  dust  particles  trapped  beneath  the  Ti  film 
prior  to  diffusion, At  the  high  temperatures  of  diffusion,  these  may  explode 
and  effect  the  removal  of  Ti  from  the  surface,  leaving  a  void  in  the  waveguide 
Ti  distribution. 

In  this  case  n^  =  n^  in  the  vicinity  of  the  scattering  center;  that 
is,  there  is  no  waveguide  (we  ignore  outdif fusion) .  We  have 

By  combining  Eq,  56  with  Eq,  28,  and  taking  Q=10,  n=2,2,  k^= 
27t/0.633  pm,  and  ~  bAn,  we  find  that  the  on-axis  dynamic  range  is  given  by 

Pp/AP  =1.46  X  10"^  um^/a^L- b^An^  (74) 

We  will  take  a  -  1  pm  as  the  largest  dust  particles  which  are  difficult  to 

locate  on  a  surface  and  difficult  to  remove  because  of  electrostatic  forces. 

Similarly,  we  will  take  bAn  0.003  as  the  largest  value  of  the  index 

inhomogeneity  associated  with  Ti  gaps  in  the  waveguide.  Then  for  a  40  dB 

-2 

dynamic  range  we  find  a  density  of  scattering  centers  of  5]  =  162  mm 


I 


i 


This  is  a  rather  dirty  surface,  probably  not  encountered  in  careful 

practice.  If  I  is  reduced  two  orders  of  magnitude  to  a  more  reasonable  level  of 

-2  -2 
1.6  mm  ,  the  dynamic  range  is  increased  to  60  dB.  If  I  =  0.2  mm  ,  while 

(1^3  um,  the  dynamic  range  is  50  dB#  We  conclude  that  scattering  from  Ti  gaps 

in  the  waveguide  is  not  likely  to  impair  spectrum-analyzer  performance  provided 

care  is  taken  in  waveguide  preparation. 

Other  volume  scattering  centers  discussed  in  Sec.  Ill  had  to  do  with 

the  occurrence  of  the  separated  phase  LiNb^Oo  and  Li-Ti-0  compounds  like 

(1A,7)  3  8 

Li2Ti^0^  during  waveguide  fabrication.  The  index  of  refraction  associated  with 
these  species  is  likely  to  differ  by  6n  =  0,1  or  more  from  that  of  the  host 
species.  This  is  generally  too  large  to  qualify  as  ROD  scattering,  and  will  be 
treated  later  as  Mie  scattering. 

However,  if  the  scattering  centers  are  confined  near  the  waveguide 
surface  to  within  a  depth  the  effective  index  is  likely  to  be  reduced  from 

6n  because  of  the  poor  overlap  with  the  waveguide  electric  field. 

From  the  discussion  following  Eqs.(2)  and  (6)  of  Ref.  1,  the  scattering 
parameter  of  that  reference  is  calculated  from 

Zo 

n 

rio  =  (k^n2/2a)Vi^I^)2n26n  /  dzE^  (z)  (75) 

o 

where  is  the  depth  of  the  scattering  inhomogeneity  and  Eq(z)  is  the  normalized 

electric  field.  We  have  been  considering  =  a>>D,  the  waveguide  depth.  In 

this  case  =  2n26n  results  from  replacing  by  in  Eq, (75),  For  the  case 

we  now  consider,  z  <<D  and 
o 

rig  =  (zo/D)(^n2An)(n^^-n^^)“^2n2i5n[Vb^^^/(Vb^'^^+l)]  (76) 

In  deriving  this  result,  we  have  approximated  the  integral  in  Eq.  (75)  as 
2  2 

z^E^  (0),  where  E^  (0)  is  given  by  Eq.(70).  For  cases  of  interest,  V  ~  3  and 
b  =  0.2,  the  quantity  in  square  brackets  above  is  about  1/2  in  value.  Use  of  this 
simplification  along  with  n2  =  2.2  and  An  =  0.007  gives  us 

n  =  8.0  X  10  ^(z  /D)2n^6n  .  (77) 

o  o  z 


56 
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In  other  words  the  scattering  perturbation  is  reduced  by  a  factor  8  x  10 

(z  /D)  from  the  case  z  »D.  We  can  write  this  result  as 
o  o 


6n 

g 


8.0  X  10  ^(z^/D)6n, 


z  <<D 
o 


(78) 


Taking  6n  =  0.1,  z^  =  0.3  ym  and  D  =  3.0  yra,  we  have  6ng  *  8  x  10  Inserting 

this  result  into  Eq.  (57),  we  predict  a  dynamic  range  of  only  14  dB.  This  re¬ 
sult  assumes  that  we  are  dealing  with  contiguous  scattering  centers  of  lO-ym 
size.  The  dynamic  range  increases  to  24  dB  for  3.2  ym  scattering  centers.  Note 
that  (Sn  should  be  taken  to  describe  fluctuations  in  the  high  index  layer  at  the 
surface,  since  a  uniform  high  index  layer  would  produce  no  scattering.  Since 
this  layer  is  presumed  to  result  from  compounds  having  6n  0.1,  we  should 
really  insert  the  variation  6(6n)  into  Eq .  (57).  The  value  5(<5n)  =  0.01  would 
increase  the  dynamic  range  by  20  dB  to  the  vicinity  of  40  dB,  depending 
on  scattering  center  size. 

Nevertheless,  this  scattering  mechanism  has  to  be  considered  important. 
Experimentally  we  have  verified  its  presence  by  measuring  improvement  in  our 
waveguides  upon  polishing  the  surface  following  the  diffusion  process. 

In  order  to  complete  our  study  of  RGB  volume  scattering  centers,  let 
us  now  consider  the  possibility  of  scattering  from  Ti  concentration  fluctuations 
tliat  occur  as  a  result  of  the  statistical  nature  of  the  Ti-film  deposition  and 
diffusion  process.  These  cannot  be  avoided  regardless  of  care  taken  in  wave- 

o 

guide  preparation.  A  200  A  Ti  film  diffused  to  make  a  2  ym  deep  waveguide  will 

20  “3 

have  an  average  concentration  in  the  material  of  c  =  5.7  x  10  cm  .  If  we 

2 

consider  a  surface  area  of  Ci  ,  the  average  number  of  Ti  atoms  within  that  area 
2 

is  <N>  -  where  D  is  the  diffusion  depth.  The  standard  deviation  in  that 

1/2  2  1/2 

number  is  about  <N>  =  (ca  D)  .  The  corresponding  fractional  variation  in 

1/9  2  -1/2 

the  number  of  Ti  atoms  in  each  cell  is  <N>  ‘'/<N>  =  (c^l  D)  .  Accordingly, 

if  An  is  the  index  change  associated  with  the  average  Ti  count,  then  6n  = 

2  ,-1/2 

Af  (cd  D)  is  the  index  change  associated  with  statistical  fluctuations  in  the 

2 

Ti  count  in  each  scattering  area  This  quantity  is  to  he  inserted  into  the 

expression  for  scattering  cross  section,  Eq.  (56). 

Although  it  appears  rather  unphysical  that  (Sn  depends  on  the  arbitrarily 
chosen  size  of  the  scattering  area,  Ci,  this  dependence  cancels  out  of  the  product 
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Ea(<t>)  upon  which  all  measured  quantities  depend*  This  occurs  because 
A  2  "^2 

o(()))  «  a  6n  and  I  *  a  for  adjacent  scattering  areas  of  size  a.  We  find 
from  Eq,  (56) 

la(0)  =  (2/ir)  n  An^/cD  .  (79) 

When  this  result  is  inserted  into  Eq*  28  for  dynamic  range,  we  find  P^/AP^  ® 
79.3  dB.  In  performing  the  calculation  we  have  utilizied  Q  =  10,  6y  =  5  iJm, 
n  =  2.2,  =  27r/0.633  ym,  An  .007,  and  T)  =  2  ym.  The  result  shows  that 

scattering  from  statistical  fluctuations  in  the  Ti  concentration  is  not  likely 
to  be  a  factor  limiting  spectrum-analyzer  dynamic  range. 

Calculations  of  Waveguide  Attenuation.  Since  we  have  calculated  the 
spectrum-analyzer  dynamic  range  associated  with  surface  roughness  and  volume¬ 
scattering  centers  in  the  RGD  limit,  it  is  interesting  to  estimate  what  those 
mechanisms  imply  regarding  waveguide  attenuation.  In  Eq .  (AO)  of  Ref.  1  we 

present  expressions  for  the  attenuation  coefficient  associated  with  surface 

2 

roughness  and  volume  scattering.  (A  factor  An  has  been  inadvertently  omitted 

from  the  expression  for  volume  scattering.)  Moreover,  these  expressions  are 

2 

proportional  to  E^(0)  ,  the  square  of  the  waveguide  field  at  the  surface.  The 

formula  for  E  (0)^  given  in  Eq.  (24)  of  Ref,  1  should  be  replaced  by  Eq .  70  of 

^  1/21/2 
this  report.  The  result  is  the  addition  of  a  factor  [Vb  / (Vb  +l)](l-b)  to 

the  formulas  for  the  attenuation  coefficient.  This  factor  shows  the  dependence 

on  waveguide  modal  characteristics,  but  for  cases  of  interest,  V  -  3,  b  ~  0.2, 

the  added  factor  is  close  to  1/2  in  value.  Incorporating  this  factor  into  the 

results  of  Ref.  1,  Eq.  (40),  we  have 


=  (2/3TT^)n,^'^^An  E  a^k  ^o^/D(k 
I  o  o 


“vol  (2/37)^)n2^'’^(n,^-n^^)~^An5:a^6n^/D(k^a)^''^  .  (80) 

-2  -2 

For  surface  roughness,  the  parameters  S  =  10  pm  ,  a  =  10  pm, 

o 

An  =  0.007,  k^  =  2tt/0.633  pm,  D  =  2.55  pm,  n2  =  2.2,  and  a  =  31  A  \jere^  found  to 
provide  a  40  dB  spectrum^analyzer  d5mamic  range.  When  inserted  into  the  first 
of  Eq.  (80),  the  same  parameters  provide 
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(81) 


a  -  “  1.1  X  10  ^dB/cm 
surf 

For  volume  scattering  centers,  parameters  for  40  dB  dynamic  range 
include  I  =  6  nun  ,  6n  =  0.003,  and  a  =  3  ym.  If  in  addition  An  =  0.007, 

1^2  ~  2.2,  and  D  =  2.55  ym,  we  find  from  the  second  of  Eqs.  (80) 

a  ,  =  1.0  X  10  ^^dB/cm  (82) 

vol 

The  tremendous  reduction  relative  to  the  attenuation  from  surface  roughness 

3  3 

comes  primarily  from  a  10  reduction  in  I  and  a  10  reduction  in  the  value  of 
2  2 

6n  relative  to  the  value  of  •  The  result  emphasizes  a  conclusion 

reported  in  Ref.  1:  that  volume  scattering  centers  are  efficient  generators  of 
in-plane  scattered  energy,  while  surface  roughr^^ss  scattering  centers  are 
efficient  generators  of  out-of-plane  scattered  energy.  However  neither  mecha¬ 
nism  produces  much  scattering  of  either  type  if  the  calculations  of  this  section 

(26) 

can  be  relied  on.  Reported  levels  of  attenuation  in  LiNbO^  waveguides  may 
well  reflect  intrinsic  absorption  losses  and  scattering  from  other  sources  than 
the  RGD  scatters  considered  here. 


Differential  Scattering  Cross  Sections  in 
the  Mie  Limit 


In  Section  II  we  described  several  types  of  scattering  centers  that 
do  not  appear  to  be  well  described  by  Rayleigh  scattering,  because  of  their 
potentially  large  size,  or  by  Rayleigh-Gans-Debye  scattering,  because  of  their 
large  refractive  index  difference  from  their  environment.  We  will  use  a  simple 
but  attractive  analysis  to  predict  the  scattering  cross  section  for  these 
Mie  scattering  centers.  We  assume  them  to  be  large  in  comparison  to  a  wave¬ 
length,  so  that  a  beam  of  light  passing  through  them  is  subject  to  Snell’s  law 

of  refraction.  If  the  typical  radius  of  curvature  seen  by  an  incident  wave- 

(2  7) 

guided  beam  is  a,  an  approximate  application  of  the  lens-makers  formula^ 
suggests  that  the  output  beam  will  focus  at  a  distance 

R  a/26n  (83) 

At  distances  much  greater  than  R,  the  half-angle  of  divergence  of  the  beam  is 
given  by 
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(84) 


0,.  =  tan  ^a/R  tan’’^2(Sn 

div 

In  the  waveguide  geometry,  most  of  the  light  emitted  into  the  upper  part  of  the 

light  cone  is  reflected  from  the  waveguide  surface  into  the  substrate,  doubling 

the  intensity  in  that  region.  At  a  distance  r>>R,a  from  the  scattering  center, 

the  area  subtended  by  the  light  cone  in  the  substrate  is  a  semicircle  of  radius 

2  2 

r  tan  ~  2r6n.  The  area  of  the  semicircle  is  A  =  2iTr  5n  .  If  the  power 

scattered  by  the  inhomogeneity  is  equal  to  the  power  intercepted, 

‘Ptot  ’ 

the  power  per  unit  area  in  the  scattered  beam  is 

AP^  ^/A  =  I  a/7Tr^5n^  .  (86) 

tot  o 

2 

The  effective  area  covered  by  a  waveguide  detector  at  a  distance  r  is  A^  =  r  Afi 
where  AJ^  is  the  solid  angle  subtended  by  the  detector,  given  in  terms  of  the 
waveguide  acceptance  angle  by  Eq.  (35),  We  find 

=  2r^(2An/n)^^^A<ti  (87) 

where  A«^  is  the  angle  subtended  by  the  detector  in  the  plane  of  the  waveguide. 

The  power  intercepted  by  the  detector  is 

1  /2  2 

AP  =  AP  A, /A  =  (l/7i)(2An/n)  ^  (I  af/6n‘')A<t)  .  (88) 

tot  d  o 

As  before,  a  factor  f  is  introduced  to  account  for  the  fact  that  not  all  rays 
scattered  within  the  acceptance  angle  will  be  trapped  in  waveguide  modes:  f  is 
the  fraction  trapped.  The  differential  scattering  cross  section  is 

o((i>)  =  AP/T^A<J)  =  (1/ti)  (2An/n)^^^rtf/(5n^  .  (89) 

Actually  light  is  scattered  only  in  the  range  “^^^iv  t)y  the 

light  cone.  Outside  this  range  0(4»)  =  0.  For  (Sn  =  0.1  the  range  is  20^^^  =  22.6*^. 
This  is  intermediate  between  Rayleigh  scattering  and  Rayleigh-Cans-Debye  scatter¬ 
ing.  However,  the  angle  Is  large  enough  relative  to  the  scattering  angles 
important  in  the  spectrum-analyzer  applic<ition  that  the  problem  may  be  treated 
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as  we  treated  Rayleigh  scattering;  that  is,  the  dynamic-range  equation  is 


P^/AP  =  [Q5;o(0)(Sy]“^  ,  (90) 

where  the  geometric  factor  Q  ^  100  as  for  Rayleigh  scattering.  Also  taking 
.Sv  =  5  urn,  and  inserting  Eq.  (89)  into  Eq.  (90),  we  find 

P^/AP  =  [(500  pm)  (I/tt)  (2An/n) ^^^Eaf /5n^]  ^  .  (91) 


This  formula  may  he  expressed  in  terms  of  the  waveguide  attenuation  coefficient. 
It  we  ignore  the  power  trapped  in  waveguide  modes  in  comparison  to  the  total 
[)owcr  scattered,  the  power  lost  to  the  waveguide  over  a  path  length  L  is 


AP,  =  21  ar.LW 
lost  o 


(92) 


The  attenuation  coefficient  is 


Substitution  of  this  result  into  Eq.  (91)  vields 


-4  2 

P  /AP  =  [1.46  X  10  af/6n  ] 
o 


(94) 


where  we  have  used  An  =  0.007  and  n  =  2.2  to  obtain  the  numerical  coefficient, 
and  where  is  given  in  dB/cm, 

The  largest  values  of  (Sn  that  we  mi gh t  be  expected  to  encounter 

experimental  Iv  are  (Sn  =  2. 2-1.0  =  1.2,  for  pits  in  the  waveguide  surface  or 

subsurface  voids,  and  (Sn  ~  0.15  for  the  separated  phase  LiNb^Co.^^^Li-Ti-O 

J  o 

compounds  formed  during  the  diffusion  process  are  also  likely  to  have  tSn^O.l. 

Smaller  values  of  are  treated  using  R(ID  theory. 

Taking  a  =  0.1  dB/cm,  f  “  1/2  and  Sn  -  0.12  in  Fq .  (94)  wo  find 

Pq/APo  =  dB.  This  is  among  the  smaller  dvnannc'  ranges  we  have  calculated 

vet.  The  condition  a  =  0.1  dB/cm  could  result  from  scattering  centers  as  small 

-2 

as  I  urn  and  densities  as  low  as  1  mm  Such  small,  dispersed  scattering 

centers  might  he  difficult  to  locate  despite  tlieir  large  index  changes.  Thus 
it  is  well  to  consider  tliem  as  relevant  to  the  spect rum-analyzer  dynamic  range 
problem. 


CALCULATIONS  OF  SCATTERING  AT  90 


i 


I 


Thus  far  we  have  emphasized  forward  seat  ter  iriK  because  it  is  this 
component  that  will  limit  spectrum-anal vzer  dynamic  range.  However,  forward 
scattering  is  difficult  to  study  experimentally  because  of  the  presence  of 
the  unscattered  beam  and  of  scattered  light  from  optical  components  of  the 
beam  forming  system.  The  observation  of  scattt'ring  at  90°  is  useful  as  a 
means  i>f  c  i  rcumvt‘nting  thest^  ctimpK-x  i  t  i  es  .  However,  it  must  be  remembered 
that  of  all  the  types  of  scattering  wi'  have  discussed  only  Rayleigh  scattering 
i>f  TM  nunU's  will  havt‘  a  lari;e  d  i  f  f  e  r^Mil  i  a  1  cross  sei't  ion  feu'  f  =  90°. 

T.et  us  ctmsider  the  expt'riment  shown  in  Fig.  18  in  whicli  in-plane 
scatteri’d  light  at  <f=90°  is  coupled  out  of  the  waveguide  and  viewed  by  a  dark- 
adapted  eye.  The  power  into  the  eye  is 

AP  =  ,  (95) 

where  P  is  the  power  in  the  waveguide,  L  is  the  length  of  the  beam  path  c'ontri- 
buting  to  the  observal  Li>n,  A^J'  is  the  angle  subtended  by  the  dilated  pupils  of 
the  ohse  ers  eyes,  and  n  is  the  output  coupl  ing  t'fficioncy  of  the  prism 
locati'd  between  tiie  scattering  area  and  the  observer. 

We  will  Lake  AP  -  10  W  to  be  the  threshold  of  vision.  This  number 

is  based  t>n  pfiys  i(^  logical  studies  aiul  its  tlerivation  is  presented  in  Appendix  A. 
TluMi  we  can  suhstituto  api^ropr  iate  values  for  P,  L,  n,  and  and  determine  the 
smallest  value  of  Xo(tt/2)  that  will  pr(^cliice  a  visual  effect.  The  experimental 
p.ir.imetors  nre  P  =  2  X  W,  I.  =  1  cm,  ti  s  0.1  and 

=  d/nn  =  7.0  X  lO"^  ,  (96) 

where  d  =  5mm  is  the  pupil  aperture,  D  =  ^0  cm  is  the  approximate  viewing 
distance,  and  2.2  is  the  refractive  correction  necessary  to  account  for  the 

Snel 1 ’ s-1 aw-assoe iated  increase  in  the  scattering  angle  upon  leaving  the  wave¬ 
guide  and  entering  free  space.  These  values  indicate 

7.0  (-ri/ 2)  •  6.6  X  10  um  ^  (97) 
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for  a  visual  effect. 


For  Rayleigh  scattering  of  TM  waveguide  modes,  Za(7r/2)  =  Zo  (3/87r), 

tot 

where  ^  is  the  total  scattering  cross  section  and  Zo  is  the  attenuation 
tot  tot 

coefficient,  a.  We  thus  find 

a  ^  2,4  X  10~^  dB/cra  =  0.024  dB/kra  (98) 

in  order  for  scattering  from  dipoJe  sources  to  he  observed  at  90°.  This  is 
two  orders  of  magnitude  better  than  the  current  attenuation  levels  found  in  the 
best  optical  fibers. 

The  reader  may  have  already  anticipated  that  in  the  experiment  we 

performed,  no  light  was  detected.  Even  if  our  crude  calculations  arc  in  error 

by  two  orders  of  magnitude,  the  result  suggests  that  a  negligible  level  of 

dipole  scattering  is  to  be  found  in  LiNbO^  waveguides.  If  the  attenuation 

coefficient  of  Eq .  (98)  is  inserted  into  Eq.  (39),  for  example,  we  predict  a 

spectrum  analyzer  dynamic  range  in  excess  of  110  dH. 

Let  us  now  consider  RCD  and  Mie  scattering  centers,  which  are  large 

in  comparison  to  wavelength.  Our  derivations  of  scattering  cross  sections  for 

these  types  of  scattering  centers  have  made  inso  of  their  refractive  properties. 

This  has  led  to  the  result  that  scattering  is  concentrated  in  the  forward  or 

near  forward  direction,  with  a(TT/2)  -  0.  i.et  us  now  consider  the  scattering 

contribution  caused  by  reflection  from  these  scattering  centers.  Wo  take  the 

total  scattered  power  as  equal  to  the  intercc'pted  power  times  an  average 

reflection  coefficient.  Tht'  later  quantity  depends  on  the  interior  index  of 

the  scattering  center,  on  its  shape  and  on  the  optical  polarization.  For 

2 

normal  incidence  the  reflection  coefficient  is  [^n/(2n  +  6n)  j  and  for  grazing 

2 

incidence  it  approaches  unity.  We  will  take  (Sn  as  an  average  value.  We 

further  assume  that  the  reflected  power  is  distributed  uniformly  over  a  solid 

angle  of  The  solid  angle  encompassed  by  the  detector  is  A0A({),  where 

1/2 

AO  =  2(2An/n)  is  the  acceptance  angle  of  the  waveguide.  If  a  fraction  f  of 
light  reflected  within  the  acceptance  angle  of  the  waveguide  is  trapped  in 
waveguide  modes,  the  detected  power  is 

AP  =  (2I^^aiSn^)  ( 1 /An)  2  (2An/n)  ^  ^^A(}>f  .  (99) 
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The  differential  scattering  cross  section  for  light  scattered  in  plane 


It  an  angle  4)  =  90®  is 


2  1/2 

o((f))  =  (AP/I^A<(>)  =  a6n  f(2An/n)  /tt 


(100) 


Inserting  this  result  in  Eq.  (97),  first  taking  f  ==  1/2,  An  -  0.007,  and  n 


2.2,  we  find  that 


Za6n^  >  0.52  X 


(101) 


must  be  satisfied  to  produce  visible  Mic  scattering  at  90®. 

For  the  case  of  weak  RGD  scatterers  confined  near  the  waveguide 

surface  we  assume  that  most  of  the  reflected  radiation  remains  trapped  in  the 

waveguide  mode.  This  results  in  an  expression  for  scattering  cross  section 

1/2 

similar  to  that  of  Eq,  (TOO),  except  that  the  product  f(2An/n)  does  not 
appear ,  The  corresponding  inequality 


2  -11  -1 
Za  6ng  ^  0.2  X  10  pm 


(102) 


must  then  be  satisfied  to  produce  a  visible  scattering  effect  at  90®. 

Consider  first  a  dense  system  of  RGD  scattering  centers,  such  as 

-2  ° 

surface  roughness,  described  by  Z  =  a  ,  a  =  10  pm  and  o  =  30  A.  This  is 

found  to  produce  a  dynamic  range  of  40  dB  and  a  value  n  =  3.7  x  10  .  We 

9  1 1  1  ^ 
find  Zafin*"  =  0.14  x  10*  pm  .  Considering  the  limitations  of  the  calcula- 
g 

tion  the  observation  of  visible  scattering  at  90®  is  marginally  possible.  For 

the  case  of  LiNb^Oo  or  Li-Ti-0  compounds  near  the  waveguide  surface,  the 
7  ° 

parameters  Z=  a”  ,  a  =  10  pm,  (z  /D)  =  0.1,  and  6(6n)  =  0.01  produce  34  dB 
dynamic  range  and  Za6n^  =  0.6  x  10  pm  .  A  visible  scattering  effect  is 
again  marginally  indicated. 

Next  consider  scattering  that  might  result  from  gaps  in  the  waveguide 

caused  by  dust  trapped  beneath  the  Ti  film  prior  to  diffusion.  We  found  that 

7  “3 

a  =3.2  pm,  Z  =  15  mm*  ,  and  6n  =10  produced  a  spectrum-analyzer  dynamic 

S  2 

range  of  40  dB.  The  same  parameters  are  associated  with  the  product  Za6n^  = 

4.8  X  10*^^  pm“^.  Again,  a  visible  scattering  effect  at  90®  is  marginally 

indicated . 


Next  consider  the  scattering  situation  where  a  =  10  ym  and  6n  ~  0.1, 

such  as  might  be  associated  with  separated  LiNbO^  phases  or  lithium  titanaf.e 

compounds.  The  scattering  centers  are  not  assumed  to  be  confined  near  the 

2 

waveguide  surface.  A  density  of  E  =  0.1  mm  is  found  to  lead  to  a  spectrum 
analyzer  dynamic  range  of  about  34  dB.  For  the  same  parameters  we  find 

=  lO"®  ura"^  .  (103) 

This  is  200  times  greater  than  the  value  from  Eq.  (101)  required  for  a  visible 
scattering  effect  at  90®.  If  Ea  is  reduced  accordingly  to  provide  consistency 
witii  our  null  scattering  observation  at  90®,  the  spectrum-analyzer  dynamic 
range  is  predicted  to  increase  to  47  dB.  However,  we  are  now  dealing  with 
scattering  center  densities  so  low  that  individual  scatterers  can  be  avoided 
in  experiments  by  careful  selection  of  optical  beam  paths.  This  is  not  found 
to  be  the  case  in  practice.  We  conclude  that  our  null  observation  of  scatter¬ 
ing  at  90°  is  an  argument  against  the  existence  of  Mie  scattering  by  islands 

associated  with  LiNb  0  and  Li-Ti-0  compounds. 

3  o 

Finally,  consider  scattering  at  90°  by  large  polishing  imperfections, 

such  as  pits  in  the  surface  as  subsurface  voids,  both  having  6n  -  1.  If  a  =  1 
-2 

and  E  =  1  mm  ,  corresponding  to  a  =  0.086  dB/cm  and  a  spectrum -analyzer 
dynamic  range  of  45  dB,  we  have 


1360^"  =  10  ^  ym  ^  .  (104) 

This  is  20,000  times  the  value  required  for  a  visible  scattering  effect  at  90°. 
This  is  a  strong  indication  that  this  type  of  scattering  is  not  important  in 
practice.  We  conclude  that  our  null  scattering  observation  is  most  consistent 
with  ROD  scattering  from  surface  roughness,  surface  compounds,  and  gaps  in  the 
waveguide  Ti  distribution.  Rayleigh  and  Mie  scattering  mechanisms  are  least 
consistent  with  our  observation. 

SUMMARY 

In  this  section  we  derived  equations  for  spectrum-analyzer  dynamic 
range  in  terms  of  spectrum-analyzer  design,  scattering-center  density,  and 
scattering-center  cross  section.  Then  we  considered  various  models  for 
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scattering  centers  and  calculated  cross  sections  based  on  those  models.  This 
has  enabled  us  to  predict  which  scattering  mechanisms  are  likely  to  be  of 
greatest  significance  in  limiting  spec t rum ^analyzer  dynamic  range.  Our  results 
are  dispersed  throughout  the  section,  so  it  is  appropriate  to  summarize  them 
here.  This  is  done  in  Table  I. 

Only  the  entries  in  the  last  column  of  Table  I  may  require  explana¬ 
tion.  These  are  ratios  of  calculated  power  scattered  at  90°  to  the  estimated 
minimum  detectable  power  for  a  human  observer,  10  W.  Values  of  this  ratio 

greater  than  unity  indicate  that  the  scattered  light  should  be  sufficiently 
intense  to  produce  a  visual  effect.  Values  less  than  unity  are  consistent  with 
our  opposite  experimental  result.  However,  the  calculation  is  so  approximate 
and  the  experiment  so  qualitative  that  we  suggest  that  ratios  as  large  as 
several  hundred  may  be  compatible  with  our  observation. 

On  this  basis,  all  scattering  mechanisms  listed  in  Table  I  are 
candidate  mechanisms  except  for  Rayleigh  scattering  from  defects  smaller  than 
a  wavelength  and  Mie  scattering  from  surface  or  subsurface  polishing  imperfec¬ 
tions.  Mie  scattering  from  islands  of  LiNb^O^  and  Li-Ti-0  compounds  that  are 
large  or  comparable  to  waveguide  thickness  are  possible  limiters  of  dynamic 
range  to  values  less  than  40  dB;  however,  the  size  and  index  of  these  islands 
would  make  it  likely  that  they  could  be  seen  with  a  microscope.  Since  we  have 
not  done  so ,  we  discount  their  importance . 

This  leaves  Rayleigh-Gans-Debye  (RGD)  scattering.  Reasonable  values 
of  surface  roughness  show  a  rather  high  dynamic  range  associated  with  this 
source  of  RGD  scattering.  Comparable  or  higher  dynamic  range  is  associated 
with  Ti  voids  in  the  waveguide  layer  and  statistical  fluctuations  in  the  Ti 
concentration.  The  only  remaining  candidate  for  significant  levels  of  RGD 
scattering  is  the  formation  of  LiNb^O^  and  Li-Ti-0  scattering  centers  near  the 
waveguide  surface.  In  fact,  this  source  of  scattering  is  verified  experimen¬ 
tally.  Burns  et  al  report  evidence  for  the  formation  of  Li-Ti-0  compounds  dur¬ 
ing  waveguide  fabrication. They  also  report  results  indicating  the  possible 
confinement  of  these  compounds  near  the  top  0.3  pm  of  the  waveguide  layer.  The 
etch-resistance  that  we  have  found  in  diffused  layers  may  be  a  further  indica¬ 
tion  of  the  presence  of  these  compounds.  Additionally  we  have  found  reduction 
of  scattering  when  the  waveguide  is  polished  following  diffusion . 
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The  results  of  this  section  are  consistent  with  the  experimental 
findings,  and  suggest  that  once  this  source  of  scattering  is  eliminated,  no 
other  of  the  potential  sources  considered  are  likely  to  limit  spectrum  analyzer 
dynamic  range  to  values  less  than  about  40  dB* 
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V.  MEASUREMENT  AND  EVALUATION  OF  WAVEGUIDE  PERFORMANCE 


THE  IN-PLANE-SCATTERING  EXPERIMENTAL  METHOD 


The  measurement  of  in-plane  scattered  energy  distributions  was 
performed  in  this  program  much  as  in  the  preceding  program  and  described  in 
Ref.  1.  The  experimental  setup  is  shown  in  Fig. 19.  A  Gaussian  beam  of 
approximately  1-mm  waist  is  prism  coupled  into  and  out  of  the  waveguide  and 
focused  with  a  lens  of  about  10-cm  focal  length  on  a  slit  of  variable  width 
s.  Typically  s  is  in  the  range  20-60  ym.  \^en  the  slit  is  displaced  from 
the  optical  axis  by  a  distance  x,  it  collects  light  scattered  in  the  wave¬ 
guide  at  an  angle  (j^  =  x/n  f,  where  n  is  the  mode  index. 

The  scattered  power  through  the  slit  is  typically  20-40  dB  down 

from  the  peak  power.  Moreover,  the  peak  power  must  be  kept  small  to  avoid  the 

(99) 

generation  of  optical  damage  by  the  photoref ractive  effect.  As  a  result, 

we  have  the  problem  of  detection  of  low  light  levels*  For  this  reason  we  use 
a  photomultiplier  tube  to  monitor  the  light  passed  by  the  slit.  In  addition, 
light  coupled  into  the  waveguide  is  amplitude  modulated  by  periodically 
varying  tlie  angle  of  incidence  of  the  input  beam,  using  an  oscillating  mirror. 
The  modulated  component  of  the  photomultiplier  tube  output  is  measured  using 
a  lock-in  amplifier.  The  log  of  the  output  from  the  lock-in  amplifier  is 
generated  by  a  log  amplifier  and  delivered  to  the  y-axis  of  an  x-y  recorder. 
The  x-axis  displacement  is  proportional  to  the  separation  x  of  the  slit  from 
the  optical  axis  of  the  beam.  When  the  x  and  y  axes  are  properly  calibrated 
we  obtain  a  plot  of  log^^CAP  (th) /P)  versus  ()),  where  AP((|))  is  the  power  through 
the  slit  from  rays  scattered  at  an  angle  ({>  in  the  waveguide  and  P  is  the  peak 
power  through  the  slit  from  the  unscattered  beam. 

Dynamic  Range 

Owing  primarily  to  the  need  to  keep  the  peak  power  in  the  waveguide 
low,  the  dynamic  range  of  our  measurement  was  limited  to  50  dB.  That  is,  the 
electronic  noise  level  was  about  53  dB  below  the  signal  generated  when  the 
maximum  power  was  passed  through  the  slit.  The  dynamic  range  could  be 
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increased  by  a  number  of  techniques:  cooling  the  photomultiplier  or  replacing 
it  with  a  lower-iioise  tube,  increasing  the  lock-in  time  constant,  increasing 
the  slit  width  at  low  power  levels,  working  at  optical  wavelengths  or  at 
elevated  temperatures  where  optical  damage  is  not  a  problem  and  higher  beam 
powers  may  be  used.  We  actually  investigated  the  simplest  of  these  methods, 
increasing  the  time  constant  and  varying  the  slit  width.  However,  both 
methods  entail  a  certain  amount  of  inconvenience.  A  longer  time  constant 
requires  a  longer  data-acquisition  time.  Varying  the  slit  width  requires  a 
more  complicated  data  analysis.  In  fact,  spectrum-analyzer  performance  is 
specified  by  the  on-axis  scattering  and  the  angular  range  of  scattering.  Since 
both  of  these  are  measurable  when  the  dynamic  range  is  50  dB,  we  made  no 
unusual  attempts  to  increase  system  dynamic  range  further.  Ultimately  this 
should  be  done  in  order  to  measure  the  scattering  at  angles  beyond  a  few 
degrees  from  the  optical  beam  axis. 

Sources  of  Spurious  Scattering 

A  perhaps  more  serious  experimental  problem  than  limited  dynamic 
range  is  the  optical  noise  generated  by  light  scattering  at  components  other 
than  the  waveguide  under  test.  Each  intersection  of  the  beam  with  a  surface 
is  a  potential  source  of  scattering.  This  includes,  in  Fig. 19,  the  Fourier- 
transfoirm  lens,  the  mirrors  used  to  direct  the  beam,  and  the  coupling  prisms. 
It  also  includes  items  left  out  of  Fig. 19  for  simplicity,  such  as  the  laser 
and  all  beam-forming  and  polarizing  optics. 

Scattering  from  items  located  well  in  advance  of  the  waveguide  can 
be  reduced  by  placing  an  aperture  slightly  larger  than  the  beam  diameter  in 
front  of  the  oscillating  mirror  in  Fig. 19.  Scattering  from  mirrors  and  the 
lens  close  to  the  waveguide  can  at  least  be  measured  to  demonstrate  that  it 
doesn't  contribute  much  to  the  total  scattering  level.  However,  we  still 
have  scattering  from  the  input  and  output  coupling  prisms  which  cannot  readily 
be  separated  from  scattering  in  the  waveguide. 
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Prism  Enhanced  Scattering 

Our  best  attempt  at  doing  this  has  been  to  measure  the  total 
scattering  from  the  waveguide-prism  system  with  the  prisms  about  1-2  mm 
apart,  and  then  with  the  prisms  about  15  mm  apart.  If  prism  scattering  is 
negligible,  the  total  scattering  should  increase  in  proportion  to  waveguide 
path  length.  Accordingly,  a  10  dB  variation  is  anticipated.  In  fact,  the 
variation  is  in  the  range  0-3  dB,  leading  to  the  conclusion  that  scattering 
at  the  prism  interfaces  is  an  important  part  of  the  total  scattering 
observed . 

Subsequent  to  the  acquisition  of  data  for  this  program,  wc  under¬ 
took  a  theoretical  study  of  scattering  associated  with  prism  coupling.  The 
results  of  this  study  are  contained  in  Appendix  C,  We  found  that  even  a 
perfect  prism  can  serve  to  enhance  the  amount  of  scattering  from  waveguide 
scattering  centers  beneath  the  prism.  The  source  of  the  enhanced  scattering 
is  the  evanescent  field  associated  with  light  from  the  input  beam  that  strikes 
the  base  of  the  coupling  prism  and  is  totally  reflected.  This  evanescent 
prism  field  can  be  more  intense  than  the  waveguide  field,  causing  the  level 
of  scattering  from  the  prism-coupling  region  to  rival  scattering  from  the 
free  waveguide  despite  the  much  longer  path  length  of  the  free  waveguide. 

If  this  path  length  is  designated  as  Lr> ,  the  effect  of  prism- 
enhanced  scattering  may  be  characterized  by  a  parameter  R  such  that 

L  =  Lg  +  RL^  (105) 

is  the  effective  path  length  over  which  scattering  may  be  said  to  occur.  If 
is  the  prism  coupling  length  and  there  is  no  prism  enhancement  of  scat¬ 
tering,  we  expect 

^  ”  ^min  “  ^c/^"g  (106) 

Ordinarily  “  1/10.  However,  the  results  of  Appendix  C  show  that 

R  1  is  possible,  while  our  experimental  results,  to  be  described,  suggest 
that  R  -  1. 

In  any  case  the  best  approach  to  data  obtained  using  input- and 
output- pr ism  coupling  is  to  subtract  the  scattered  signal  obtained  when  the 
free  path  length  iias  two  different  values,  say  L.  and  L^.  The  difference 
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signal  then  represents  scattering  from  the  path  length  L^-Lg.  Assuming 
that  R  remains  the  same  in  the  tvo  measurements,  the  effects  of  prism 
enhanced  scattering  are  subtracted  away. 

Polished  Edges  for  End  Fire  Coupling 

There  must  be  some  concern  as  to  whether  the  above  cited  proced  ire 
works  xcell,  in  view  of  the  possibility  that  R  varies  with  each  new  placement 
of  the  prism  on  the  waveguide.  This  concern  is  justified  by  occasional 
experimental  results  that  show  greater  scattering  when  the  prism  separation 
is  reduced. 

To  avoid  this  source  of  experimental  error,  we  thorouglily  investi¬ 
gated  tile  possibility  of  replacing  prism  coupling  with  end-f ire  coupl ing 
tiirough  polisiied  waveguide  edges,  Tiie  procedures  used  to  generate  polisfied 
edges  were  as  follows. 

Two  LiNbO'j  slabs  obtained  from  Crystal  Technology  were  clamped 
tightly  together  between  two  steel  2*1. 4- mm  cubes  that  were  ground  to  an 
angular  precision  of  +  30”  (angle  between  faces  ==  90°  4*  30”).  The  LiNhO^ 

slabs  were  25mm  x  2  5mm  x  3mm,  \^ith  both  25mm  x  25mm  faces  polisiied  to  a 
flatness  of  \/10.  The  samples  were  cleaned  to  permit  intimate  contact  when 

clamped  between  the  stainless  steel  Mocks. 

Tile  entire  assembly  was  ground  using  1  5-pin  alumina  abrasive  on  a 
glass  lap  until  the  faces  t(''  be  polisiied  were  flusii.  At  tills  stage  the 
25mm  x  3xm  LiNbO^  faces  iiave  a  ligiitly  fros'td  appearance.  Polishing  is 
accomplished  through  the  subsequent  use  of  increasingly  smaller  abrasives, 
including  3  pm  diamond,  1  pm  diamond  on  solder  laps,  and  finally  colloidal 
silicii  in  tlie  form  of  Syton  lIT-30.  All  operations  arc  performed  in  an  auto-* 
mated  fnsliion  using  a  Logitech  PM2  polishing  machine.  For  the  final  polisii 
with  Syton,  the  solder  lap  is  replaced  by  an  expanded  polyurethane  lap. 

In  nearly  all  samples  fabricated  in  this  manner,  there  exist 
millimeter-length  regions  along  the  polished  edge  which  appear  to  be  defect- 
free  under  SEM  examination.  Mot  'over,  curvature  of  the  edge  is  not  apparent 
even  undei  10,000X  magnification,  indicating  that  the  radius  of  curvature  of 
the  edge  is  well  under  0.1  pm.  These  characteristics  bode  well  for  the 
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replacement  of  prism  input  coupling  with  end-fire  coupling.  However,  cer¬ 
tain  problems  are  inherent  in  the  technique. 

One  serious  problem  results  from  the  intimate  contact  achieved 
between  the  two  LiNbO^  slabs  butted  together  for  polishing.  When  polishing 
is  completed,  the  two  slabs  are  held  together  by  strong  attractive  forces 
and  cannot  be  separated  without  introducing  at  least  a  small  degree  of 
translat ion«'il  motion  of  one  sample  across  the  face  of  the  other.  This  usually 
results  in  some  damage  to  the  two  faces  in  contact,  which  are  those  on  which 
the  waveguide  is  to  be  fabricated.  The  damage  often  appears  to  be  caused  by 
an  abrasive  particle  or  a  LiNb03  chip  that  is  introduced  to  the  interface 
region  during  the  grinding  operation. 

This  problem  is  of  tolerable  severity  when  only  one  polished  edge 
is  required  per  waveguide.  One  can  usually  select  sample  areas  of  good  sur¬ 
face  quality  for  study.  The  problems  associated  with  end-fire  coupling  are 
componded,  however,  when  one  tries  to  obtain  two  polished  edges  for  both  in¬ 
put  and  output  end-fire  coupling.  We  expended  significant  energy  in  this 
direction  owing  to  our  desire  to  mimic  the  spectrum-analyzer  configuration  as 
closely  as  possible. 

A  second  polished  edge  per  sample  requires  greater  effort  than  the 
first  polished  edge.  Ideally,  the  assembly  consisting  of  steel  cubes  con¬ 
tacting  LiNbO'j  slabs  should  simply  be  inverted  in  order  to  polish  the 
opposite  3mm  x  25mm  crystal  faces  using  the  identical  procedure  outlined 
above.  In  practice  adjacent  LiNb03  slabs  will  not  be  in  intimate  contact 
over  their  entire  area  owing  to  a  slight  bow  in  one  or  both  crystals.  The 
resulting  polished  edge  is  rounded  to  an  unsatisfactory  degree  and  the  crystals 
have  to  be  remounted  for  successful  polishing  to  occur.  Additionally,  the 
bow  affords  an  entry  way  for  particulate  matter  that  can  damage  the  surface 
when  the  crystals  are  separated  prior  to  remounting,  and  then  after 
remounting  and  polishing  of  the  second  edge. 

We  attempted  to  surmount  these  problems  hv  using  a  nonabrasive 
cement  between  the  two  adjacent  LiNb03  surfaces.  This  would  serve  to  protect 
the  surfaces.  At  the  same  time,  it  was  hoped  that  the  cement  would  be  suf¬ 
ficiently  thin  and  close  to  LlNb03  in  hardness  tliat  little  edge  rounding  would 
occur. 
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The  cements  that  we  employed  were  Lakeside  cement  and  dental  resin* 
The  former  polished  away  too  rapidly,  causing  severe  edge  rounding.  The 
dental  resin  produced  a  satisfactory  edge  sharpness  except  in  those  locations 
where  a  bubble  had  formed  during  setting*  In  these  regions,  chipping  or  dis¬ 
tortion  of  the  polished  edge  took  place*  Possibly  handling  procedures  could 
be  developed  for  producing  highly  uniform  layers  of  dental  resin,  but  we  did 
not  pursue  the  matter  further. 

To  those  who  may  wish  to  develop  and  later  experiment  with  end-fire 
input  and  output  coupling,  we  note  two  additional  problems  with  polished 
edges  that  can  occur.  Each  edge  may  have  regions  of  good  quality  extending 
for  several  millimeters;  however,  these  regions  may  not  both  be  in  line  with 
a  beam  travelling  in  the  preferred  direction.  In  that  case,  one  may  have  to 
bear  with  scattering  from  one  of  the  two  edges.  A  further  problem  with  end- 
fire  coupling  is  the  separation  of  substrate  light  from  waveguide  light. 

Prisms  accomplish  this  quite  nicely,  a  fact  that  helps  to  compensate  for  the 
disadvantage  of  prism-enhanced  scattering. 

In  our  opinion,  the  optimum  experimental  configuration  for  studying 
in-plane  scattering  makes  use  of  end-fire  input  coupling  and  prism  output 
coupling.  A  single  polished  edge  can  be  generated  without  too  much  difficulty 
or  damage  to  the  waveguide  surface.  This  allows  one  to  avoid  prism-enhanced 
scattering  by  use  of  end-fire  input  coupling.  The  use  of  a  prism  output 
coupler  permits  good  discrimination  between  waveguide  and  substrate  modes. 
While  prism  enhanced  scattering  does  occur  during  output  coupling  owing  to 
the  evanescent  field  of  the  output  coupled  beam,  the  enhancement  from  this 
source  of  scattering  is  low  enough  to  be  ignored.  Much  of  our  discussion  of 
experimental  data  will  refer  to  an  in-plane  scattered  energy  distribution  that 
was  obtained  using  end-fire- input  and  prism-output  coupling. 


NEThODS  FOR  ANALYSIS  OF  EXPERIMENTAL  RESULTS 


Development  of  the  Basic  Formulas 

In  Sec.  IV  we  derived  expressions  for  the  scattered  power  de¬ 
tected  by  diode-array  elements  in  a  spectrum  analyzer*  These  expressions 
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do  not  relate  directly  to  the  experiment  shown  in  Fig,  19  because  of  signi^ 
ficant  variations  from  the  geometry  of  that  experiment  to  the  spectrum 
analyzer  geometry  shown  in  Fig,  13,  Therefore  it  is  necessary  to  derive 
new  expressions  which  will  enable  us  to  relate  the  detected  signal  in  our 
in-plane  scattering  experiments  to  the  waveguide-scattering  parameters  dis¬ 
cussed  in  Sec.  IV. 

The  starting  point  is  the  equation 

AP  =  EWLIorr(4))  Acf*  ,  (107) 

where  AP  is  the  power  scattered  in-plane  at  angle  c})  to  the  optical  path,  and 
within  an  angular  sector  A(j),  measured  in  the  waveguide  material,  Z  is  tiie 
density  of  scattering  centers,  Iq  is  the  power  per  unit  beam  width,  W  is  the 
beam  width,  a((>)  is  the  differential  scattering  cross  section,  and  L  is  the 
effective  path  length  in  the  waveguide,  which  may  include  the  contribution 
from  prism  enhancement,  if  appropriate,  as  described  by  Eq,  105.  l^en  this 
power  is  coupled  out  of  the  waveguide,  and  Aif)  change  as  required  by  the 
Snell’s  law  of  refraction.  For  the  small  angles  of  interest,  we  have 

^ext  '  ’  (108) 

=  np(A4')  , 

measured  external  to  the  waveguide,  where  n  is  the  mode  index.  In  the  focal 
plane  of  the  Fourier  transform  lens  in  Fig,  19,  the  light  energy  scattered  at 
in  the  waveguide  appears  at  a  distance 

X  =  f'f’ext  =  (109) 

from  the  optical  axis.  It  spans  an  interval 

Ax  =  f(A({')ext  =  fn^A.|)  (110) 

where  f  is  the  lens  focal  length,  A  slit  of  width  s  in  the  focal  plane  there¬ 
fore  collects  light  scattered  over  an  angular  range 

H  =  s/n  f  (111 

g 

in  the  waveguide.  Inserting  this  result  into  Eq,(107),  we  obtain 

AP  =  KWl.Ioa((t))s/n^f  (112) 

as  the  detected  power  at  the  slit  location  x  - 


This  equation  is  not  valid  on-axis,  where  the  slit  collects  light 
from  the  unscattered  incident  beam.  Let  us  consider  a  beam  having  a  peak 
power  per  unit  width  of  and  a  total  power  of  I^W.  If  the  intensity  pro¬ 
file  of  the  beam  is  uniform,  W  is  the  actual  beam  width.  If  the  profile  is 
Gaussian,  W  =  W^(tt/8)^^^,  where  is  the  diameter  of  the  beam  measured 
between  the  points  at  which  the  intensity  is  I^e"*^.  If  the  slit  width  is 
small  in  comparison  to  fX^/W,  the  power  passed  by  the  slit  is 

AP(0)  =  P  =  qlo(w2/fXp)s  ,  (113) 

where  q  =  1  for  a  uniform  beam  and  q  ~  2  for  a  Gaussian  beam.  If  the  slit 
width  is  larger  than  fX^/W,  the  power  passed  by  the  slit  is  the  total  power 
Ij^W,  Through  the  use  of  the  log  amplifier  shown  in  Fig.  19,  the  quantity 
measured  in  the  experiment  and  plotted  on  the  x-y  recorder  is  the  negative 
o  f  the  dynamic  range 

-D.R.  =  10  log^^lAP ((}))/?]  .  (lU) 

For  a  small  slit  width,  we  obtain 

AP((j))/?  =  .  (115) 

Th"s  is  the  smallest  value  of  AP/P  that  is  consistent  with  a  given  value  of 
)'a(4)).  As  the  slit  width  increases,  AP((|))  increases  linearly  with  s,  while  P 
saturates  above  s  =  fX^/W.  When  this  happens, 

APCct*)/?  =  ^La((|))  (s/n^f)  .  (116) 

Since  in  our  experiments  s  -  fX^/W,  it  appears  more  appropriate  to  use  this 
result  than  Eq.  115.  The  answer  is  the  same,  however  if  Eq .  115  is  used  with 
q  =  1,  as  for  a  uniform  beam: 

AP((j.)/P  =  E(L/W)a(<^)(Xo/n^)  (117) 

Taking  W  =  1  mm,  L  =  15  mm,  n^  =  2.2  and  Xq  -  0.633  ym,  we  obtain 

AP(<^)/P  =  4.32  pm  5:a(<^)  (118) 

Recall  that  the  spectrum-analyzer  dynamic  range  was  described  by  a  formula 

AP(4>)/P  =  [Q6yla(4.)l"’  (119) 


78 


where  10  -  Q  ^  100  depending  on  the  spectrum-analyzer  geometry  and  scattering 
angular  range,  and  6y  -  5  ym  was  the  detector  aperture.  This  indicates  that 
the  spectrum- analyzer  dynamic  range  will  be  between  10.6  and  20.6  dB  lower  (worse) 
than  the  dynamic  range  measured  using  the  system  of  Fig.  19.  If  we  were  to 
evaluate  a  waveguide  sufficient  to  produce  a  spectrum  analyzer  dynamic  range 
of  40  dB,  we  would  want  the  measurement  system  of  Fig. 19  to  possess  a  dynamic 
range  of  about  60  dB.  Moreover,  scattering  from  the  various  optical  com¬ 
ponents  that  intersect  the  beam  should  have  their  scattered  intensity  reduced 
at  least  70  dB  in  order  not  to  rival  the  scattering  from  the  waveguide.  The 
best,  and  perhaps  only  way  to  achieve  this  is  to  eliminate  all  optical  com¬ 
ponents  on  the  detector  side  of  the  waveguide,  and  to  carefully  spatially 
filter  the  light  incident  on  the  waveguide  to  remove  as  much  of  the  scattered 
component  as  possible. 

We  did  not  do  this  during  the  program  because  the  scattering  from 
the  waveguide-prism  combination  was  in  the  range  20-30  dB,  and  was  much 
greater  than  the  background  scattering  from  components.  Since  the  observed 
scattering  falls  off  by  10  dB  within  1®  of  the  optical  axis,  it  appears  to  be 
Rayleigh  Cans  Debye  scattering  from  scattering  centers  of  about  10  pm  diameter. 
Thus  Q  -  10  is  appropriate  in  Eq. 119,  and  a  30  dB  dynamic  range  in  the  experi¬ 
ment  of  Fig. 19  translates  into  a  20  dB  dynamic  range  for  the  spectrum  analyzer. 
This  is  20  dB  lower  than  desired  and  compatible  with  the  results  of  Table  I 
of  Sec.  IV  only  for  the  case  of  scattering  by  LiNb30g,  Li-Ti-0  compound  for¬ 
mation  near  the  waveguide  surface. 


Mechanics  of  Data  Analysis 


The  result  of  performing  the  in-plane  scattering  experiment  shown 
in  Fig.  19  is  a  recorder  tracing  of  10  log^Q[AP(<j))/P]  versus  (j).  Figure  20 
shows  the  experimental  result  that  will  be  the  basis  for  much  of  our  later 
discussion.  It  is  desired  to  analyze  data  of  this  form  to  obtain  information 
about  scattering  parameters  of  the  waveguide.  According  to  Eq.  117, 

A?U)f?  =  [Ka(<l))](L/W)(Xo/n^) 


(120) 


0 


10  log  (AP/P) 
(dB) 


Fig.  20.  In  plane  scattered  energy  distribution  for  sample 
154  (TMq  mode,  kj^c  obtained  using  end-fire  input 
coupling.  Scattering  angles  are  those  which  would 
be  measured  in  the  waveguide  material. 
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Known  parameters  are  W,  and  n^,  L  may  or  may  not  be  known  accurately. 

If  prism  input  coupling  is  used, 

L  =  Lg  +  Le  (121) 

where  Lg  is  the  known  separation  between  coupling  prisms  and  is  the  path 
length  associated  with  prism  enhancement,  which  may  be  comparable  to  Lg. 

The  result  of  Fig.  20  was  obtained  using  end-fire  input  coupling, 
so  that  Lg  =  0  and  Eq.  14  applies  with  the  product  Icr(<)>)  being  the  only 
unknown  quantity.  The  rapid  fall-off  of  the  scattered  intensity  with 
increasing  (p  suggests  that  the  observed  scattering  is  of  the  Rayleigh-Gans- 
Debye  type  discussed  in  Sec.  IV.  The  formula  for  the  differential  scattering 
cross  section  given  in  Eq,  56  is 

a((fi)  =  (2/TT)nkoV<5ng^[l+(kona()))^]  [1+R(<(>)]  (122) 

In  this  expression  R((f^)  is  a  random  function  that  is  to  be  associated  with 
the  fluctuations  in  scattered  intensity  that  are  observed  in  Fig. 20.  The 
average  curve  drawn  through  the  data  is  that  corresponding  to  <R ((}))>  =  0, 
and  therefore  has  a  Lorentzian  fall-off  with  <f>.  By  inserting  the  average 
value  of  a(()))  into  Eq.  120,  we  obtain 

AP(<(i)/P  =4(L/N)k-^5:a^6n  ^[l+(k  na(}>)^]  (123) 

^  g  o 

The  theoretical  curve  in  Fig.  20  is  consistent  with  the  parameter  values 

2  2 

L  =  15  mm,  W  =  1  mm,  k  =  27t/0.633  pm,  n  =  2.2,  a  =  11.7  ym  and  la  6n  = 

-10  ° ^ 

2.25  X  10  .  The  on -axis  dynamic  range  of  37.4  dB  of  the  theoretical  curve 

suggests  that  the  dynamic  range  of  a  spectrum  analyzer  fabricated  using  this 

waveguide  would  be  37.4  -  10,6  =  26.8  dB.  This  is  in  accord  with  Eqs.  118,  119 

and  the  discussion  following  those  equations. 

These  conclusions  are  based  on  the  fitting  of  a  Lorentzian  curve 

to  the  data.  We  have  generally  done  this  in  the  following  manner.  The 

quantity  AP/P  is  measured  at  2N  different  scattering  angles 

(J>n,  =  ±mA(l),  m  =  1,2,...N  .  (124) 

2 

From  Eq.  123  the  quantity  P/AP  is  a  linear  function  of  (p  ,  having  slope  = 

(k^na) ^A<L/W)ko^?^a^6n  ^  and  intercept  =5  Pf(L/W)k  2?^a^5n  The  fit  to  the 
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experimental  curve  is  obtained  by  plotting  the  2N  data  points  vs 

and  determining  the  best  straight  line  using  the  method  of  least*squares 
deviation. 

Typically  At))  =  0,08®  and  2N  10.  Points  are  not  usually  taken 
for  angles  larger  than  (t>  =  NA(|)  ~  0,4®  because  the  fluctuations  in  scattered 
energy  with  (()  become  severe  and  the  contribution  from  electronic  noise 
becomes  a  consideration.  Often  the  fluctuations  in  scattered  energy  are  even 
severe  at  smaller  angles,  causing  the  data  points  to  deviate  significantly 
from  the  best- fitting  straight  line.  However,  when  one  determines  the 
standard  deviation  in  slope  and  intercept  and  plots  the  ’^highest"  and 
"lowest**  Lorentzian  fit  to  the  experimental  distribution  using  this  infor¬ 
mation,  one  finds  that  the  curves  tend  to  lie  along  the  top  and  bottom  enve¬ 
lopes  determined  by  the  local  maxima  and  minima  of  the  fluctuations.  This 
situation  would  continue  to  exist  even  if  A(()  were  decreased  and  N  increased, 
because  the  fluctuations  correspond  to  uncertainties  that  are  inherent  in 
the  experimental  data  and  not  the  statistical  analysis  of  the  data. 


Complications  Associated  with  Prism-Enhanced  Scattering 


This  situation  is  more  complicated  when  prism-input  coupling  is 
used,  since  the  effective  scattering  path  length  L  is  unknown,  owing  to  the 
contribution  from  prism-enhanced  scattering.  The  unknown  quantities  in  Eq,  120 
are  now  ra(()))  and  L,  and  two  measurements  are  needed  to  determine  these  quan¬ 
tities  separately. 

As  we  have  indicated  previously,  the  two  measurements  involve  the 
determination  of  the  function  AP(()))/P  at  two  different  waveguide  prism  sepa¬ 
rations,  say  Lg  and  Lg.  The  effective  path  lengths  L  and  L'  used  in  Eq.  14 
are 

L  =  Lg  +  Le  (125) 


where  Lg  and  Le  are  the  prism  enhancement  path  lengths  appropriate  to  the 
two  measurements.  The  difference  in  scattered  power  obtained  by  subtracting 


the  two  results  is 


[AP((t))-AP'((^)]/P  =  [^a((t))](Ao/Wng)a^~L’g+Le-L^)  (126) 

If  the  prisTO-enhancement  path  length  is  the  same  in  the  two  measurements, 

Lg-Lg  =  0  and  Ia(((>)  may  be  determined  from  the  experimental  data  as  described 
above • 

A  difficulty  in  the  analysis  of  Eq«L26  is  posed  by  the  fact  that 
AP(4))/P  and  Ap'(<t>)/P  are  both  fluctuating  functions  of  angle.  Their  difference 
also  fluctuates.  It  is  easier  to  analyze  the  data  by  first  obtaining  the 
best  fitting  Lorentzian  curve  to  each  experiment,  and  then  subtracting  the 
two  Lorentzians  to  obtain  a  measure  of  waveguide  performance  over  the  path 
length  Lg-Lg. 

A  weakness  in  the  use  of  Eq.  126  to  analyze  in~plane  scattering  data 
is  the  fact  that  and  Lg  may  not  cancel.  The  derivations  of  Appendix  C 
suggest  that  Le  and  Lg  are  more  likely  to  be  equal  if  the  gap  separating  the 
input  coupling  prism  and  the  waveguide  is  constant  in  the  two  measurements. 
This  is  best  done  by  leaving  the  input  prism  constant  and  moving  the  output 
prism  to  accomplish  the  desired  change  in  Lg.  However,  our  experiments  were 
performed  prior  to  the  derivation  of  prism-enhancement  theory,  and  this  was 
not  done.  For  completeness  we  note  the  theoretical  prediction  of  Appendix  C 
that  Lg  itself  may  be  made  small  by  utilizing  a  broad,  highly  collimated  input 
beam.  This  too  was  not  done  in  our  experiments,  but  should  be  done  in  future 
experiments  of  this  type. 

INTERPRETATION  OF  EXPERIMENTAL  RESULTS  FOR  SAMPLE  154 

In  view  of  the  uncertainties  regarding  the  data  obtained  using 
input  coupling  prisms,  we  will  concentrate  on  the  interpretation  of  Fig.  20 
obtained  using  end-fire  coupling.  We  found 

d  =  11.7  pm 

^  =  2.25  X  (127) 

r, 

from  the  theoretical  fit  to  the  data.  We  would  now  like  to  correlate  these 
results  with  some  of  the  potential  sources  of  scattering  identified  in 
Secs.  II  and  IV. 
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The  quantity  Id  in  Eq.  127  is  the  scattering-center  surface 

density  times  scattering-center  area.  It  corresponds  to  the  fraction  of 

waveguide  surface  covered  by  scattering  centers.  For  contiguous  scattering 

centers,  such  as  surface-roughness  undulations,  Z  =  and  has  its 

largest  possible  value  of  unity.  6x1^  then  has  its  smallest  possible  value 
-5  * 

I. 5  X  10  .  The  size  of  the  scattering  centers,  a  =  11.7  ym,  suggests  that 
we  are  dealing  with  Rayleigh-Gans-Debye  (RGD)  scattering.  Although  Mle 
scattering  is  a  possibility,  we  have  discounted  its  importance,  in  Sec.  IV, 
because  the  scattering  particles  should  be  visible  by  microscope. 

Data  Interpretation  Based  on  Scattering  from  Surface  Compounds 

Restricting  our  attention  to  RGD  scattering  centers  having  a  = 

II. 7  ym  and  6n  ^  1.5  x  10  ^  we  see  from  Table  I  of  Sec,  III  that  the  candi- 
dates  are  surface  roughness,  Ti  voids,  and  LiNb30g  and  Li-Ti-0  surface 
compounds.  The  most  striking  correspondence  between  our  results  and  the 
entries  of  that  table  is  in  the  case  of  scattering  by  surface  compounds.  For 
that  scattering  mechanism  we  combine  Eqs.  76,  77,  and  78  to  find  an  expression 
for  the  mode  index  perturbation  in  terms  of  waveguide  parameters.  It  is 

5n  =  4nAn(n^-l)"^[Vbl/2/(vbl/2+2^)]5^(2„/D)  (128) 

g  ° 

where  n  is  the  LiNbO^  substrate  index,  An  is  the  waveguide  index  perturbation, 

V  =  kQD(2nAn)l/2  is  the  normalized  waveguide  depth  parameter,  b  -  (ng-n)/An 
is  the  normalized  mode  index,  6n  is  the  index  perturbation  associated  with  the 
surface  compounds,  ±s  the  depth  of  the  surface  compounds,  and  D  is  the  wave¬ 
guide  depth. 

It  is  instructive  to  estimate  An  and  D  from  waveguide  fabrication 
parameters  t,  the  Ti  film  thickness,  t^^,  the  diffusion  time,  and  T  the 
diffusion  temperature.  The  equations  to  be  used  are 

An^  =  1.08  t/D 

D  =  2[D(T)t^]l/2  ^29) 

where  the  subscript  e  refers  to  the  extraordinary  index  and  P(T)  is  the 
diffusion  coefficient.  The  numerical  factor  1.08  is  obtained  from  the 
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From  the  same  work 


experimental  results  of  Burns,  Klein,  West,  and  Plew 
we  have  P(T)  =  3.6  x  10"*^^  cm^/sec  at  T  =  950® C,  the  fabrication  temperature 
for  our  waveguide.  The  waveguide  was  also  formed  using  a  film  thickness 

o 

T  =  315A  and  a  diffusion  time  tj^  =  3h.  Appendix  D  presents  some  considera¬ 
tions  on  the  inaccuracy  of  this  analysis  based  on  uncertainty  regarding  the 
correct  value  of  t.  Nevertheless,  we  find 

Ang  =  0.027 

D  =  1.25  pm  .  (130) 

These  results  provide  a  normalized  diffusion  depth  V  =  4.28.  From  the  uni- 

(24) 

versal  curves  for  diffused  waveguides  having  a  Gaussian  index  profile,  we 
find  the  normalized  mode  index  corresponding  to  this  value  of  V  to  be  b  = 

0.44.  Inserting  these  results  into  Eq.  128  we  obtain 

5n^  =  0.046  5n(z^/D)  (131) 

The  results  of  Burns  et  al^^^  shows  that  Zq/D  -  0.1  for  the  case  when  D  - 

2  pm.  If  this  ratio  holds  for  our  smaller  D  -  1.25  pm,  we  find 

5n  =  0,0035  (132) 

required  to  obtain  the  smallest  mode  index  consistent  with  our  experiment, 

=  1.5  X  10~^. 

Some  of  the  surface  compounds  that  may  act  as  the  source  of  (Sn  are 
LiNb30g,  having  the  three  refractive  indices  2.28,  2.36  and  2.40,  and  Li2Ti03, 
having  an  index  of  2,09.  11211307  is  another  possibility  discussed  in  Ref,  7 

probably  having  a  similar  value  for  refractive  index.  It  would  appear  that 
the  largest  value  of  5n  associated  with  these  compounds  is  on  the  order  0.1, 

but  Ref.  7  indicates  that  the  compounds  are  a  reasonably  dilute  mixture  with 

the  host  LiNb03.  The  authors  estimate  Li-Ti-0  compound  concentration  to  be 
less  than  10  mol%.  This  would  probably  reduce  the  effective  value  of  (Sn  to 
about  0.01,  close  to  what  we  have  calculated  in  Eq.  132, 

Another  fact  that  indicates  the  correctness  of  these  ideas  is  that, 
if  iSn  had  the  value  0.1,  as  for  segregated  compounds,  the  compounds  would 
probably  be  visible  with  a  microscope.  In  accord  with  Eq.  131,  6n  =  0.1 
implies  that  6n^  =  4,6  x  10"^.  Substituting  this  result  into  Eq,  127  we  find 
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that  1=8  inin"^  is  the  surface  density  consistent  with  the  data  of  Fig.  20. 

It  is  difficult  to  imagine  that  this  concentration  of  particles,  12  ym  in 
size  and  having  an  index  differential  of  0,1,  would  be  able  to  escape  optical 
detection. 

Data  Interpretation  Based  on  Scattering 
From  Surface  Roughness 


Although  we  have  certinly  indentified  one  possible  source  for 
scattering  of  the  magnitude  observed  in  Fig.  20,  let  us  consider  the  possi¬ 
bility  that  surface-roughness  scattering  is  also  compatible  with  that  result. 

For  surface-roughness  scattering,  the  effective  mode  index  is  found 
using  Eq.  58  and  64, 

5n  =  a(?n  /DD) 
g  g 

3n^/aD  =  k.oAn(2nAn)l/2ot,/3v)  (I33) 

With  <Sng  =  1.5  X  lO"^  ko  =  2it/0.633  pm.  An  =  0.027,  n  =  2.2,  and  3b/3V  =  0.10 
in  the  vicinity  of  V  =  4.3,  as  measured  from  the  universal  curves,  we  find 

o 

=  16  A. 

This  result  shows  that  a  reasonably  smooth  surface  can  still  lead 
to  significant  waveguide  scattering  levels.  This  possibility  was  not  pre¬ 
dicted  in  Sec.  IV  owing  to  our  use  of  An  =  0.007  in  calculations,  rather 
than  An  =  0.027,  as  obtained  for  the  waveguide  in  Fig. 20.  The  scattered 
intensity  varies  as  An^  for  the  surface-roughness  mechanism,  so  the  large- 
An  waveguide  employed  here  was  not  optimum  from  this  point  of  view.  It  is 
therefore  possible  that  the  observed  scattering  was  indeed  caused  by  surface 
roughness.  This  mechanism  could  operate  either  in  addition  to  or  instead  of 
the  surface  compound  formation  mechanism  discussed  above. 

The  possibility  that  surface  roughness  was  the  cause  of  the  scat¬ 
tering  in  Fig.  20  occurred  to  us  well  in  advtince  of  the  realization  that  sur¬ 
face  compound  formation  was  also  consistent  with  that  data.  As  a  result, 
emphasis  was  placed  on  an  effort  to  measure  surface  roughness  in  a  variety 
of  samples  and  correlate  the  results  with  in-plane  scattering  levels  for 
the  same  samples.  This  was  done  to  further  demonstrate  (or  disprove)  the 
correspondence  between  surface  roughness  and  observed  waveguide  scattering. 
The  nature  and  the  results  of  these  experiments  is  described  later. 


Data  Interpretation  Based  on  Scattering  from  T1  Voids 


In  Secs.  II  and  IV  we  considered  the  possibility  of  scattering  from 
regions  where,  for  one  reason  or  another,  Ti  did  not  enter  the  waveguide. 
Particulate  matter  trapped  beneath  the  evaporated  film  was  cited  as  one  source 
of  this  type  of  scattering.  The  mode-index  change  is  approximately  given  by 


6n^  =  -  bAn 

In  our  experiment  6n  =  0.44  x  0.027  =  0.012.  In  order  to  satisfy  Eq.  127, 

^  -*2 
we  require  the  surface  density  of  scattering  centers  to  be  E  =  ,012  nvm  ,  The 

total  number  of  scattering  particles  would  be  2LW  ^  1  since  L  =  15  mm  and 
W  =  1  mm.  It  is  conceivable  that  one  scattering  particle  could  escape 
detection.  Yet  the  observed  scattered  energy  distribution  would  be  in 
reality  a  diffraction  pattern  having  regular  features  associated  with  the 
geometry  of  the  particle.  This  is  in  contrast  to  what  we  observe.  There¬ 
fore  we  discount  this  as  an  important  scattering  mechanism  in  our  experiment. 

Potential  for  Improved  Scattering  Performance 

Later  we  will  discuss  possible  methods  for  reducing  scattering  in 
LiNb03  waveguides  by  reducing  the  strengtli  or  number  of  scattering  centers 
such  as  surface  roughness  and  surface  compounds.  At  this  point  let  us  con¬ 
sider  the  scattering  level  that  might  have  been  achieved  simply  by  fabricating 
the  waveguide  of  Fig. 20  differently,  without  any  attempt  at  reducing  the 
strength  of  scattering  centers.  Had  increased  the  diffusion  temperature 
to  1000“C,  at  which  P(T)  =  9.4  x  10“^^  cm^/sec,^^^  and  the  diffusion  time  to 
8h,  the  diffusion  depth  would  have  been  3.29  pm  instead  of  1.25  pm.  Had  we 

O  o 

decreased  the  Ti  film  thickness  to  150A  from  315A,  the  value  of  An^  would  have 
been  0.0049,  down  from  0.027.  The  normalized  depth  of  the  waveguide  would 
have  been  V  =  4.79,  corresponding  to  a  normalized  mode  index  b  =  0.45.  We 

o 

continue  to  use  db/dV  =  0.1,  6n  =  0.0035  and  o  =  16  A. 


MICROCOPY  RISOLUTION  ILSl  CHART 


In  the  case  of  scattering  from  surface  compounds,  Eq.  128  now  yields 

6n  =*  3.0  X  10”^,  compared  to  1.5  x  for  the  present  waveguide.  This  Implies 

8 

a  factor  25  reduction  in  scattered  intensity,  A  spectrum  analyzer  utilizing 
the  proposed  waveguide  would  have  a  dynamic  range  of  40,8  dB,  up  from  26,8  dB 
predicted  for  the  waveguide  of  Fig  20. 

In  the  case  of  scattering  from  surface  roughness,  Eq,  133  now  yields 

6n  =  1,1  X  10“^,  down  from  1.5  x  10  This  implies  a  factor  184  reduction  in 

8 

scattered  intensity.  A  spectrum  analyzer  utilizing  the  proposed  waveguide 
would  have  a  dynamic  range  of  49.4  dB,  up  from  26,8  dB. 

These  results  are  highly  encouraging  from  the  point  of  view  of 
ultimately  fabricating  waveguides  suitable  for  a  spectrum  analyzer,  A  note 
of  caution  is  that  geodesic  lens  radiation  losses  may  increase  with  the 
introduction  of  waveguides  that  are  optimum  from  a  scattering  point  of  view. 

Also  note  that  prism-enhanced  scattering  increases  with  smaller  An  and 
larger  D.  This  means  that  the  superior  scattering  performance  of  such  wave¬ 
guides  may  not  be  apparent  if  they  are  tested  using  prism  input  coupling. 

The  careful  reader  may  have  noticed  that  the  data  of  Fig.  20  is 
for  TM^  propagation,  while  the  data  is  analyzed  using  a  formalism  developed 

for  TE  modes.  We  justify  this  on  the  basis  that  data  taken  for  the  TE 
o  o 

mode  on  the  same  sample  (154)  was  quite  similar  in  appearance  to  the  data 
of  Fig.  20  (see  Fig.  29  on  page  107). 


OTHER  IN-PLANE  SCATTERING  EXPERIMENTS 


The  in-plane  scattered-energy  distribution  of  Fig,  20  was  obtained 
approximately  midway  through  the  program.  Earlier  data  was  taken  on  a  variety 
of  samples  to  test  the  relative  merits  of  a  substrate  polish  performed  at 
Battelle  and  to  compare  scattering  from  indiffused  and  outdiffused  waveguides. 
Sample  154  of  Fig,  20  was  prepared  to  test  the  effects  of  a  new  sample  cleaning 
procedure,  but  the  evidence  for  surface-roughness  scattering  from  that  waveguide 
had  a  strong  influence  on  later  experiments  performed  during  the  program.  These 
wore  performed  primarily  to  demonstrate  the  correlation  or  lack  of  correlation 
between  measured  levels  of  surface  roughness  and  waveguide  scattering. 

We  now  examine  these  experimental  results  to  see  what  additional  light 
they  shed  on  the  understanding  of  waveguide-scattering  phenomena.  Keep  in  mind 
that  the  data  for  each  waveguide  were  acquired  with  two  different  coupling-prism 
separations  to  cancel  the  effects  of  prism-enhanced  scattering.  The  results  are 
meaningful  only  to  the  extent  that  prism-enhanced  scattering  was  the  same  in  the 
two  measurements. 

Sample  Histories 

Ti-Film  Deposition 

Table  II  identifies  the  samples  used  in  the  program  and  presents 
relevant  details  concerning  their  fabrication.  The  substrates  were  purchased 
from  Crystal  Technology.  They  were  25mm  x  25mm  x  3mm  in  volume  and  polished  to 
A/10  flatness  on  both  25mm  x  25mm  surfaces.  HP-6  was  from  a  high-purity  batch 
while  the  others  were  selected  acoustic  grade.  It  is  likely  that  at  least  four 
different  boiiles  are  represented  by  the  samples  used  in  the  program. 

Ti-film  thicknesses  were  measured  with  a  crystal  monitor  during  their 
deposition  by  e-beam  evaporation.  The  monitor  was  calibrated  using  an  ang¬ 
st  rometer.  No  account  was  taken  of  Che  fact  that  the  measured  thickness  is  not 
that  of  pure  Ti  but  is  that  of  an  unknown  oxide,  depending  on  the  purity  of  the 
vacuum  system.  This  problem  is  addressed  In  Appendix  D.  Ti  films  deposited  on 
samples  154  and  155  probably  had  a  greater  proportion  of  Ti  than  did  other 
samples  documented  in  Table  II,  owing  to  the  use  of  a  tighter  vacuum  system. 


89 


Thermal  Diffusion 


All  diffusions  were  carried  out  in  flowing  oxygen  rather  than  in 

argon,  as  is  done  in  some  laboratories. This  procedure  was  developed  at 

Battelle  to  simplify  fabrication  by  eliminating  the  need  for  a  post-diffusion 

(31) 

heat  treatment  in  oxygen.  However,  the  fact  that  oxygen  plays  an  important 

role  in  the  diffusion  process  means  that  the  diffusion  atmosphere  could  have  an 
impact  on  waveguide  quality.  We  did  not  test  this,  though  it  would  have  been 
worthwhile  to  do  so. 

Following  the  diffusion  heat  treatment,  samples  were  rapidly  quenched 
to  600®  by  moving  them  to  the  edge  of  the  hot  zone  of  the  tubular  oven.  Then 
the  samples  were  brought  to  room  temperature  in  a  time  comparable  to  or  longer 
than  30  min.  The  quench  to  600®C  is  performed  to  minimize  the  time  spent  in  the 
800-850®C  temperature  range,  where  surface  defects  can  be  produced  as  a  result  of 
the  generation  of  the  separated  phase  LiNb^Og . We  have  observed  this  deg¬ 
radation  in  samples  allowed  to  cool  slowly  overnight,  and  in  samples  that  were 
intentionally  allowed  to  remain  at  800-S50®C  for  times  on  the  order  of  ten 
minutes.  A  rapid  quench  from  the  diffusion  temperature  to  600®C  is  probably  not 
required,  although  we  generally  performed  the  quench  in  a  few  seconds. 

Cleaning  Procedure 

The  method  used  to  clean  LiNbO^  substrates  for  waveguide  fabrication 
has  been  subject  to  considerable  evolution  over  the  years  that  Battelle  has  had 
experience  in  the  area.  Two  approaches  are  now  dominant.  One  that  minimizes 
contact  with  the  LiNbO^  surface  has  as  its  principle  features  a  soak  in  "Micro" 
detergent,  follwed  by  a  long  (5  min)  rinse  in  hot,  filtered  tap  water,  and  con¬ 
cluded  with  a  short  rinse  in  hot  distilled  tap  water.  Variations  on  the  main 
theme  include  ultrasonic  agitation  of  the  "Micro"  and  a  concluding  bake  out  at 
250®C.  If  the  sample  is  not  baked,  it  is  blown  dry  with  filtered  air  or  N2. 
Slight  swabing  of  the  surface  with  a  cotton  swab  is  permitted  during  the  early 
rinsing  stages.  Otherwise  there  is  no  contact  with  the  surface. 

A  problem  with  minimizing  surface  contact  has  to  do  with  the  very 
strong  electrostatic  attraction  that  LlNbOg  exerts  on  airborne  particulates. 


Surface  dust  is  not  easily  rinsed  away.  An  approach  that  we  have  found  success¬ 
ful  has  been  to  first  subject  the  sample  to  a  “final  polish”  with  colloidal 
silica.  This  physically  removes  firmly  attached  dust  and  other  residue.  Then 
the  sample  is  wiped  once  with  a  “Kleenex”  or  a  lintless  wipe  that  has  been 
wetted  with  methanol.  If  the  sample  has  been  recently  sawed,  care  must  be 
taken  to  insure  that  loose  chips  are  not  dragged  across  the  surface  in  this 
step.  Nomarskl  microscopy  reveals  a  clean,  nearly  dust-free  surface. 

This  cleaning  procedure  was  used  for  samples  154  and  155,  while  other 
samples  cleaned  during  the  program  were  cleaned  using  the  minimum-contact  ap¬ 
proach.  We  somewhat  prefer  the  contact  appraoch  for  its  simplicity  and  effec¬ 
tiveness,  but  we  also  believe  that  it  doesn’t  make  much  difference  which  method 
is  used  as  far  as  waveguide  scattering  is  concerned.  Our  results  suggest  that 
the  scattering  level  associated  with  cleaning-related  scattering  centers  is  not 
very  great. 

Sample  HP- 6 


HP-6  was  from  a  high-purity  boule  of  LiNbO^  that  was  left  over  from  a 
previous  program.  The  sample  was  ground  and  polished  to  the  extent  that  the 
effects  of  previous  diffusions  and  heat-treatments  were  in  all  likelihood  elimi¬ 
nated.  Then  a  waveguide  was  prepared  on  the  substrate  according  to  the  descrip¬ 
tion  in  Table  II. 

It  was  believed  that  the  scattering  level  from  this  sample,  when 
compared  with  that  from  a  waveguide  formed  on  a  Crystal- Techno logy  prepared 
substrate,  would  indicate  the  potential  for  scattering  from  surface  roughness 
and/or  subsurface  polishing  imperfections  that  might  exist  in  one  or  the  other  to 
a  greater  degree. 

Sample  135 


Waveguide  sample  135  was  fabricated  In  a  similar  fashion  to  HP-6  in 
order  to  make  the  comparison  indicated  above.  Also,  sample  135  contained  a  sur¬ 
face  area  approximately  25  mm  x  6  mm  that  was  not  diffused  with  Ti,  and  hence 
supported  an  outdif fused  waveguide.  The  reason  for  this  was  to  allow  comparison 
of  scattering  from  indiffused  and  outdlffused  waveguides  formed  on  the  same 
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substrate.  Unfortunately,  the  outdif fused  portion  was  only  tested  without  cor¬ 
rection  for  prism-enhanced  scattering,  so  this  comparison  could  not  be  made*  An 
outdif fused  waveguide  was,  however,  carefully  made  and  tested  later  in  the 
program. 

Samples  154  and  155 

These  samples  were  initially  fabricated  to  test  the  influence  of  sur¬ 
face  cleaning  procedures  on  LiNbO^  waveguide  scattering.  However,  following  the 
analysis  of  the  data  for  154,  shown  in  Fig.  20,  it  was  felt  that  surface  rough¬ 
ness  was  a  more  significant  source  of  scattering.  Subsequent  samples  used  in 
the  program  were  fabricated  with  this  possibility  in  mind. 

Samples  166 ,  168 ,  and  169 

These  samples  were  designed  for  a  series  of  measurements  in  which  wave¬ 
guide  scattering  was  to  be  measured  as  a  function  of  surface  roughness.  Rough¬ 
ness  was  to  be  varied  by  post-diffusion  polishing  of  the  waveguide  surface. 

o 

Sample  166  was  fabricated  from  a  relatively  thin,  175-A  Ti  film.  A  thin  film 
was  employed  to  minimize  the  upheaval  of  the  topography  of  the  LiNbO^  surface 
caused  by  the  diffusion  process.  This  upheaval  was  surmised  from  measurements 
made  on  sample  154,  in  which  the  scattered  component  of  light  reflected  from 

o 

the  surface  at  about  90°  suggested  a  roughness  of  about  30  A  (see  Appendix  E) 
or  an  order  of  magnitude  greater  than  the  NBS-  measured  roughness  for  an  undif¬ 
fused  substrate  (see  Appendix  A) . 

o 

Sample  169  was  diffused  from  a  relatively  thick,  600  A  Ti  film.  This 
was  done  to  provide  contrast  with  sample  166,  and  also  to  permit  checking  of  the 
results  of  the  Phase  I  program,  which  showed  the  best  scattering  performance 
to  be  provided  by  a  heavily  diffused  waveguide  that  was  subjected  to  a  long  post¬ 
diffusion  polish. 

Sample  168  was  an  outdiffused  waveguide  formed  to  verify  that  topo¬ 
graphical  upheaval  of  the  surface  was  indeed  related  to  Ti  indiffusion,  and  also 
to  provide  the  comparison  of  scattering  from  indiffused  and  outdiffused  wave¬ 
guides  that  was  mentioned  earlier  in  connection  with  sample  135. 
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Experimental  Results  for  Samples  Other  Than  154 


Scattered-»Energy  Distributions  for  Separated  and 
Adjacent  Input  and  Output  Prism  Couplers 

In  Figs.  21-28  we  present  in-plane  scattered-energy  distributions  for 
the  samples  whose  histories  have  just  been  described.  In  each  experimental  dis¬ 
tribution,  the  solid  curve  shows  scattering  observed  when  the  prism  separation 
was  about  15  mm.  The  dashed  curve  shows  the  scattering  observed  when  the  prism 
separation  was  as  close  as  could  be  conveniently  obtained,  about  1-2  iron.  The 
distributions  have  been  truncated  at  -20  dB  so  the  main  peaks  associated  with 
unscattered  light  are  off-scale.  The  scattering  angles  shown  along  the  horizon¬ 
tal  axis  are  those  which  would  be  measured  in  the  waveguide  material. 

The  smooth  curves  shown  in  the  figures  are  the  Lorentzian-shaped 
theoretical  fits  to  the  data  obtained  by  methods  discussed  earlier  in  this 
section.  These  curves  are  solid  or  dashed  in  correspondence  to  the  experi¬ 
mental  data  that  they  fit.  Table  III  summarizes  the  theoretical  results.  This 
table  presents  the  scattering  level  at  0*^  and  1°  measured  in  air  [(1/n)®  measured 
in  the  waveguide],  as  provided  by  the  theoretical  fit  for  both  separated  and 
adjacent  coupling  prisms. 

Interpretation  of  Results 

It  is  desired  to  evaluate  waveguide  scattering  by  subtracting  the  scat¬ 
tering  level  obtained  using  adjacent  prisms  from  that  obtained  using  separated 
prisms.  The  results  of  this  calculation  are  shown  in  i^ow  5  of  Table  III. 

Three  entries  are  omitted:  Samples  HP-6  mode,  ^  _L  c)  and  135  (TF^  mode, 

k  c)  show  a  negative  scattering  contribution  from  the  waveguide,  while  sample 
155  (TM^  mode,  k  | j  c)  shows  a  very  large  scattering  that  belies  the  close  match 
between  experimental  curves  seen  in  Fig.  2A. 

Results  for  other  samples  show  that  166,  168,  and  169  have  the  lowest 
scattering  levels,  near  -36  dB  and  comparable  to  that  obtained  for  sample  154. 
Samples  135  and  155  have  scattering  levels  approximately  10  dB  greater.  In  view 
of  the  fact  that  the  fluctuations  in  scattered  energy  for  many  of  the  samples 
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-0,5  -0.25  0  0.25  0.50 


Scattering  Angle  <p  (Degrees) 

Fig.  21.  In-plane  scattering  energy  distribution  for  sample  HP-6  (TE  mode, 
kj^c),  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves).  Bell-shaped  curves  are 
theoretical.  Scattering  angles  are  referenced  to  waveguide  material. 
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Fig.  24. 
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Scattering  Angle  (Degrees) 

In-plane  scattering  energy  distribution  for  sample  155  (TM^  mode, 
k((c),  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves) .  Bell-shaped  curves  are 
theoretical.  Scattering  angles  are  referenced  to  waveguide  material 
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Scattering  Angle  (Degrees) 

25.  In-plane  scattering  energy  distribution  for  sample  155  (TM^  mode, 
kj_c)  ,  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves) .  Bell-shaped  curves  are 
theoretical.  Scattering  angles  are  referenced  to  waveguide  material. 
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Fig.  26.  In-plane  scattering  energy  distribution  for  sample  166  (TM^, 
kj[c)  ,  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves).  Bell-shaped  curves  are 
theoretical.  Scattering  angles  are  referenced  to  waveguide  material 
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Tn-plane  scattering  energy  distribution  for  sample  168  (TEq, 
k[c),  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves).  Bell-shaped  curves  are 
theoretical.  Sc'attering  angles  are  referenced  to  waveguide  material 
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Fig.  28.  In-plane  scattering  energy  distribution  for  sample  169  (TMq, 
kj^c)  ,  taken  with  separated  coupling  prisms  (solid  curves)  and 
adjacent  coupling  prisms  (dashed  curves) .  Bell-shaped  curves  are 
theoretical.  Scattering  angles  are  referenced  to  waveguide  material 
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(row  6) was  comparable  to  the  changes  in  scattering  produced  using  separated  and 
adjacent  coupling  prisms,  it  is  difficult  to  know  what  significance  to  attach 
to  these  results.  Let  us,  therefore,  consider  the  entries  of  Table  III  from  a 
different  point  of  view. 

Taken  as  a  whole,  the  data  show  the  dominance  of  scattering  from  the 
prism-coupling  region  over  that  from  the  free  waveguide.  Since  prism  scattering 
results  from  an  enhancement  of  scattering  from  waveguide  imperfections  in  the 
input-coupling  region,  those  waveguides  with  the  lowest  inherent  scattering  will 
also  have  the  lowest  prism-enhanced  scattering.  We  can  thus  determine  the  best 
samples  in  Table  III  by  looking  for  those  with  the  lowest  scattering  as  measured 
with  adjacent  coupling  prisms.  A  critical  assumption  is  that  the  proportionality 
constant  relating  inherent  waveguide  scattering  to  prism  scattering  is  the  same 
for  all  vv'aveguides.  Exceptions  to  this  will  have  to  be  discussed  individually. 

An  examination  of  row  3  in  Table  III,  for  on-axis  scattering  and 
adjacent  coupling  prisms,  shows  HP- 6  to  be  the  worst  waveguide  tested,  having  a 
scattering  level  of  -22  dB.  Recall  that  this  was  the  only  sample  not  polished 
by  the  supplier.  Crystal  Technology,  Waveguides  153  and  166  are  now  found  to  be 
comparable  in  performance  and  superior  to  135.  In  row  5,  166  was  superior  to 
both  155  and  135,  which  were  comparable.  We  are  thus  inclined  to  rank  the  three 
samples  in  order  of  increasing  qualitv  as  follows:  135  (poorest),  155  (inter¬ 
mediate),  and  166  (best).  This  is  significant  because  the  waveguides  were  fabri¬ 
cated  with  slightly  different  procedures:  135  was  diffused  from  a  thicker  Ti 

o  o 

film  than  wt're  155  and  166  (270  A  instead  of  175  A),  while  155  was  formed  from  a 
film  evaporated  in  a  tighter  vacuum  system  than  was  the  film  used  for  166.  Our 
results  suggest  that  the  use  of  a  thin,  oxygen-rich  Ti  film  will  improve  wave¬ 
guide  quality,  presumably  by  impeding  the  formation  of  Li-Ti-0  surface  compounds 
and  by  minimizing  any  surface  roughness  that  may  be  induced  by  the  diffusion 
process . 

Oiitdif fused  waveguide  sample  168  presents  a  paradox  in  that  it  shows 
quite  low  scattering  obtained  as  the  difference  between  data  for  separated  and 
adjacent  coupling  prisms,  yet  the  absolute  scattering  level  measured  with  adja¬ 
cent  coupling  prisms  in  higher  than  for  any  sample  in  Table  III  except  HP-6, 
Outdif fused  waveguides  are  expected  to  have  less  scattering  than  indif fused 
waveguides  because  they  are  immune  to  scattering  from  Ll-Ti-0  surface  compound 
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formation  and  from  diffusion-induced  surface  roughness.  This  can  be  reconciled 

with  the  high  level  of  scattering  observed  for  adjacent  coupling  prisms  only  if 

we  assume  a  large  degree  of  prism-enhanced  scattering.  This  is  in  agreement 

with  the  predictions  of  Appendix  C.  Equation  C-17  indicates  that  prism-enhance- 

2 

ment  varies  in  proportion  to  (D/An)  ,  where  D  is  the  waveguide  depth  and  An  is 
the  surface  index  change.  This  ratio  can  be  two  orders  of  magnitude  larger  in 
an  outdif fused  waveguide  than  in  an  indlf fused  waveguide.  We  are  thus  confi¬ 
dent  in  asserting  that  sample  168  is  one  of  the  better -quality  waveguides  tested. 

The  last  sample  to  consider  in  Table  III  is  sample  169.  This  wave¬ 
guide  shows  the  best  scattering  performance  of  any  waveguide  in  the  table,  by 
either  of  the  two  criteria  that  we  have  employed.  This  is  surprising,  because 
the  fabrication  conditions  were  deliberately  chosen  to  produce  an  inferior  wave- 

o 

guide.  A  600  A-thick  Ti  film  was  diffused  to  produce  Li-Ti-0  surface  compound 
formation  and  diffusion  induced  surface  roughness.  It  was  planned  to  monitor 
waveguide  quality  improvement  by  performing  post  diffusion  polishing  of  the 

o 

sample,  as  was  done  in  Ref.  1.  In  that  work,  720  A-thick  Ti  films  diffused  for 
three  hours  showed  considerable  surface  granularity  which  we  associated  with 
surface  compound  formation.  In  a  waveguide  formed  similarly  to  169,  only  very 
slight  granularity  was  observed.  Our  interpretation  is  that  the  higher  diffu¬ 
sion  temperature  (1000°C  instead  of  950®C)  and  longer  diffusion  time  (8  h  instead 
('f  3  h)  cause  the  surface  compounds  to  break  up,  forming  a  more  homogeneous  wave¬ 
guide  layer.  However,  caution  must  be  exercised  in  interpreting  the  results  of 
Table  III  for  this  sample,  since  the  extent  of  prism-enhanced  scattering  is  not 

known  and  can  not  be  easily  predicted.  Both  An  and  D  are  larger  for  this  wave- 

2 

guide,  and  the  ratio  (D/An)  describing  prism  enhancement  may  be  greater  than, 
less  than  or  comparable  to  that  for  other  indif fused  samples  in  the  table. 

Connection  with  Previous  Results 

In  summary,  it  appears  that  the  best  waveguides  are  obtained  by  fabri¬ 
cation  procedures  that  minimize  disruption  of  the  sample  surface,  either  by  sur¬ 
face  compound  formation  or  diffusion  induced  roughness.  We  find  from  experiment 
that  this  can  be  done  by  the  use  of  thin  Ti  films  having  a  high  oxygen  content 
and/or  by  the  use  of  long  diffusion  times  and  high  temperatures.  Both  tech- 
r?  iques  should  lead  to  smaller  values  of  An  and  larger  values  of  D,  which  the 
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theory  of  Sec*  IV  shows  to  be  conducive  to  low  scattering  cross  sections.  Thus 
our  experimental  results  are  qualitatively  supportive  of  the  theory.  Unfortu¬ 
nately,  they  do  not  go  far  enough  toward  determining  the  best  waveguide  scat¬ 
tering  that  may  be  obtained.  That  is  the  case  because  we  concentrated  in  our 
experiments  on  sources  of  scattering  that  were  related  to  waveguide  handling 
procedures  such  as  polishing  and  cleaning,  rather  than  to  waveguide  diffusion 
conditions,  such  as  diffusion  time,  temperature  and  Ti-film  thickness.  This 
choice  was  made  because  the  results  of  the  Phase  I  program^^^  suggested  that  it 
was  appropriate.  We  now  believe  that  many  of  these  results  were  improperly 
influenced  by  prism-enhanced  scattering. 

The  last  samples  examined  in  the  current  program,  166,  168,  and  169, 
show  an  appropriate  shift  in  emphasis  to  the  study  of  diffusion  conditions  and 
their  effects  on  scattering.  Unfortunately,  significant  time  was  lost  in 
attempting  to  develop  end-fire  coupling  techniques  and  an  in-depth  study  was 
not  possible.  The  most  serious  loss  was  a  set  of  experiments  to  evaluate  the 
effects  of  post-diffusion  polishing  of  various  waveguide  samples.  These  experi¬ 
ments  would  have  provided  information  regarding  the  relative  importance  of  dif¬ 
fusion-induced  surface  roughness  and  surface  compound  formation,  and  could  have 
resulted  in  a  more  definitive  recipe  for  low  scattering  waveguides.  As  the 
current  program  ends,  we  know  the  desirability  of  fabricating  waveguides  with 
small  An  and  large  D,  but  we  do  not  know  the  optimum  way  to  achieve  these  con¬ 
ditions,  whether  by  heavy  post-dif fusion  polishing  of  waveguides  formed  using 
thick  Ti  films  or  by  light  post-diffusion  polishing  of  waveguides  formed  using 
thin  Ti  films  and  longer  hotter  diffusion  treatments. 

Our  initial  attempts  to  answer  these  questions  involved  the  development 
of  a  reflection  scattering  technique  to  measure  surface  roughness.  This  tech¬ 
nique  was  mentioned  previously  in  connection  with  sample  154,  where  it  was  used 
i:o  verify  that  the  surface  roughness  of  that  sample  was  enough  to  cause  the 
observed  level  of  scattering.  We  planned  to  make  repeated  use  of  this  tech¬ 
nique  following  post-diffusion  polishes  of  waveguides  in  the  160  series  to  more 
precisely  correlate  waveguide  scattering  with  surface  roughness.  Some  initial 
results  with  these  waveguides  are  described  in  Appendix  E, 

In  order  to  complete  our  presentation  of  In-plane  scattering  data  we 
include  as  Fig.  29  the  scattered  energy  distribution  for  sample  154  taken  using 
end-fire  input  coupling  of  TE^  polarized  light. 
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VI.  SUMMARY  AND  CONCLUSIONS 


This  report  has  described  an  intensive  investigation  of  sources  of 
scattering  in  LlNbO^,  Tl-dif fused  optical  waveguides.  The  investigation,  pri¬ 
marily  experimental  in  nature,  has  been  complicated  by  experimental  difficul¬ 
ties  of  the  type  which  escape  Immediate  detection,  and  which,  when  discovered, 
frustrate  correction.  I  speak  primarily  of  the  problems  associated  with  prism- 
enhanced  scattering  and  the  unsuccessful  attempt  to  eliminate  them  by  the  polish¬ 
ing  of  edges  for  end-fire  coupling.  These  problems  have  limited  the  scope  of 
the  program  to  the  extent  that  waveguides  other  than  Ti-diffused  LiNbO^  were 
not  investigated .  We  have,  however ,  managed  to  interpret  the  data  obtained 
for  LiNbO^  through  the  use  of  extensive  (and  initially  unanticipated)  calcula¬ 
tions,  so  that  important  conclusions  are  drawn  with  regard  ro  the  suitability 
of  the  material  for  integrated-optical  signal-processing  applications  such  as 
the  spectrum  analyzer.  While  we  have  not  experimentally  studied  other  candi¬ 
date  waveguides  for  these  applications,  we  note  that  our  calculations  are 
relevant  to  all  waveguides  having  a  small  surface-to -substrate  index  change  An. 

In  addition,  our  treatment  of  the  scattering  associated  with  prism  coupling  and 
our  findings  regarding  the  polishing  of  edges  for  end-fire  coupling  should  prove 
useful  to  those  whose  concern  is  with  silicon-substrate  waveguides  as  well  as 
with  LiNbO^. 

A  listing  of  the  achievements  of  the  program  follows: 

•  Identification  of  potential  sources  of  scattering  in  Ti- 
dif fused  LiNbO^  waveguides. 

•  Comprehensive  listing  of  experimental  methods  useful  for 
the  study  of  waveguide  scattering  sources,  including 
experimental  results  for  those  we  employed. 

•  Calculation  of  spectrum-analyzer  dynamic  range  in 
terms  of  scattering  cross  sections,  scattering-center 
densities,  and  spectrum-analyzer  geometry. 

•  Derivation  of  formulas  for  in-plane  scattering  cross 
sections  for  Mie,  Rayleigh-Gans-Dehye ,  and  Rayleigh 
scattering  centers 
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•  Analysis  of  scattering  associated  with  prism  coupling. 

•  Investigation  of  edge-polishing  techniques  for  end-fire 
coupling. 

•  Development  of  procedures  for  data  acquisition  and  analysis 
in  the  presence  of  prism-enhanced  scattering. 

•  Development  of  a  qualitative  but  highly  sensitive  experi¬ 
ment  for  detection  of  in  plane  scattering  at  90®. 

•  Measurement  of  surface  roughness  associated  with  Ti 
diffusion  using  a  reflection-scattering  experiment. 

•  Testing  and  evaluation  of  in-plane  scattering  in  LiNbO^ 
waveguides  formed  using  a  wide  variety  of  fabrication 
procedures . 

The  research  summarized  by  these  achievements  has  led  us  to  conclude 

that  LiNbO^  is  an  eminently  useful  substrate  for  the  spectrum-analyzer  appli- 
(2)-^ 

cation. 

We  find  that  the  initial  substrate,  as  polished  by  the  supplier,  has 

one  of  the  smoothest  surfaces  known  to  !JBS.  The  process  of  Ti  diffusion  roughens 

the  surface  slightly  and  produces  the  formation  of  non-LiNbO^  compounds  within 

(7)  ^ 

the  top  fraction  of  a  micron  of  the  waveguide  surface  .  Most  likely  these 
are  compounds  of  lithium,  titanium,  and  oxygen,  such  as  Li2Ti0^  or  Li2Ti^02^^^ 

Vie  find  that  the  index  inhomogeneity  produced  by  these  compounds  or  by  the 
surface  roughness  associated  with  them  are  each  sufficient  to  explain  the  scat¬ 
tering  levels  observed  in  our  best  experimental  waveguides.  These  waveguides, 
however,  were  fabricated  in  such  a  way  as  to  produce  a  relatively  large  value 
of  the  surface  index  change  An  and  a  small  value  of  the  diffusion  depth  D.  We 
now  know  that  these  conditions  are  decidedly  non-optimum  from  the  point  of  view 
of  reducing  waveguide  scattering,  a  fact  that  was  not  obvious  experimentally 
because  they  turn  out  to  be  optimum  conditions  for  minimizing  the  significant 
contribution  from  prism-enhanced  scattering.  On  the  basis  of  theory  developed 
during  the  program  we  predict  that  employing  waveguide-diffusion  treatments  to 
produce  An  -  0.005  and  D  -  3.3  um,  to  give  a  specific  but  not  all-inclusive  ex¬ 
ample,  will  reduce  the  scattering  contributed  by  surface  roughness  to  negligible 
levels . 
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Residual  scattering  would  then  be  associated  with  surf ace-compound 
formation.  This  contribution,  while  remaining  greater  than  that  from  surface 
roughness,  is  also  reduced  by  fabricating  waveguides  having  small  An  and  large 
D.  This  reduction  occurs  simply  because  the  modal  distribution  of  guided  light 
produces  a  smaller  field  at  the  surface  under  conditions  of  small  An  and  large  D, 
However,  there  is  some  evidence  that  the  longer,  hotter  diffusion  treatment 
leading  to  these  conditions  will  also  increase  the  material  homogeneity  of  the 
surface. This  produces  a  further  reduction  in  scattering  quite  apart  from 
that  associated  with  the  modal-field  distribution.  On  this  basis  we  conclude 
that  waveguides  suitable  for  the  spectrum  analyzer  should  result  from  simple 
adjustments  of  the  diffusion  conditions.  However,  we  emphasize  that  the  wave¬ 
guides,  when  evaluated  using  prism  coupling,  may  not  appear  to  be  of  high 
quality  because  of  the  increase  in  prism-enhanced  scattering  when  An  is  reduced 
and  D  is  made  larger. 

Finally,  we  note  that  our  expectation  that  low  waveguide  scattering 
results  from  longer,  hotter  diffusion  treatments  does  not  guarantee  the  success 
of  a  spectrum  analyzer  fabricated  using  such  a  treatment.  Reference  32  shows 
that  geodesic  lenses  formed  in  waveguides  having  a  small  An/D  ratio  may  have 
through-puts  intolerably  low  for  applications.  If  this  happens  to  be  true,  it 
may  be  necessary  to  employ  waveguides  having  An  and  D  values  that  are  non-optimum 
from  the  point  of  view  of  in-plane  scattering  characteristics.  Scattering  from 
these  waveguides  may  still  be  reduced  to  tolerable  levels  by  post-diffusion 
polishing  of  the  surface.  This  was  investigated  in  the  first  phase  of  this  two 
year  scattering  program. Figure  12  of  Ref.  1,  reproduced  as  Fig.  C2 
of  Appendix  C,  shows  the  best  scattering  performance  observed  during  either  phase 
of  the  scattering-reduction  program.  It  was  obtained  by  post-diffusion  polishing 
of  a  sample  that  initially  showed  evidence  of  significant  surface-compound  forma¬ 
tion. 

In  summary,  we  believe  that  Ti  diffused  waveguides  formed  in  LiNbO^ 
using  suitable  diffusion  conditions  and  post-diffusion  polishing  as  required 
will  suffice  for  the  production  of  integrated  optical  spectrum  analyzers  having 
dynamic  range  values  in  excess  of  40  dB,  and  perhaps  as  large  as  60  dB. 
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APPENDIX  A 


NATIONAL  BUREAU  OF  STANDARDS 
EVALUATION  OF  SURFACE  ROUGHNESS  OF  LiNbO^ 

Attached  is  a  report  of  test  results  obtained  by  personnel  of  the 
National  Bureau  of  Standards.  The  object  of  the  test  was  to  obtain  roughness 
and  autocorrelation-length  information  for  LiNbO^  substrate  164,  This  sample 
was  sent  to  NBS  in  as-  received  condition .  As  for  other  samples  used  in  the 
program,  the  supplier  was  Crystal  Technology,  The  sample  dimensions  were 

mm  X  25  mm  X  3  mm  and  the  polish  of  both  broad  faces  was  specified  to  be  ^ 

' /4  in  flatness.  The  sample  was  selected-acoustic  grade  in  quality. 

In  summ  irv  of  the  attached  report,  the  workers  described  the  sample 
surface  as  i>ne  oi'  the  smoothest  they  had  ever  examined.  The  dates  referred  to 
in  the  report  are  tor  the  year  1980 . 
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To:  Battelle  Memorial  Institute 
Columbus  Laboratories 
505  King  Avenue 
Columbus ,  Ohio  43201 
Attn;  David  Vahey 


This  report  covers  NBS  measurements  of  the  surface  texture  parameters  of 
a  polished  LiNb03  specimen.  The  specimen  was  measured  on  two  occasions: 

Jan.  25  and  Apr.  7  of  this  year.  Between  these  dates  the  data  analysis 
program  was  changed  so  that  the  precision  of  the  Rq  calculation  was  improved 
on  April  7.  The  results  for  the  rms  roughness  and  autocorrelation 
length  a  are  tabulated  below. 


Table  I 

a 

(nn) 

(ym) 

Jan.  25 

<0.3 

36 

% 

Apr.  7 

<0.35 

23 

The  above  values  have  estimated  uncertainties  of  approximately  30%  due 
primari ly  to  uncertainty  in  the  calibrating  step  height,  uncertainty  in 
the  measurement  of  horizontal  displacement,  digitization,  and  the  variation 
in  the  surface  proper t ies  themselves  f ron  place  to  place.  The  30% 
uncertainty  represents  a  95%  confidence  interval. 

Measurement  Procedure 

Surface  topography  is  measured  at  the  NBS  by  means  of  a  minicomputer /stylus 
instrument  system.  Using  an  interferometrically  measured  step,  the  system 
was  calibrated  on  each  value  of  magnification  employed  during  a  measurement. 
Profiles  of  the  calibrating  step  and  the  roughness  area  under  test  were 
stored  in  the  minicomputer  memory  using  12  bit  analog  to  digital  conversion. 
Each  profile  contained  4000  digitized  points. 

The  stylus  tip  was  chisel-shaped  with  approximate  dimensions  0.1  pm  x  1  pm. 

The  response  of  the  system  was  limited  to  a  band  of  surface  v^avelengths 
determined  by  the  stylus  width  and  various  electronic  filters  described  below. 
Rg  values  were  calculated  from  three  successive  traverses  of  each  position. 

The  average  of  these  was  calculated  and  is  shown  in  the  data  sheets.  Then 
the  digitized  profile  data  from  the  third  traverse  at  each  position  were 
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stored  in  the  computer.  The  other  parameters  and  the  statistical  functions 
were  subsequently  calculated  f rom  the  stored  profiles. 

The  experimental  parameters  are  shown  in  Table  II. 

Table  II 


Date 

Trace  Length 
(pm) 

Approximate 
Point  Spacing 
(pm) 

Positions 

Bandwidth 

(pm) 

Jan . 

2S 

750 

0.1875 

1-6 

4  -  200 
(low  pass  - 
high  pass  filters) 

Apr , 

7 

80 

0.02 

1-6 

1-19 

(stylus  width  - 
high  pass  filters) 

Apr. 

7 

1500 

0.375 

7-12 

A  -  200 
(low  pass  - 
high  pass  filters) 

Please  note  that  there  were  two  modes  of  operation  on  Apr.  7,  The  results 
quoted  in  Table  T  were  taken  for  the  lonj’-trace  mode  only.  The  results  for 
the  sliort  trace  mode  are  consistent  with  these. 

The  value  for  was  calculated  according  to  the  formula 


where  yj  is  the  hei<^ht  of  the  filtered  profile  with  respect  to  the  mean 
line  at  position  1, 

The  autocorrelation  length  was  calculated  from  the  average  autocorrelation 
function  each  day.  It  is  defined  as  the  lag  distance  at  which  the 
function  drops  to  10%  of  its  value  at  zero  lag.  Since  the  noise  of  the 
instrument  was  appreciable  at  the  high  magnification  (1  million  X)  required 
to  make  the  measurement,  the  spike  at  ^.ero  lag  was  ignored  and  the  ACF  was 
replotted  on  semilog  paper  as  several-point  averages.  The  values  for  t;  were 
taken  from  these  plots. 

Add i t i ona  1  Da  ta 

A  number  of  other  surface  parameters  and  functions  were  ca leu la  ted  from  the 
surface  profiles  and  the  results  are  i  nc  1  tided  with  this  report .  These 
inc 1 ude  the  parameters  of  average  slope ,  average  wavelength ,  peak-count 
wavelength,  peak-to-mean  line  height,  skewness,  and  kurtosis,  and  the  power 
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spectral  density  function  (PSD),  Results  were  calculated  for  each  profile 
and  averages  of  these  were  calculated  for  each  surface. 

The  PSD  represents  the  breakdovn  of  the  profile  into  its  component  spatial 
frequencies  f.  For  most  surfaces  the  PSD  monotonically  decreases  with 
increasing  spatial  frequency.  The  cutoffs  shown  at  low  frequency  are  due 
to  the  low-pass  electronic  filters.  The  sharp  spikes  are  probably  not  real 
spatial  components  of  the  surface  but  rather  components  of  periodic  noise 
picked  up  during  the  measurement. 

Measurements  made  by: 


Supervisor 


For  the  Director, 

D,  R,  Flynn,  Yrogram  Manager 
Mechanical  Production  Metrology 
Center  for  Mechanical  Engineering 
and  Process  Technology 
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APPENDIX  B 


90°  SCATTERING  EXPERIMENTAL  ANALYSIS:  DETERMINATION 
OF  THE  THRESHOLD  OF  VISION 

In  Sec .  IV  wo  described  an  experiment  to  test  for  scattered  light  at 
90°.  Owing  to  the  low  level  of  scattering  encountered  in  LiNbO^  waveguide 
s^iTTiple  135  at  that  scattering  angle,  we  found  it  necessary  to  use  the  human  eye 
as  the  detector.  Since  the  eye  has  a  nonlinear  response  to  light  we  limited 
consideration  to  deciding  whether  or  not  any  light  could  be  detected,  rather 
than  how  much  light  could  be  detected.  To  convert  our  qualitative  observation 
to  quantitative  information,  it  is  necessary  to  specify  the  threshold  of  human 
vision . 

Reference  B1  indicates  that  for  a  point  source  of  red  light,  the  thres- 

2 

hold  of  vision  is  about  0.15  vilni/m  for  foveal  vision  and  about  0.004  ulm/m*" 

for  extra-fovea]  vision,  Foveal  vision  is  described  as  that  employed  under 

ordinary  lighting  conditions,  when  we  "look  right  at"  an  object.  For  the  dark- 

adapted  eye,  there  is  a  gain  in  sensitivity  of  about  40X  obtained  by  looking  at 

a  slight  angle  away  from  the  object.  Assuming  that  this  condition  applied  at 

2 

least  during  part  of  our  experiment,  the  threshold  of  vision  is  (0.004  ulm/m  )  x 

2 

(lW/682  Im)  X  the  aperture  of  the  eye.  Taking  the  latter  quantity  to  be  5  mm  , 

we  find  P^^  =  3  x  10  W.  Even  if  foveal  vision  is  employed  throughout  the 

experiment,  which  seems  unlikely,  the  threshold  is  still  comparable  to  the 
-15 

10  W  value  employed  in  the  analysis  of  Sec.  IV. 

If  the  scattering  source  is  viewed  as  an  extended  source  positioned, 

as  far  as  the  eye  can  tell,  at  the  output  coupling  spot,  it  is  probably  more 

appropriate  to  calculate  the  threshold  of  vision  using  the  minimum  perceptible 

—6  2 

luminance,  indicated  to  be  5  x  10  cd/m  in  Ref.  Bl.  The  power  associated  with 

this  luminance  is  obtained  as  its  product  with  the  solid  angle  subtended  by  the 

eye  in  the  observation  and  the  area  of  the  illuminating  region.  If  the  aperture 
2 

of  the  eye  is  5  mm  and  the  viewing  distance  is  300  mm,  the  solid  angle  subtended 
2^2  A 

is  5  mm  /3  x  10  mm  =  1.7  x  10  ,  The  area  of  the  illuminating  region  is  the 

2 

area  of  the  output  coupling  spot,  about  10  mm  x  0 . 5  mm  =  5  mm  .  Since  1  cd  = 

1  Im/steradian  =(1/682)  W/steradian,  we  find  the  power  associated  with  the 
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threshold  of  vision  to  be 


P  =  (5  X  10  cd/m  )(1.7  x  10  ^  sterad)(5  x  10  ^  m^)  (1/682)  (W/steradian) 

—  1  ft 

=  6  X  10  W  (Bl) 

This  is  slightly  lower  than  the  value  calculated  for  a  point  source, 

A  potential  problem  is  that  our  experiment  could  not  be  done  in  a  com¬ 
pletely  darkened  room,  owing  to  the  glow  of  the  laser  discharge  and  scattering 
from  beam  forming  optics  and  light  in  the  waveguide  substrate.  Nevertheless  the 
value  =  10  W  used  in  Sec.  IV  is  two  orders  of  magnitude  larger  than  that 
which  we  have  calculated  here,  and  this  could  compensate  somewhat  for  the  error 
introduced  by  utilizing  visual  sensitivity  data  presumed  to  apply  in  a  totally 
darkened  environment.  We  note  that  1  nW  of  HeNe  laser  light  incident  on  a  white 
card  is  visible  to  the  eye  under  reduced  but  not  totally  darkened  lighting  con¬ 
ditions.  Assuming  that  100%  of  this  power  is  reflected,  an  upper  limit  of 
-14 

10  W  is  incident  on  the  eye  at  a  distance  of  30  cm.  This  is  only  one  order 
of  magnitude  larger  than  the  value  we  have  employed  for  the  threshold  of  vision. 

Re^erenc^' 

Bl.  J.V>/.T.  Walsh,  Photometry,  Constable  and  Company  Ltd.,  London,  1953, 

('hapter  III. 
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SCATTERING  ASSOCIATED  WITH  PRISM  COUPLING 

In  this  appendix  we  present  a  theoretical  analysis  of  what  is 

referred  to  in  the  text  as  **prism-enhanced  scattering."  Our  major  conclusions, 

cited  often  in  the  text,  are  the  following:  (1)  An  ideal  prism  can  effectively 

magnify  the  scattering  contribution  from  scattering  centers  in  the  waveguide 

beneath  the  input-coupling  prism,  relative  to  the  scattering  contribution 

from  the  free  waveguide.  (2)  The  enhancement  is  greater  when  the  angular 

divergence  of  the  input  beam  acts  to  produce  a  low  coupling  efficiency,  and 

should  be  reduced  to  tolerable  levels  by  using  a  well-collinated  input  beam. 

Figure  Cl  shows  a  diagram  of  the  prism  waveguide  geometry  in  the 

vicinity  of  the  input-coupling  region.  In  order  to  be  specific,  we  consider 

scattering  centers  associated  with  Li-Ti-0  compound  formation  just  beneath 

the  waveguide  surface.  We  will  take  that  view  of  scattering  in  which  the 

source  of  the  scattered  field  is  the  polarization  vector  6P  -  P  -  P .  ,  i 

~  —actual  —ideal, 

where  polarization  in  the  presence  of  scattering  centers  and 

-ideal  polarization  that  would  result  if  there  were  no  scattering  centers. 

From  the  definition 

P  =  (e  -  e^)E  ,  (Cl) 


we  determine 


6P  =  6eE  (C2) 

wheri  F  is  the  total  field  at  the  surface  and  is  the  change  in  per- 

2 

mittivitv  associated  with  the  scattering  mechanism.  Since  =  n  r  in 

o 

the  optical  regime,  2n6ne^  may  be  substituted  in  Eq,  2. 

The  scattered  field  from  a  given  scattering  center  will  be 

proportional  to  the  dipole  moment,  which  is  the  integral  of  the  polarization 

6P  over  the  volume  of  the  center.  For  the  case  of  islands  of  Li-Ti-0  compounds 

2 

formed  near  the  waveguide  surface,  we  use  a  typical  area  a  ,  measured  in  the 
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plane  of  the  waveguide,  and  a  depth  z  Thus 


2 

E  «  6nEa  z  , 
s  o 


(C3) 


where  E  is  the  (scalar)  field  from  one  scattering  center.  The  intensity 

®  2 
from  N  scattering  centers  is  proportional  to  NE^.  For  this  discussion  we 

are  interested  only  in  the  fact  that  N  is  proportional  to  the  beam  path 

length  L.  For  the  scattered  intensity,  we  write 


2  2  4  2 

I  =  K6n  E  a  z  L 
o 


(C4) 


where  K  is  a  constant  of  proportionality. 

Our  interest  is  in  comparing  the  scattered  intensity  from  the 
prism-coupling  region  to  that  from  the  free  waveguide.  We  assume  that  all 
parameters  in  Eq.  C4  except  for  E  and  L  are  unaffected  by  the  presence  of 
the  prism.  Thus  we  write 

I  =  K'e\  ,  K'  =  KSn^a'^z  (C5) 

g  g  g  o’ 

for  the  scattered  intensity  from  the  free  waveguide,  where  is  the  wave¬ 
guide  field  at  the  surface,  and  is  the  path  length  between  input  and 
output  coupling  prisms.  Similarly,  the  scattered  intensity  from  the  input- 
prism  coupling  region  is 

I  =  K’e\  ,  (C6) 

c  c  c 

where  E^  is  the  total  field  at  the  waveguide  surface  and  is  the  coupling 
length. 

The  field  E^  contains  contributions  from  the  modal  waveguide 

field  E  and  from  the  evanescent  field  associated  with  light  that  is  totally 
6 

reflected  at  the  base  of  the  input-coupling  prism.  If  this  field  is  called 
Ep,  then  the  largest  value  of  E^  at  the  waveguide  surface  is 
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E  +  E 
P  g 


(C7) 


If  Ep>>Eg,  the  light  intensity  scattered  in  the  coupling  region  relative  to 
that  scattered  in  the  free  waveguide  is 


2  2 

I  /I  =  E  L  /EL 
eg  P  c  g  g 


(C8) 


This  ratio  can  be  unity  or  greater  even  if  L^<<L  ,  since  we  assume  E  >>E  , 

^  ^  g  P  g 

To  see  if  this  inequality  can  hold  in  practice,  a  calculation  based  on  the 

(Cl) 

prism-coupling  theory  of  Tien  and  Ulrich  is  in  order. 

Equations  11  and  12  of  the  latter  reference  provide  an  expression 
for  the  field  E^  in  terms  of  the  amplitude  of  a  plane  wave  incident  on 
the  base  of  the  coupling  prism.  If  A,,  is  taken  to  be  the  amplitude  of  the 


plane  wave,  consistent  with  the  notation  of  Tien  and  Ulrich, 


(Cl) 


we  have 


^321 


(C9) 


where  ^-^21  amplitude  transmission  across  the  gap  of  the  three  layer 

system  consisting  of  prism,  air  gap,  and  waveguide  layer.  Equation  12  of 
(C2) 

Ulrich  contains  a  general  expression  for  ^221* 

,,  ,  2  2.1/2.  2  ,.-1/2  -2i(t).-l 

4h(n  -n  )  (n  -1)  e  ^(1+h  e  ) 

321  p  p 


(CIO) 


^  -1,.2  ,,1/2,  2  2, -1/2, 

4;  =  tan  -1)  ('Ip"''  )  1  « 


where  n^  is  the  prism  index,  n  is  the  LiNbO^  index,  S  is  the  gap  width,  and 
h<<l  describes  the  weak-coupling  limit. 

To  obtain  the  corresponding  expression  for  the  field  from  the 
waveguide  mode  at  the  surface,  we  use  the  differential  equation  derived 
by  Tien  and  Ulrich  to  describe  the  interaction  between  the  waveguide  field 
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and  the  prism  field.  This  is  Eq,  (42)  of  Reference  Cl.  In  terms  of  real 
variables,  rather  than  normalized  variables,  and  in  terms  of  our  notation, 
the  differential  equation  is 


dA^/dx  =  aA^  -  yA^ 

2  2  1/2 

a  =  (i _  /2nD)(n^-n  ) 

321  1  g 

(Cll) 

/I  /  2  2.  .  2  2.1/2.  2  1.-1.2.,  . 

Y  ="  (1/nD)  (n^-n^)  (Up-n  )  ^ 

2  7  7  2  7  7  2  -1 

f  (h)  =  4h  [(1+h^)  -4h  (n  -l)(np-l)  ]  . 

In  this  equation  n^  =  n  +  An  is  the  waveguide  surface  index, 
n^  =  n  +  Anb  is  the  mode  index,  and  x  is  the  coupling  distance.  The  result 
can  be  expressed  in  terms  of  the  notation  of  Tien  and  Ulrich  by  using 


2  2  1/2 

(n^-n^)  =  n^cos9^  -  n^/tan0^ 


(ri2) 


valid  for  waveguides  like  LiNbO^  having  An<<n.  In  this  equation  0^  -  7r/2  is 
the  angle  of  incidence  on  the  waveguide  surface  of  the  bounce  modes  described 
by  Tien  and  Ulrich. 

2  2  1/2  2  2  1/2 

In  Eq . (C 11)  we  have  also  used  ''^^21  /(n^-n  )  in  place  of 

the  transmission  T  of  Tien  and  Ulrich.  Additionally,  the  coefficient  y  in 

Eq,  (Cll)  is  the  expression  that  results  when  the  quantity  of  Reference  Cl 

is  expressed  in  terms  of  the  strength  parameter  h  and  refractive  index  parameters. 

This  is  done  using  Eq.  (24)  of  Reference  Cl.  Finally  in  Eq .  (Cll),  D  is  used  to 

describe  waveguide  depth  instead  of  the  W  used  by  Tien  and  Ulrich, 

Note  that  in  the  differential  Eq.  (Cll)  is  the  amplitude  of  one 

of  the  two  equivalent  bounce  waves  into  which  the  waveguide  mode  may  be 

resolved.  The  peak  field  of  the  mode  is  2A^ .  The  field  at  the  surface  is 

found  in  terms  of  the  peak  field  by  solving  the  waveguide  boundary  conditions. 

(C3) 

For  a  strongly  assymetrlc  waveguide  like  LlNbO^,  the  analysis  of  Marcuse 
indicates  that 
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E  (x)  =  2A  (x)(nj-n^)^'^^(n^-l)"’'^^  .  (C13) 

g  J-  1  8  8 

Taking  =  n  +  An  and  n^  =  n  +  Anb,  with  An<<n,  this  result  is  expressed 
approximately  as 

E  (x)  =  2A^  (x)  [  2nAn(l~b)  1  .  (C14) 

It  is  instructive  to  substitute  this  result  and  the  result  of 
I>|s,  (C9)  and  (CIO)  into  Eq .  (ClI)  to  arrive  at 

Ji:,,(x)/dx  =  (2An(l-b)/D]f(n“-l)~^^’E  -  E  (x)  ]  .  (C15) 

p  p  p  g 

i’h i s  differential  equation  relates  the  waveguide  field  at  the  surface  to  the 
evanescent  prism  field  at  the  surface.  For  the  case  of  weak  coupling,  the 
'.ocond  term  in  Eq .(CIj)  i '  ly  be  neglected,  and  the  waveguide  surface  field  is 
found  to  be 


E  (x)  =  [2Ai-i(l-b)/D](n^~l)~^'^^E  x  .  ('C16) 

P 

11  X  =  a  is  the  effective  length  over  which  input  coupling  occurs,  we  find 

E  (f)  =  [1.021An(l-b)f/D]E  ,  (C17) 

P 

whore  we  have  used  n  =  2.2  to  obtain  the  numerical  coefficient  1.021.  If 

An  =  0.005,  b  =  0.5,  D  =  3.0ym,  and  £  =  =  1mm,  the  beam  width  and  coupling 

spot  size  we  find  E  (£)  ==  0.851  E  .  The  waveguide  surface  field  is  then 
8  P 

comparable  to  the  evanescent  prism  Field.  In  accord  with  the  discussion 
associated  with  Eqs . (C6)-(( 8), the  scattering  from  the  input-coupling  region 
is  enhanced  by  the  contribution  of  the  prism  field,  but  it  remains  small  in 
comparison  to  scattering  from  the  free  waveguide,  owing  to  disparity  in  the 
path  lengths . 

This  conclusion  is  based  on  the  assumption  that  Z  “  1mm.  A  very 
different  conclusion  can  result  if  we  suppose  that  effective  coupling  occurs 
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over  a  smaller  distance.  This  could  occur,  for  example,  if  the  incident  beam 
has  a  slight  divergence,  so  that  the  correct  coupling  angle  is  achieved  only 
across  part  of  the  beam.  This  would  be  the  part  nearest  the  termination  of  the 
coupling  spot,  in  order  to  avoid  coupling  out  part  of  the  waveguide  beam. 

To  get  an  indication  of  what  value  £  might  have,  we  begin  by 
deriving  a  relation  for  coupling  efficiency.  The  defining  expression  is  (Cl) 


(c/47r)  n  D  Ai(£)A*(£) 

(c/Stt)  n  W  A^A* 
p  J  3 


(C18) 


where  the  modal  field  A^(£)  may  be  determined  from  the  field  at  the  surface 
E  (£)  or  by  integrating  the  differential  Eq.(Cll).  For  low  coupling 

O 

efficiencies  such  as  are  encountered  in  practice,  we  use  A^(£)  =  aA^£,  where 
a  is  given  in  Eq.(Cll).  We  also  note  the  relation  between  beam  width  W  and 
coupling  spot  size  L^,  ^  cos0^,  where  9^  =  sin  ^(n^/n^)  is  the  angle  of 

incidence  on  the  base  of  the  coupling  prism.  Substituting  these  results  in 
Eq.  18,  we  obtain 

n  =  I  T  |^(nJ-n^)(n^-n^)"^^V/2nL  D  .  (C19) 

'  3zl  1  g  p  c 

2  2  2  2  2 

From  the  defining  Eq.  10,  we  find  I Tooi I  “  (h) (n  -n  )/(n  -1),  where 

2  2  2  ^  ^ 
f  (h)  is  given  in  Eq.(Cl]).  f  (h)  varies  from  4h  for  small  values  of  h  to 

approximately  2  for  h=l. 

We  are  interested  in  the  possibility  of  attaining  reasonable 
coupling  efficiencies  with  small  values  of  £.  This  is  best  accomplished  if 
I  ^321 1  ^  has  its  largest  value,  l*’^32ll^  ~  ^^”p  ”  ^  )/(np-‘l)  =  3.755.  Eq.  19 
then  becomes, 

n  =  2.037An(l-b)£^/L^D  ,  (C20) 

2  2 

where  we  have  used  n^-n^  -  2nAn(l-b)  .  For  An  =  0.005,  b  *  0.5,  £  =  200vim, 

=  Imm,  and  D  »  3.0vim,  we  obtain  n  =  6.8%  for  the  input  coupling  efficiency. 
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We  typically  observe  waveguide  through-put  in  the  range  5-10%.  This  is 

marginally  consistent  with  the  above  value  for  n  if  output  coupling  efficiency 

is  in  the  range  30-100%.  This  high  efficiency  is  possible  because  there  is 

no  restriction  on  output  coupling  length  Imposed  by  beam  size  or  divergence, 

as  is  the  case  with  input  coupling. 

Using  t  *  200ym  in  Eq.(C17),  we  determine  E  (£)  »  0.170  E  . 

8  P 

That  is,  the  peak  waveguide  field  is  approximately  one-sixth  the  magnitude  of 

the  prism  field  at  the  waveguide  surface.  The  prism  field  is  not  effected 

much  by  slight  beam  divergence,  and  may  be  regarded  as  constant  over  the 

coupling  length  L  .  The  maximum  possible  scattered  intensity  from  this 
^  2 

field  is  then  proportional  to  >  which  is  to  be  compared  to  scattering 

frtnn  the  free  waveguide,  proportional  to  E^L  .  The  ratio  of  prism-coupling 
scattering  to  wavegixide  scattering  is 

R  =  I  /I  =  e\  /e’l  ,  (C21) 

eg  P  c  g  g 

A  numerical  value  R  =  2.3  is  found  for  the  example  E  =  0.170  E  ,  L  = 

g  p  c 

1,0  mm,  and  L  =  15.0  mm. 
g 

A  more  accurate  analysis  would  employ  an  effective  coupling  length 
in  Eq.(C21)  to  account  for  the  fact  that  some  of  the  light  scattered 
into  waveguide  modes  by  the  prism  field  is  out-coupled  before  it  reaches  the 
free  waveguide.  The  quantity  e  the  probability  that  a  photon 

scattered  at  x  will  reach  x  ==  without  being  lost  to  the  waveguide,  where 
Y  is  the  amplitude  attenuation  coefficient  given  in  Eq.(Cll).  The  power  scattered 
by  the  prism  field  E^  which  reaches  the  free  waveguide  is 


E^L-  = 
pc  p*'  o 


(C22) 


from  which  we  obtain 


l’  =  (l/2Y)(l-e~^^^C) 


(^2  3) 
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For  the  numerical  parameters  used  in  the  calculation  so  far,  we  find  = 

481ym.  The  light  scattered  from  the  prism  coupling  region  then  reduced  to 
111%  of  the  total  light  scattered  from  the  free  waveguide. 

While  it  must  be  acknowledged  that  the  derived  result  R  =  1.11 
was  obtained  through  the  choice  of  favorable  parameters,  the  calculation 
does  add  support  to  the  experimental  observation  that  prism-associated  scattering 
is  a  significant  part  of  the  total  in  plane  scattering  measured.  It  is  instruc¬ 
tive  to  consider  that  the  total  effective  scattering  path  length  is  given 
by 


L  =  L  +  L 
g  e 


(C24) 


where  =  RL^  is  the  effective  additional  scattering  path  length  associated 
with  prism  enhancement.  From  Eqs . (Cl 7) , (C2l) , (C23) ,  and  (Cll)  we  obtain 


L  =  0.96  [D/£An(l-b)]^(l/2Y)  (1-e 
e 

Y  =  0.51[An(l-b)/D]f^(h) 


(025) 


In  obtaining  this  result,  we  have  used  from  Eq.(C22)  in  place  of  in 
Eq.  (C21). 

If  we  perform  scattering  measurements  with  the  coupling  prisms 

separated  first  by  L  and  then  by  l’ ,  the  ratio  of  scattered  intensity  in 

g  g 

the  two  measurements  is  given  by 


I  /I' 
s  s 


(L  +  L  )/(L'  +  L  ) 
g  ^  g  e 


(C26) 


where  it  is  assumed  that  the  effective  length  of  the  coupling  region, 

is  the  same  in  the  two  experiments.  We  have  used  =  L  /lO  in  experi- 

^  g  g 

ments.  This  leads  to 


I  /I'  =  (1  +  L  /L  )/(0.1  +  L  /L  ) 

s  s  eg,  eg 


(C27) 
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If  the  ratio  is  approximately  2.  In  other  words,  a  factor  of  10 

^.'hange  in  free  waveguide  path  length  is  accompanied  by  only  a  factor  of 
2  change  in  scattering.  This  is  characteristic  of  our  observations.  We 

have  even  obtained  results  in  which  I’  >  I  ,  as  can  occur  if  L  increases 

s  '\^  s  e 

when  the  path  length  is  reduced  owing  to  variations  In  coupling  efficiency. 

The  main  assumption  required  to  obtain  values  of  L^/L^  as  large 

as  uniLy  is  that  of  a  small  effective  coupling  longtl^  P.  Supporting  evidence 

!*m-  a  small  value  of  f  is  the  fact  tluit  coupling  to  a  waveguide  mode  occurs 

v>v<‘r  a  range  of  tuning  angles  that  happens  to  be  about  1  inrad  in  our  experiments. 

This  r.inge  includes  the  natural  tuning  range  for  tl\e  mode  and  a  contribution 

\  Vim  external  factors  sucli  as  beam  divergence.  We  will  use  AO  ^  to  describe 

^  nat 

t  lu‘  natural  angular  turning  range  and  AO,.  to  describe  the  tuning  range 

u  1 V 


.,ssi*riartd  with  beam  divergence.  II 
-  ‘i’et  elticient  coupling,  to  occur  over  a  distanci 


i  .s  cimu)ared  to  A^"''  ,  wo 

d  Lv  '  nat 


(*  !.  Ai'  /AO,, 

c  na  t  d  IV 


(C28) 


A  value  for  mav  be  tvileulated  frc>m  formulas  presented  in 

((:2) 

{  lu'  paper  by  Ulrich.  It  n  ^^in  ^  a  del  i  m  s  the  ..inc,Le  (M  incidenci^  - 

•  ^  ^  P  ^  g  3 

(er  coupling,  to  tlie  nu^de  havinc,  tlit*  td!icli\'e  indi^y  n^,  a  slight  variation 

ill  sav  .^0,  causes  a  change  in  n  bv  an  amv>unt 

3  ■  g  • 


6n  =  n  cos  0..  60 
g  P  3 


.2  ’1/2 
(n^  -  a  )  -0 


(C29) 


Ulrich  shows  that  the  energy  density  of  light  in  the  waveguide  region  falls 
off  in  boron tzian  fashion  with  Increasing,  6n  ,  or  equivalently,  with  increasing 
60.  'fhe  ful  l-width-at-half-maximiim  ener^g/  Lk-nsitv  is 


2  2  1/2 

An  =  (n  -n  )  AO 
g  P  ■  nat 


2K 


(C30) 


(C2) 


in  Ulrich’s  notation,  where  K  is  given  bv  Kq .  (37)  of  his  paper. 

m 

Transferring  our  notation  to  his  result  we  have,  in  the  limit  h<<l. 
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A9  «  16h^[2n  An(l~b) (n^-1) 
nat  p 

X  *  *0  b)  /An 

\l>  =  2k^D[2nAn(l-b)]^^^  -  2(Ji^q  -  24,^^^  (C31) 

=  tan"^[b/(l-b)]^^^ 

((>^2  “  tan  ^  [  (n^-1)  /  (2nAn)  (1-b)  ] 

In  evaluating  we  may  neglect  the  contribution  since 

^12  ~  values  of  b  that  are  of  interest.  We  find 

■;/.’ib  =  -[b(l-b)]^'^^[k^D(2nAnb)^^^+l]  .  (C32) 


By  inserting  this  result  into  Eq,  (C27),  we  obtain 

AG^at  "  [16h^/(np  -  1) ] [ An (1-b) /k^D] [Q/ (Q+1 ) ] 
Q  =  k^D(2nAnb)^'^^ 


(C33) 


For  the  parameter  values  used  in  this  section,  including  h  =  1,  An  =  0.005, 

b  =  0.5,  D  =  3.0pm  and  k  =  2TT/0.633ym,  A9  ^  is  in  the  vicinity  of  0.1  mrad, 

o  na  L 

This  is  about  one-tenth  the  total  tuning  range  that  we  observe  for  our  waveguides, 
suggesting 


A0  /A6  ,  .  =  0.1  .  (C34) 

nat  div 

In  this  case,  only  one-tenth  of  the  extent  of  the  incident  wavefront  is  capable 
of  coupling  efficiently  to  the  waveguide,  or  £  =  L^/10,  This  is  comparable 
to  the  values  found  in  our  earlier  calculations  to  give  large  prism  enhancement 
of  scattering. 
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One  possible  alternative  to  the  euncept  of  prism-enlianced  waveguide 
scattering  holds  that  the  scattering  results  from  inperfecti^ ns  in  the  coupling 
prisms.  We  have  several  observations  that  suggest  that  this  is  not  the  case. 
First,  the  angular  distribution  of  inplane  scattering  does  not  change  much 
when  the  prism  separation  is  reduced.  This  suggests  that  the  size  of  the 
scattering  centers  is  the  same  in  the  waveguide  region  as  in  the  prism  region. 
Barring  coincidence,  this  suggests  that  tin,*  scattorers  are  Localized  in  Lite 
waveguide.  Second,  we  have  observed  in  our  best  waveguide  a  very  low  level 
of  scattering,  even  though  no  special  precautions  were  taken  with  the  prism 
couplers.  This  wavt»guide  was  discussi*d  in  d(*tai]  in  Uei  .  1.,  and  its  sc  at¬ 
tering  performance  is  reproduced  in  Figure  C2.  We  interpret  the  result  as 
indicating  that  the  reduction  of  waveguide  scattering  centers,  as  was  done  by 
polsihing  in  this  case,  is  accompanied  l)y  the  reduction  of  prism-enlianced 
scattering  as  well.  This  is  consistent  with  our  model,  but  not  consistent 
with  the  idea  that  the  prisms  have  fixed  stcittcring  centers  associated  with 
them. 

Finally,  we  observe  that  pri  s!n^ct;up  li  ng  seat  tiring  is  greater  in  the 
case  of  input  coupling  than  output  I'oupjJng.  Tliis  would  not  be  the  case  for 
scattering  centers  fixed  to  the  prism,  hut  it  is  tlie  c<ise  lor  prism  onltanced 
waveguide  scattering.  The  mechanism  that  our  model  employs  for  pri sm-enhanred 
scattering  requires  the  evanescent  field  of  the  input  beam  to  act  as  the  source 
of  scattering.  This  source  is  not  aval  1  able  in  output  coup! ing .  Hence,  less 
enhancement  is  predicted  in  the  case  of  output  coupling. 

Wc  observed  that  this  prediction  was  satisfied  by  an  experiment  in 
wl\ich  end- fire  coupling  was  used  along  with  a  single-prism  coupler.  Total 
scattering  was  less  when  the  prism  was  used  for  output  coupling  and  input 
coupling  was  accomplished  by  end  fire. 

We  conclude  that  tlie  model  tjf  pr  i  sm-en  ha  need  waveguide  scattering 
exhibits  many  of  the  features  that  wt‘  have  encountered  in  experiments.  Let  us 
consider,  then,  how  this  effect  may  bo  red uc oil  in  future  observations.  The 
main  contribution  leading  to  a  large  prism  eidianc'ement  is  the  existance  of  a 
small  effective  length  (  =  lOOpm  over  whic'h  input  coupling  can  occur.  If  f 
he  increased  ten  fold  to  the  value  associated  with  the  coupling  spot  size  or 
beam  width,  prism  enhancement  of  the  scattering  is  calculated  to  be  negligible. 
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ANGLE 


In-plane  scattered-energy  distribution 
for  the  TEq  mode  propagating  parallel 
to  the  optic  axis,  after  11  h  total 
polishing  time. 


We  determined  above  that  I  =  JOOum  was  compatible  with  a  beam 
divergence  of  about  1  mrad.  This  is  the  natural  divergence  angle  of  many 
low  power  He-Ne  lasers,  consistent  with  tiie  typical  beam  width  W  =  lirnn. 
Increasing  t  should  simply  be  a  matter  of  carefully  controlling  beam  collination 
and  insuring  that  the  Gaussian  waist  coincide  with  the  input  coupling  region* 
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APPENDIX  D 


DETERMINATION  OF  WAVEGUIDE  PARAMETERS  FROM  DIFFUSION  CONDITIONS 

It  may  be  noted  that  the  determination  of  scattering  cross  sections 
using  the  various  formulas  of  the  text  Inevitably  requires  knowledge  of  the 
waveguide  index  parameters  An  and  D,  These  are  generally  determined  using 
the  relations 


An  =  AT/D 

D  = 


(Dl) 


where  A  is  an  empirically  determined  coefficient,  T  is  the  Ti-film  thickness, 

T  is  the  diffusion  temperature,  P(T)  is  the  diffusion  coefficient  at  that 
temperature,  and  t^  is  the  diffusion  time. 

The  purpose  of  this  appendix  is  to  point  out  the  caution  that  must 
be  taken  in  employing  this  type  of  analysis.  The  largest  source  of  error  is 
in  the  film  thickness  T.  This  thickness  is  intended  to  correspond  to  that  of 
an  oxygen-free  Ti  film  prior  to  diffusion.  However,  Ti  is  an  effective  getter 
of  oxygen  and  our  e-beam-cvaporated  films  invariably  show  more  transparency 
than  they  would  be  expected  to  for  a  pure  Ti  film,  owing  to  the  presence  of 
residual  oxygen  in  the  vacuum  system.  Measured  values  of  T  are  therefore  too 
large  by  a  factor  that  depends  on  the  oxygen  content  of  the  film.  Our  own 
experience  is  that  this  content  depends  strongly  on  the  vacuum  system  employed. 
Therefore,  in  applying  the  empirical  relationship  An  =  AT/D,  the  waveguide  in 
question  should  be  fabricated  from  the  same  vacuum  system  used  in  experiments 
to  determine  the  constant  A. 

In  our  recent  work,  LiNbO^  waveguides  have  been  diffused  from  e-beam 
evaporated  films  formed  in  two  different  vacuum  systems.  In  each  case,  the  e-beam 
gun  was  the  same.  The  sample  of  Figure  20,  emphasized  in  this  report  because 
of  its  examination  using  end-fire  coupling,  was  formed  using  the  more  efficient 

o 

vacuum  system.  It  was  diffused  from  a  Ti  film  of  thickness  250  A  as  determined 
by  a  crystal  thickness  monitor  and  315  H  as  determined  using  an  angstrometer . 


i)-2 

Optical  transmission  was  estimated  to  be  in  the  range  10-25%.  Most  of  the 
other  waveguides  of  the  program  were  diffused  using  slightly  more  trans¬ 
missive  films,  say  33%  transmissive.  These  films  were  evaporated  in  a  less 
efficient  vacuum  system.  For  a  given  thickness  these  films  would  have  had 
less  Ti  and  therefore  would  produce  waveguides  with  a  smaller  /In.  We  also 
note  that  the  waveguides  having  a  greater  oxygen  content  entered  the  substrate 
more  rapidly.  This  would  indicate  a  larger  diffusion  depth  1).  Since  most 
scattering  mechanism  decrease  in  strength  with  smaller  An  and  larger  D,  it 
appears  appropriate  to  diffuse  films  that  luive  been  evaporated  in  relatively 
impertect  vacua.  It  also  appears  reasonable  to  carry  out  the  diffusion  in  a 
oxygen- rich  atmosphere,  as  we  do,  ratlier  than  in  argon. 

It  is  ironic  that  the  waveguide  of  Figure  20,  which  shows  the  best 
scattering  performance  of  any  waveguide  fabricated  during  the  current  program, 
was  fabricated  under  non-optimum  conditions  leading  to  a  large  An  and  small 
0.  Wc  believe  that  the  good  performance  resulted  from  the  elimination  of 
prism-enhanced  scattering  and  that  much  better  performance  may  be  expected 
from  waveguides  fabricated  under  the  more  optimum  conditions  outlined  in 
Sec.  V. 


APPENDIX  E 


MEASUREMENT  OF  SURFACE  ROUGHNESS  FROM 
SCATTERING  OF  REFLECTED  LIGHT 


This  appendix  describes  a  series  of  measurements  that  we  made  in  order 
to  determine  the  surface -roughness  parameters  of  our  diffused  LiNbO^  waveguides. 
These  measurements  were  made  with  the  idea  of  correlating  surface  roughness  with 
measured  in-plane  scattering  levels,  using  the  theoretical  results  of  Sec.  IV. 
The  measurements  proved  to  be  exceedingly  complicated  and  of  limited  precision, 
so  this  plan  was  not  fully  implemented.  However,  those  measurements  that  were 
made  are  worth  discussing. 


EXPERIMENTAL  CONFIGURATION 

The  experimental  configuration  that  was  employed  is  shown  in  Fig.  El. 

A  gently  focused  laser  beam  was  incident  on  the  sample  at  approximately  45®.  It 
is  critical  that  the  beam  be  highly  filtered  spatially  so  that  the  only  scatter¬ 
ed  light  observed  in  the  far  field  is  that  associated  with  the  sample  surface. 
This  filtering  is  most  readily  accomplished  if  the  beam  is  focused  on  or  near 
the  sample.  However,  the  size  of  the  beam  incident  on  the  sample  should  be 
large  in  comparison  to  the  roughness  autocorrelation  length  in  order  to  justify 
a  statistical  analysis  of  the  results. 

Note  from  Fig.  El  that  scattered  light  is  monitored  in  the  actual  far 
field.  A  lens  is  not  used  to  generate  the  far-field  condition  because  of  the 
scattering  which  would  be  introduced  by  surface  roughness  and  imperfections  in 
the  lens.  Another  problem  is  that  of  second-surface  reflections  from  the  sample. 
These  were  either  eliminated  by  blocking  with  a  razor  edge,  or  they  were  allowed 
to  contribute  to  the  scattered  signal.  Blocking  is  difficult  because  light  dif¬ 
fracted  from  the  razor  is  often  comparable  to  or  larger  than  the  scattered 
signal.  On  the  other  hand,  allowing  the  second  surface  to  contribute  to  the 
scattering  introduces  an  unknown  quantity  into  the  measurement  because  the 
second  surface  will  not  contain  a  diffused  waveguide,  generally,  and  may  also 
have  a  greater  incidence  of  scratches  and  imperfections  introduced  during  sample 
handling. 
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Add  to  this  the  fact  that  the  surfaces  of  LiNbO^  are  intrinsically 
very  smooth  (see  Appendix  A)  and  one  can  begin  to  appreciate  the  difficulty  of 
the  measurement.  We  will  be  describing  the  measurement  of  scattering  levels  in 
the  range  -60  to  -70  dB. 


THEORETICAL  ANALYSIS 


Our  analysis  of  surface  scattering  parallels  that  of  Sec.  IV  for 
Rayleight-Gans-Debye  scattering  from  waveguide  index  inhomogeneities.  We  con¬ 
sider  in  Fig.  E2  a  single  roughness  element  in  an  otherwise  smooth  surface.  The 
roughness  height  is  a  and  the  lateral  size  (autocorrelation  length)  is  Cl,  Light 
reflected  from  the  surface  has  its  wavefront  locally  perturbed  as  a  result  of 

the  phase  perturbation.  The  strength  of  the  perturbation  is  Acj)  =  k  o  cos0,  and  it 

o 

occupies  a  region  of  size  acosO, in  the  near  field. 

In  the  far  field,  the  amplitude  distribution  is  that  of  the  unperturbed 
beam  plus  the  contribution  from  the  roughness  element,  as  described  by  the  laws 
of  diffraction.  For  the  perturbed  wavefront  in  the  near  zone  has  the  peak 

amplitude 


A  =  iA  k  ocosG 
no  o  0 


(FI) 


If  the  roughness  element  is  circular  in  the  plane  of  the  surface  with  radius  a, 
the  wavefront  perturbation  will  be  elliptical,  with  major  axis  a  perpendicular 
to  the  plane  of  incidence  and  minor  axis  acosO  in  the  plane  of  incidence. 

For  convenience  we  will  assume  a  Gaussian  profile 


A(x,y) 


A  +  y^) 

no 


(E2) 


where  x  and  y  are  coordinates  attached  to  the  reflected  beam:  the  z  axis  is  the 
direction  of  propagation,  the  x  axis  is  in  the  plane  of  incidence,  and  the  y  axis 
is  normal  to  the  plane  of  incidence. 

In  the  Fraunhofer  zone,  the  amplitude  distribution  is  the  Fourier- 
transform  of  Eq.  (E2) ,  or 


A,(X,Y)  =  -dk  /2T,z)e^‘'o"  dxdyA^(x.y)e 
t  o 


i(k  /z)(xX  +  VY) 

D 


(E3) 
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Fig.  E2.  Geometry  for  analysis  of  scattering  caused  by 
reflection  from  a  rough  surface. 


where  X,Y  are  the  coordinates  in  the  observation  plane  that  correspond  to  x,y 
in  the  near  zone,  and  z  is  the  distance  of  the  observation  plane  from  the  plane 
of  the  sample.  The  result  of  carrying  out  the  integration  is 


A,(X,Y)  =  -A  e^‘^o^(k^oa^cos^0/8z)expf-x\^a^cos^e/16z 
too  o 


2.  2.2  2, 


Hbz^] 

o 


The  power  from  this  scattered  field  intercepted  by  a  detector  of  area  S  at  a 
distance  z  from  the  sample  is 

6P  =  (c/8ti)  IAjI^S  .  (E5) 

As  for  an  in-plane  scattering  measurement,  this  power  is  referenced  to  the  peak 
power  in  the  unscattered  beam.  This  power  is  determined  by  a  Fourier-transform 
similar  to  Eq,  (E3)  for  the  incident  field 


A(r)  =  A  e 
o 


-4r^/W^ 


where  W  is  the  width  of  the  incident  beam  on  the  surface,  measured  between 
2 

opposite  1/e  intensity  points.  By  inserting  Eq ,  (E6)  in  place  of  A^(x,y)  in 
Eq.  (E3)  we  obtain 

A  ,'D%  -A  /I  /o  ^  ikR  /,,2,,,  -(R^/z^)kV/16 

Au(R)  =  xA^(k^/27Tz)e  /4)e  o  (E7] 

as  the  unscattered  amplitude  in  the  far  zone.  Note  that  r  and  R  are  the  radial 
coordinates  corresponding  to  x,y  and  X,Y,  respectively. 

The  peak  power  collected  by  the  detector  of  area  S  is 


P  =  (c/Stt)  |A^(0)  I  S.  (E8) 

The  relative  scattered  power  is  obtained  by  dividing  Eq .  (E8)  into  Eq .  (E5)  and 
substituting  the  results  of  Eqs.  (E7)  and  (E4) .  The  result  is 
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[x  lx  'I 

(6P/P)  =  (k  a  cos  9/W  )  exp(-((5  cos^0/16 

O  X  o 

2  2  2 


(E9) 


where  ((>^  =  X/z  and  =  Y/z  are  Introduced  to  indicate  the  angle  of  scattering. 

Thus  far  we  h^ive  considered  only  scattering  from  a  single  roughness 
element  positioned  at  the  most  intense  position  of  the  beam.  Other  elements  see 
less  incident  power  and  their  scattering  strength  is  reduced  proportionately. 

Lot  x*,y*  be  cartesian  coordinates  of  the  plane  of  the  surface,  with  the  origin 
at  beam  center  and  the  x’  axis  in  the  plane  of  incidence.  Owing  to  the  oblique 
approach  of  the  beam,  the  intensity  at  x*,y’  is 


I(x'.y')  =  1(0)  exp[-8x^cos‘^e/W  -  8y'^/W]  . 


(ElO) 


Within  an  increment  of  area  dx’dy’  the  total  number  of  roughness  elements  is 
2 

dx’dyV^X  ,  and  their  contribution  to  scattered  power  is  6P  times  the  exponential 
term  in  Eq .  (10).  The  total  scattered  power  is 


AP  =  (SPa  ^  JJ  dx’dy’  e 


,  ,  ,  -8x' “cos‘e/W^  -Ky'^/W*^ 


(Ell) 


=  6P(TTW"/8a^cos0)  . 

The  total  relative  scattered  power  is  obtained  from  Eq .  (E9)  by  multiplying  by 
(AP/6P).  The  result  is 

AP/P  =  (tt/8)  (k^a^Cl^cos^6/W^)  exp(-(t)^k^a^cos^6/16  (E12) 

o  X  o 

-4.\^a^/16)  . 

y  o 

Experimentally,  we  scan  along  the  X  axis,  that  is,  in  the  plane  of  in¬ 
cidence.  Thus  (|)  =  0  in  Eq .  E12,  and  we  replace  (j)  by  <(>  for  simplicity.  The 

y  ^ 

equation  to  use  in  analyzing  the  data  is 

2  2  2 

/D  t  /Q\n2  2.2  3„,,,2,.  -ij)  k  a  cos‘'0/16  .  /■im\ 

AP(<}i)/P  =  ('Tt/8)(k  0  d  cos  0/W  )  e  o  (.EIJ) 

o 
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Note  that  W  is  the  waist  of  the  incident  beam,  since  we  assumed  in  the  analysis 
that  the  incident  wavefronts  were  planar  in  the  vicinity  of  the  sample.  However, 
the  experiment  may  be  performed  with  a  diverging  or  converging  beam.  The  scat¬ 
tered  power  will  remain  constant  since  the  number  of  scattering  centers  goes  up 
to  compensate  for  the  lower  intensity  at  the  sample  surface.  One  must  take 
care  to  use  the  beam  waist  for  W  in  Eq.  E14  rather  than  the  beam  diameter  at  the 
sample - 


RESULTS  AND  ANALYSIS 


Figure  E3  shows  the  first  experimental  results,  obtained  for  sample 
154,  In  plotting  the  normalized  scattered  intensity  versus  scattering-angle 
squared  we  expect  to  find  a  linear  variation  except  near  the  orgin  where  the 
unscattered  intensity  is  large.  These  data  are  not  plotted  in  Fig.  E3,  so  a 
good  straight  line  is  obtained  down  to  the  noise  level  of  the  system.  The  slope 
and  intercept  of  this  line  provide  information  regarding  the  autocorrelation 
length  a  and  the  surface  roughness  a.  The  relevant  formulas  obtained  using 
Eq .  E13  are 

slope  -  (10  log  e) (k  acos0/4)^  (E15) 

o 

intercept  =  10  log  f  (tt/S)  (k^a^a^cos^e/W")  ] 

o 

The  first  expression  is  used  to  find  a  value  for  CL,  This  value  is  then  inserted 
into  the  second  expression  to  enable  the  determination  of  a  value  for  o.  Even 
if  the  data  shows  a  good  straight-line  behavior,  as  does  Fig.  E3,  we  caution  the 
reader  not  to  place  significance  on  the  values  so  obtained  beyond  a  factor  of 
about  2.  Our  derivation  has  not  followed  a  mathematically  rigorous  path  (K2)  . 
and  even  if  it  had  we  would  ultimately  have  to  make  assumptions  regarding  the 
surface  roughness  statistics  valid  for  our  sample.  Thus  we  take  some  liberty 
when  we  refer  to  (X  as  the  autocorrelation  length  and  o  as  the  rms  surface  rough- 

o 

ness.  Nevertheless  the  data  of  Fig,  E3  indicates  that  ^  52  pm  and  o  =  29  A. 

To  obtain  these  values,  we  used  k^  =  (2ti/0 .633)  pm  ^0  =  47®,  and  W  =  142  pm. 

It  was  this  experiment  that  led  us  to  believe  that  surface  roughness 

<• 

was  an  Important  source  of  the  observed  scattering  level  In  sample  154,  as  shown 
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in  Fig.  20.  Although  we  later  found  that  surface-compound  formation  could  also 
explain  our  in-plane  scattering  results,  we  embarked  on  a  program  to  measure 
surface  roughness  for  various  waveguide  samples  and  to  correlate  the  results  with 
in-plane  scattering  levels  for  the  same  waveguides.  This  led  to  refinements  in 
the  Experiment  of  Fig.  El  and  made  us  aware  of  the  experimental  difficulties  of 
the  measurement  which  we  had  not  appreciated  when  the  data  of  Fig.  E3  were  taken. 

First  of  all,  we  had  assumed  that  all  scattering  observed  was  that  from 
roughness  in  the  LiNbO^  surface .  In  fact,  it  is  impossib le  to  filter  out  all 
background  scattering  from  optics  used  in  the  measurement,  so  an  initial  deter¬ 
mination  of  this  background  scattering  must  be  undertaken.  This  was  done  by  re¬ 
moving  the  sample  from  the  optical  path  and  measuring  the  scattered  energy  in 
the  incident  beam^  Then  a  similar  measurement  was  made  with  the  crystal  in  place. 

These  measurements  were  made  for  LiNhO^  substrate  165.  The  NBS  Taly- 
surf  examination  of  the  similar  substrate  164  had  shown  the  surface  to  have  a 

o 

smoothness  less  than  or  comparable  to  3  A  .  One  does  not  therefore  expect  to 
see  much  additional  scatter  produced  by  the  LiNbO,^  surface.  The  experimental 
results  of  Fig,  E4  show  this  to  be  true.  Data  taken  with  and  without  the  crystal 
in  place  are  identical  to  within  experimental  uncertainty.  Tiie  error  bars  in 
Fig.  E4  pertain  to  the  data  taken  without  tlie  crystal  in  place,  shown  as  solid 
data.  Several  runs  were  made  and  the  variations  i  {i  scattered  intensity  nt  a 
given  angle  were  used  to  determine  the  error  bars.  fhe  positi(>n  of  the  dots  are 
determined  by  averaging  the  results  for  two  runs. 

Best-fitting  straight  lines  for  the  data  are  also  plotted  in  Fig.  E4, 

The  dashed  line  pertains  to  scattering  observed  with  tlie  sample  in  place.  The 
solid  line  refers  to  scattering  from  the  incident  beam.  Tlie  proximitv  of  cho  two 
lines  in  comparison  to  the  size  of  the  error  bars  sliows  that  the  LiNbO..  roughness 

o  ^ 

cannot  be  measured  by  this  technique.  In  fact,  if  a  -  3  A  as  determined  by  NBS 
for  a  similar  sample,  we  expect  the  intercept  for  surface  roughness  scattering  to 
fall  between  -70  and  -75  dB,  almost  three  orders  of  magnitude  below  the  measured 
background  scattering  level. 

Figure  E5  shows  the  results  of  a  comparable  experiment  performed  after 

a  Ti-dif fused  waveguide  was  fabricated  in  sample  165.  The  fabrication  conditions 

0 

were  similar  to  those  for  sample  169;  that  is,  a  600  A  Ti  film  was  diffused  for 
8  h  at  1000°C.  The  surface  roughness  prod\ired  by  the  diffusion  was  sufficient  to 
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cause  a  measurable  scattered  signal  above  the  background  level.  The  scattering 
associated  with  the  waveguide  alone  was  determined  by  subtracting  the  value  of 
AP/P  associated  with  the  background  from  the  value  associated  with  the  back¬ 
ground  plus  waveguide.  When  this  was  done  as  a  function  of  scattering  angle, 
using  the  straight  lines  in  Fig.  E5  to  determine  AP/P  values,  a  curve  was  con¬ 
structed  which  is  virtually  a  straight  line  having  intercept  -45.67  dB  and 
4  2 

slope  6.45  X  10  dB/radian  .  This  straight  line  is  plotted  in  Fig.  E6  along 
with  the  straight  line  obtained  from  Fig.  F3  for  sample  154.  Despite  signifi¬ 
cant  differences  in  slope  and  intercept,  the  two  lines  are  found  to  be  consis¬ 
tent  with  fai rly  similar  values  for  d  and  o .  We  find 

a (154)  =  52  pm 
a(165)  =  72  pm 

o 

0(154)  =  29  A 
0(165)  =  29  A 

This  is  surprising  in  view  of  the  very  different  fabrication  conditions  rv>r  the 
two  waveguides.  However  it  is  consistent  with  (uir  (H>servation  tliat  waveguide 
169,  fabricated  similarly  to  165,  had  tn-plano  scattering  characteristics  ('ompa- 
table  to  134,  It  may  be  that  diffusion  induced  roughness  is  smoothed  away  for 
sufficiently  long  diffusion  times. 
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